DEMO_febio_0034_sphere_cone_slide_body_force
Below is a demonstration for:
- Building geometry for a spherical blob with tetrahedral elements which is being aspirated into a tube. This demo consists off:
- Defining the boundary conditions
- Coding the febio structure
- Running the model
- Importing and visualizing the displacement results
Contents
Keywords
- febio_spec version 2.5
- febio, FEBio
- indentation
- contact, sliding, friction
- rigid body constraints
- tetrahedral elements, tet4
- triangular elements, tri3
- shell elements
- sphere
- static, solid
- hyperelastic, Ogden
- displacement logfile
- stress logfile
clear; close all; clc;
Plot settings
fontSize=15; faceAlpha1=0.8; faceAlpha2=0.3; markerSize=40; lineWidth=3; cMap=[1 0.5 0.4; 0.9 0.3 0.27; 0.8 0.2 0.18; 0.7 0.1 0.09; 0.6 0 0; 0.5 0 0; 0.4 0 0;]; [cMap]=resampleColormap(cMap,250);
Control parameters
% Path names defaultFolder = fileparts(fileparts(mfilename('fullpath'))); savePath=fullfile(defaultFolder,'data','temp'); % Defining file names febioFebFileNamePart='tempModel'; febioFebFileName=fullfile(savePath,[febioFebFileNamePart,'.feb']); %FEB file name febioLogFileName=fullfile(savePath,[febioFebFileNamePart,'.txt']); %FEBio log file name febioLogFileName_disp=[febioFebFileNamePart,'_disp_out.txt']; %Log file name for exporting displacement febioLogFileName_strainEnergy=[febioFebFileNamePart,'_energy_out.txt']; %Log file name for exporting strain energy density % Sphere parameters sphereRadius=3;% numElementsMantel=6; % Ground plate parameters tubeRadius=sphereRadius.*[1 0.1]; tubeAngle=3*(pi/180); tubeLength=abs(diff(tubeRadius))/tan(tubeAngle); % Material parameter set c1=1e-3; %Shear-modulus-like parameter MPa m1=2; %Material parameter setting degree of non-linearity k_factor=10; %Bulk modulus factor k=c1*k_factor; %Bulk modulus d=1e-9; %Density % FEA control settings max_refs=25; %Max reforms max_ups=0; %Set to zero to use full-Newton iterations opt_iter=15; %Optimum number of iterations max_retries=25; %Maximum number of retires symmetric_stiffness=0; min_residual=1e-20; analysisType1='dynamic'; timeTotal1=1; %Analysis time numTimeSteps1=100; %Number of time steps desired step_size1=timeTotal1/numTimeSteps1; dtmin1=(timeTotal1/numTimeSteps1)/100; %Minimum time step size dtmax1=timeTotal1/10; %Maximum time step size analysisType2='dynamic'; timeTotal2=100; %Analysis time numTimeSteps2=10; %Number of time steps desired step_size2=timeTotal2/numTimeSteps2; dtmin2=(timeTotal2/numTimeSteps2)/100; %Minimum time step size dtmax2=timeTotal2/2; %Maximum time step size timeTotal=timeTotal1+timeTotal2; %Contact parameters contactPenalty=20; laugon=0; minaug=1; maxaug=10; fric_coeff=0.1; %Specifying load sphereVolume=4/3*(pi*sphereRadius^3); %Sphere Volume in mm^3 sphereMass=sphereVolume.*d; %Sphere mass in tone sphereSectionArea=pi*sphereRadius^2; bodyLoadMagnitude=(9.81*1000)*75; forceBodyLoad=sphereMass.*bodyLoadMagnitude; stressBodyLoad=forceBodyLoad/sphereSectionArea;
Creating model geometry and mesh
%Control settings cPar.sphereRadius=sphereRadius; cPar.coreRadius=cPar.sphereRadius/2; cPar.numElementsMantel=numElementsMantel; cPar.numElementsCore=round(numElementsMantel*1.5); cPar.outputStructType=2; cPar.makeHollow=0; cPar.cParSmooth.n=25; %Creating sphere [meshOutput]=hexMeshSphere(cPar); % Access model element and patch data Fb_blob=meshOutput.facesBoundary; Cb_blob=meshOutput.boundaryMarker; V_blob=meshOutput.nodes; E_blob=meshOutput.elements;
Visualize blob mesh
hFig=cFigure; subplot(1,2,1); hold on; gpatch(Fb_blob,V_blob,Cb_blob,'k',0.8); patchNormPlot(Fb_blob,V_blob); axisGeom(gca,fontSize); colormap(gjet); icolorbar; camlight headlight; hs=subplot(1,2,2); hold on; title('Cut view of solid mesh','FontSize',fontSize); optionStruct.hFig=[hFig hs]; gpatch(Fb_blob,V_blob,'kw','none',0.25); meshView(meshOutput,optionStruct); axisGeom(gca,fontSize); drawnow;
Creating tube model
pointSpacingBlob=mean(patchEdgeLengths(Fb_blob,V_blob)); pointSpacingTube=pointSpacingBlob/2; rEnd=sphereRadius+(sphereRadius.*((sphereRadius-tubeRadius(2))/tubeLength)); V_curve_tube=[sphereRadius rEnd 0; -tubeLength tubeRadius(2) 0;]; nResample=ceil(max(pathLength(V_curve_tube))./pointSpacingTube); V_curve_tube=evenlySampleCurve(V_curve_tube,nResample,'pchip',0); cPar.closeLoopOpt=1; cPar.numSteps=[]; %If empty the number of steps is derived from point spacing of input curve cPar.w=[1 0 0]; [F_tube,V_tube]=polyRevolve(V_curve_tube,cPar); center_of_mass_tube=mean(V_tube,1);
Join model node sets
V=[V_blob; V_tube; ]; F_tube=F_tube+size(V_blob,1);
Visualizing model
cFigure; hold on; gtitle('Model components',fontSize); hl(1)=gpatch(Fb_blob,V,'rw','k',0.8); hl(2)=gpatch(F_tube,V,'kw','k',0.5); legend(hl,{'Blob','Tube'}); clear hl; axisGeom(gca,fontSize); camlight headlight; drawnow;
Get contact surfaces
F_contact_blob=Fb_blob;
Visualize contact surfaces
cFigure; hold on; title('Tube blob contact pair','fontsize',fontSize); hl(1)=gpatch(F_tube,V,'rw','k',0.8); patchNormPlot(F_tube,V); hl(2)=gpatch(F_contact_blob,V,'kw','k',0.5); patchNormPlot(F_contact_blob,V); legend(hl,{'Master','Slave'}); clear hl; axisGeom(gca,fontSize); camlight headlight; drawnow;
Defining the FEBio input structure
See also febioStructTemplate and febioStruct2xml and the FEBio user manual.
%Get a template with default settings [febio_spec]=febioStructTemplate; %febio_spec version febio_spec.ATTR.version='2.5'; %Module section febio_spec.Module.ATTR.type='solid'; %Control section stepStruct1.Control.analysis.ATTR.type=analysisType1; stepStruct1.Control.time_steps=numTimeSteps1; stepStruct1.Control.step_size=step_size1; stepStruct1.Control.time_stepper.dtmin=dtmin1; stepStruct1.Control.time_stepper.dtmax=dtmax1; stepStruct1.Control.time_stepper.max_retries=max_retries; stepStruct1.Control.time_stepper.opt_iter=opt_iter; stepStruct1.Control.max_refs=max_refs; stepStruct1.Control.max_ups=max_ups; stepStruct1.Control.symmetric_stiffness=symmetric_stiffness; stepStruct1.Control.min_residual=min_residual; stepStruct2.Control.analysis.ATTR.type=analysisType2; stepStruct2.Control.time_steps=numTimeSteps2; stepStruct2.Control.step_size=step_size2; stepStruct2.Control.time_stepper.dtmin=dtmin2; stepStruct2.Control.time_stepper.dtmax=dtmax2; stepStruct2.Control.time_stepper.max_retries=max_retries; stepStruct2.Control.time_stepper.opt_iter=opt_iter; stepStruct2.Control.max_refs=max_refs; stepStruct2.Control.max_ups=max_ups; stepStruct2.Control.symmetric_stiffness=symmetric_stiffness; stepStruct2.Control.min_residual=min_residual; %Add template based default settings to proposed control section [stepStruct1.Control]=structComplete(stepStruct1.Control,febio_spec.Control,1); %Complement provided with default if missing [stepStruct2.Control]=structComplete(stepStruct2.Control,febio_spec.Control,1); %Complement provided with default if missing %Remove control field (part of template) since step specific control sections are used febio_spec=rmfield(febio_spec,'Control'); febio_spec.Step{1}.Control=stepStruct1.Control; febio_spec.Step{1}.ATTR.id=1; febio_spec.Step{2}.Control=stepStruct2.Control; febio_spec.Step{2}.ATTR.id=2; %Material section febio_spec.Material.material{1}.ATTR.type='Ogden unconstrained'; febio_spec.Material.material{1}.ATTR.id=1; febio_spec.Material.material{1}.c1=c1; febio_spec.Material.material{1}.m1=m1; febio_spec.Material.material{1}.c2=c1; febio_spec.Material.material{1}.m2=-m1; febio_spec.Material.material{1}.cp=k; febio_spec.Material.material{1}.density=d; febio_spec.Material.material{2}.ATTR.type='rigid body'; febio_spec.Material.material{2}.ATTR.id=2; febio_spec.Material.material{2}.density=1; febio_spec.Material.material{2}.center_of_mass=center_of_mass_tube; %Geometry section % -> Nodes febio_spec.Geometry.Nodes{1}.ATTR.name='nodeSet_all'; %The node set name febio_spec.Geometry.Nodes{1}.node.ATTR.id=(1:size(V,1))'; %The node id's febio_spec.Geometry.Nodes{1}.node.VAL=V; %The nodel coordinates % -> Elements febio_spec.Geometry.Elements{1}.ATTR.type='hex8'; %Element type of this set febio_spec.Geometry.Elements{1}.ATTR.mat=1; %material index for this set febio_spec.Geometry.Elements{1}.ATTR.name='Blob'; %Name of the element set febio_spec.Geometry.Elements{1}.elem.ATTR.id=(1:1:size(E_blob,1))'; %Element id's febio_spec.Geometry.Elements{1}.elem.VAL=E_blob; febio_spec.Geometry.Elements{2}.ATTR.type='quad4'; %Element type of this set febio_spec.Geometry.Elements{2}.ATTR.mat=2; %material index for this set febio_spec.Geometry.Elements{2}.ATTR.name='Tube'; %Name of the element set febio_spec.Geometry.Elements{2}.elem.ATTR.id=size(E_blob,1)+(1:1:size(F_tube,1))'; %Element id's febio_spec.Geometry.Elements{2}.elem.VAL=F_tube; % % -> Surfaces febio_spec.Geometry.Surface{1}.ATTR.name='contact_master1'; febio_spec.Geometry.Surface{1}.quad4.ATTR.lid=(1:1:size(F_tube,1))'; febio_spec.Geometry.Surface{1}.quad4.VAL=F_tube; febio_spec.Geometry.Surface{2}.ATTR.name='contact_slave1'; febio_spec.Geometry.Surface{2}.quad4.ATTR.lid=(1:1:size(F_contact_blob,1))'; febio_spec.Geometry.Surface{2}.quad4.VAL=F_contact_blob; % -> Surface pairs febio_spec.Geometry.SurfacePair{1}.ATTR.name='Contact1_tube_blob'; febio_spec.Geometry.SurfacePair{1}.master.ATTR.surface=febio_spec.Geometry.Surface{1}.ATTR.name; febio_spec.Geometry.SurfacePair{1}.slave.ATTR.surface=febio_spec.Geometry.Surface{2}.ATTR.name; %Boundary condition section % -> Fix boundary conditions % febio_spec.Boundary.fix{1}.ATTR.bc='x'; % febio_spec.Boundary.fix{1}.ATTR.node_set=febio_spec.Geometry.NodeSet{1}.ATTR.name; % febio_spec.Boundary.fix{1}.ATTR.bc='y'; % febio_spec.Boundary.fix{1}.ATTR.node_set=febio_spec.Geometry.NodeSet{1}.ATTR.name; % febio_spec.Boundary.fix{2}.ATTR.bc='z'; % febio_spec.Boundary.fix{2}.ATTR.node_set=febio_spec.Geometry.NodeSet{1}.ATTR.name; febio_spec.Loads.body_load{1}.ATTR.type='const'; febio_spec.Loads.body_load{1}.x.VAL=bodyLoadMagnitude; febio_spec.Loads.body_load{1}.x.ATTR.lc=1; % -> Prescribed boundary conditions on the rigid body febio_spec.Boundary.rigid_body{1}.ATTR.mat=2; febio_spec.Boundary.rigid_body{1}.fixed{1}.ATTR.bc='x'; febio_spec.Boundary.rigid_body{1}.fixed{2}.ATTR.bc='y'; febio_spec.Boundary.rigid_body{1}.fixed{3}.ATTR.bc='z'; febio_spec.Boundary.rigid_body{1}.fixed{4}.ATTR.bc='Rx'; febio_spec.Boundary.rigid_body{1}.fixed{5}.ATTR.bc='Ry'; febio_spec.Boundary.rigid_body{1}.fixed{6}.ATTR.bc='Rz'; %Contact section febio_spec.Contact.contact{1}.ATTR.surface_pair=febio_spec.Geometry.SurfacePair{1}.ATTR.name; febio_spec.Contact.contact{1}.ATTR.type='sliding-elastic'; febio_spec.Contact.contact{1}.two_pass=1; febio_spec.Contact.contact{1}.laugon=laugon; febio_spec.Contact.contact{1}.tolerance=0.2; febio_spec.Contact.contact{1}.gaptol=0; febio_spec.Contact.contact{1}.minaug=minaug; febio_spec.Contact.contact{1}.maxaug=maxaug; febio_spec.Contact.contact{1}.search_tol=0.01; febio_spec.Contact.contact{1}.search_radius=0.1; febio_spec.Contact.contact{1}.symmetric_stiffness=0; febio_spec.Contact.contact{1}.auto_penalty=1; febio_spec.Contact.contact{1}.penalty=contactPenalty; febio_spec.Contact.contact{1}.fric_coeff=fric_coeff; %LoadData febio_spec.LoadData.loadcurve{1}.ATTR.id=1; febio_spec.LoadData.loadcurve{1}.ATTR.type='linear'; febio_spec.LoadData.loadcurve{1}.point.VAL=[0 0; timeTotal1 1; timeTotal 1]; %Output section % -> log file febio_spec.Output.logfile.ATTR.file=febioLogFileName; febio_spec.Output.logfile.node_data{1}.ATTR.file=febioLogFileName_disp; febio_spec.Output.logfile.node_data{1}.ATTR.data='ux;uy;uz'; febio_spec.Output.logfile.node_data{1}.ATTR.delim=','; febio_spec.Output.logfile.node_data{1}.VAL=1:size(V,1); febio_spec.Output.logfile.element_data{1}.ATTR.file=febioLogFileName_strainEnergy; febio_spec.Output.logfile.element_data{1}.ATTR.data='sed'; febio_spec.Output.logfile.element_data{1}.ATTR.delim=','; febio_spec.Output.logfile.element_data{1}.VAL=1:size(E_blob,1);
Quick viewing of the FEBio input file structure
The febView function can be used to view the xml structure in a MATLAB figure window.
febView(febio_spec); %Viewing the febio file
Exporting the FEBio input file
Exporting the febio_spec structure to an FEBio input file is done using the febioStruct2xml function.
febioStruct2xml(febio_spec,febioFebFileName); %Exporting to file and domNode
Running the FEBio analysis
To run the analysis defined by the created FEBio input file the runMonitorFEBio function is used. The input for this function is a structure defining job settings e.g. the FEBio input file name. The optional output runFlag informs the user if the analysis was run succesfully.
febioAnalysis.run_filename=febioFebFileName; %The input file name febioAnalysis.run_logname=febioLogFileName; %The name for the log file febioAnalysis.disp_on=1; %Display information on the command window febioAnalysis.disp_log_on=1; %Display convergence information in the command window febioAnalysis.runMode='internal';%'internal'; febioAnalysis.t_check=0.25; %Time for checking log file (dont set too small) febioAnalysis.maxtpi=1e99; %Max analysis time febioAnalysis.maxLogCheckTime=3; %Max log file checking time [runFlag]=runMonitorFEBio(febioAnalysis);%START FEBio NOW!!!!!!!!
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
--- STARTING FEBIO JOB --- 04-Jun-2019 13:00:53
===========================================================================
________ _________ _________ __ _________
| |\ | |\ | |\ | |\ / \\
| ____|| | ____|| | __ || |__|| | ___ ||
| |\___\| | |\___\| | |\_| || \_\| | // \ ||
| || | || | || | || __ | || | ||
| ||__ | ||__ | ||_| || | |\ | || | ||
| |\ | |\ | \\ | || | || | ||
| ___|| | ___|| | ___ || | || | || | ||
| |\__\| | |\__\| | |\__| || | || | || | ||
| || | || | || | || | || | || | ||
| || | ||___ | ||__| || | || | \\__/ ||
| || | |\ | || | || | ||
|___|| |________|| |__________|| |__|| \_________//
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 . 9 . 0 ---
http://febio.org
FEBio is a registered trademark.
copyright (c) 2006-2019 - All rights reserved
===========================================================================
Success loading plugin NeoHookeanPi.so (version 0.0.0)
Reading file /mnt/data/MATLAB/GIBBON/data/temp/tempModel.feb ...SUCCESS!
*************************************************************************
* WARNING *
* *
* 230 isolated vertices removed. *
*************************************************************************
]0;(step 1/2: 0%) tempModel.feb - FEBio 2.9.0.0
===== beginning time step 1 : 0.01 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.01
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.599144e-10 1.697853e-08 0.000000e+00
energy 1.507741e-04 8.729688e-18 1.507741e-06
displacement 1.328960e+02 1.328960e+02 1.328960e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.01
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.599144e-10 8.788726e-11 0.000000e+00
energy 1.507741e-04 7.072029e-08 1.507741e-06
displacement 1.328960e+02 1.072782e-03 1.328306e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.01
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.599144e-10 5.701843e-12 0.000000e+00
energy 1.507741e-04 6.750887e-10 1.507741e-06
displacement 1.328960e+02 7.341109e-05 1.326744e-04
------- converged at time : 0.01
]0;(step 1/2: 1%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.028
===== beginning time step 2 : 0.038 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.038
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.902608e-09 9.778511e-06 0.000000e+00
energy 2.084601e-03 2.555929e-04 2.084601e-05
displacement 8.828945e+02 8.828945e+02 8.828945e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.038
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.902608e-09 1.114335e-06 0.000000e+00
energy 2.084601e-03 7.739712e-05 2.084601e-05
displacement 8.828945e+02 4.935043e-01 9.083835e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.038
stiffness updates = 0
right hand side evaluations = 5
stiffness matrix reformations = 3
step from line search = 0.563221
convergence norms : INITIAL CURRENT REQUIRED
residual 7.902608e-09 1.262631e-06 0.000000e+00
energy 2.084601e-03 4.965149e-04 2.084601e-05
displacement 8.828945e+02 2.073865e+02 1.953588e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.038
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 4
step from line search = 0.338180
convergence norms : INITIAL CURRENT REQUIRED
residual 7.902608e-09 3.972903e-06 0.000000e+00
energy 2.084601e-03 1.011565e-03 2.084601e-05
displacement 8.828945e+02 6.452235e+03 1.549308e-02
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.038
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 5
step from line search = 0.367730
convergence norms : INITIAL CURRENT REQUIRED
residual 7.902608e-09 2.458404e-06 0.000000e+00
energy 2.084601e-03 1.695133e-04 2.084601e-05
displacement 8.828945e+02 3.140316e+03 3.257357e-02
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.038
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 6
step from line search = 0.613408
convergence norms : INITIAL CURRENT REQUIRED
residual 7.902608e-09 4.203341e-07 0.000000e+00
energy 2.084601e-03 4.551389e-05 2.084601e-05
displacement 8.828945e+02 2.755608e+03 5.427285e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.038
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.902608e-09 2.689395e-09 0.000000e+00
energy 2.084601e-03 1.142567e-05 2.084601e-05
displacement 8.828945e+02 9.745975e+02 6.979038e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.038
stiffness updates = 0
right hand side evaluations = 15
stiffness matrix reformations = 8
step from line search = 0.500000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.902608e-09 6.711444e-03 0.000000e+00
energy 2.084601e-03 3.779249e-01 2.084601e-05
displacement 8.828945e+02 1.791939e+03 7.138767e-02
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
*************************************************************************
* ERROR *
* *
* Negative jacobian was detected. *
*************************************************************************
------- failed to converge at time : 0.038
Retrying time step. Retry attempt 1 of max 25
AUTO STEPPER: retry step, dt = 0.0269231
===== beginning time step 2 : 0.0369231 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.0369231
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.677452e-09 8.991032e-06 0.000000e+00
energy 2.015493e-03 2.213148e-04 2.015493e-05
displacement 8.495391e+02 8.495391e+02 8.495391e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.0369231
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.677452e-09 1.026379e-06 0.000000e+00
energy 2.015493e-03 6.991100e-05 2.015493e-05
displacement 8.495391e+02 4.922748e-01 8.741511e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 5
stiffness matrix reformations = 3
step from line search = 0.573360
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 1.201253e-06 0.000000e+00
energy 2.015493e-03 4.672237e-04 2.015493e-05
displacement 8.495391e+02 2.034385e+02 1.891690e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 4
step from line search = 0.356696
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 3.595495e-06 0.000000e+00
energy 2.015493e-03 9.899549e-04 2.015493e-05
displacement 8.495391e+02 6.384704e+03 1.521359e-02
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 5
step from line search = 0.382736
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 2.232301e-06 0.000000e+00
energy 2.015493e-03 1.141005e-04 2.015493e-05
displacement 8.495391e+02 2.876204e+03 3.131187e-02
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 6.246737e-08 0.000000e+00
energy 2.015493e-03 1.797287e-04 2.015493e-05
displacement 8.495391e+02 5.681936e+03 6.366285e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 1.455872e-10 0.000000e+00
energy 2.015493e-03 2.285142e-07 2.015493e-05
displacement 8.495391e+02 1.500138e+00 6.307343e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 1.231105e-09 0.000000e+00
energy 2.015493e-03 3.121425e-08 2.015493e-05
displacement 8.495391e+02 2.977469e+00 6.288499e-02
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 13
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 6.591845e-12 0.000000e+00
energy 2.015493e-03 1.935031e-08 2.015493e-05
displacement 8.495391e+02 2.024370e-01 6.295747e-02
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 15
stiffness matrix reformations = 10
step from line search = 0.065176
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 7.996353e-12 0.000000e+00
energy 2.015493e-03 9.661843e-09 2.015493e-05
displacement 8.495391e+02 7.412381e-02 6.294917e-02
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 18
stiffness matrix reformations = 11
step from line search = 0.500000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 5.732985e-06 0.000000e+00
energy 2.015493e-03 7.159359e-04 2.015493e-05
displacement 8.495391e+02 7.912762e+01 6.311145e-02
Reforming stiffness matrix: reformation #12
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
12
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 19
stiffness matrix reformations = 12
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 5.459216e-07 0.000000e+00
energy 2.015493e-03 3.435724e-05 2.015493e-05
displacement 8.495391e+02 4.595587e+02 5.304395e-02
Reforming stiffness matrix: reformation #13
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
13
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 20
stiffness matrix reformations = 13
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 1.092752e-08 0.000000e+00
energy 2.015493e-03 4.970191e-05 2.015493e-05
displacement 8.495391e+02 5.771681e+02 6.466135e-02
Reforming stiffness matrix: reformation #14
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
14
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 21
stiffness matrix reformations = 14
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 4.063209e-09 0.000000e+00
energy 2.015493e-03 1.498220e-06 2.015493e-05
displacement 8.495391e+02 8.664238e+00 6.332280e-02
Reforming stiffness matrix: reformation #15
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
15
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 22
stiffness matrix reformations = 15
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 3.755215e-11 0.000000e+00
energy 2.015493e-03 1.545779e-07 2.015493e-05
displacement 8.495391e+02 6.506702e-01 6.320109e-02
Reforming stiffness matrix: reformation #16
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
16
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 23
stiffness matrix reformations = 16
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 8.850617e-11 0.000000e+00
energy 2.015493e-03 6.370318e-08 2.015493e-05
displacement 8.495391e+02 9.927141e-01 6.305812e-02
Reforming stiffness matrix: reformation #17
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
17
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 24
stiffness matrix reformations = 17
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 1.640773e-11 0.000000e+00
energy 2.015493e-03 1.347894e-08 2.015493e-05
displacement 8.495391e+02 2.172922e-01 6.300701e-02
Reforming stiffness matrix: reformation #18
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
18
Nonlinear solution status: time= 0.0369231
stiffness updates = 0
right hand side evaluations = 25
stiffness matrix reformations = 18
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.677452e-09 6.357802e-13 0.000000e+00
energy 2.015493e-03 1.003658e-09 2.015493e-05
displacement 8.495391e+02 2.909546e-02 6.299103e-02
------- converged at time : 0.0369231
]0;(step 1/2: 4%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: decreasing time step, dt = 0.024586
===== beginning time step 3 : 0.0615091 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.0615091
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.989823e-08 9.207456e-07 0.000000e+00
energy 1.471545e-03 2.768945e-05 1.471545e-05
displacement 1.624454e+02 1.624454e+02 1.624454e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.0615091
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.989823e-08 6.556115e-07 0.000000e+00
energy 1.471545e-03 3.070859e-05 1.471545e-05
displacement 1.624454e+02 1.222823e-01 1.699668e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 3
step from line search = 0.094373
convergence norms : INITIAL CURRENT REQUIRED
residual 5.989823e-08 3.726886e-06 0.000000e+00
energy 1.471545e-03 9.456940e-04 1.471545e-05
displacement 1.624454e+02 2.724524e+02 7.920517e-04
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 4
step from line search = 0.076486
convergence norms : INITIAL CURRENT REQUIRED
residual 5.989823e-08 3.930388e-06 0.000000e+00
energy 1.471545e-03 1.342804e-03 1.471545e-05
displacement 1.624454e+02 2.165336e+02 1.826800e-03
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 5
step from line search = 0.298560
convergence norms : INITIAL CURRENT REQUIRED
residual 5.989823e-08 6.836010e-06 0.000000e+00
energy 1.471545e-03 5.077918e-03 1.471545e-05
displacement 1.624454e+02 8.229481e+03 1.749226e-02
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 13
stiffness matrix reformations = 6
step from line search = 0.445313
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 4.805400e-06 0.000000e+00
energy 5.077918e-03 1.374434e-03 5.077918e-05
displacement 1.624454e+02 9.984073e+03 5.360917e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 15
stiffness matrix reformations = 7
step from line search = 0.491043
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 1.597281e-06 0.000000e+00
energy 5.077918e-03 2.760012e-04 5.077918e-05
displacement 1.624454e+02 2.901666e+03 8.130017e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 17
stiffness matrix reformations = 8
step from line search = 0.500742
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 4.958983e-07 0.000000e+00
energy 5.077918e-03 4.108847e-05 5.077918e-05
displacement 1.624454e+02 6.768011e+02 9.666269e-02
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 19
stiffness matrix reformations = 9
step from line search = 0.509298
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 1.760353e-07 0.000000e+00
energy 5.077918e-03 3.581089e-06 5.077918e-05
displacement 1.624454e+02 1.634505e+02 1.045795e-01
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 21
stiffness matrix reformations = 10
step from line search = 0.480176
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 6.892976e-08 0.000000e+00
energy 5.077918e-03 4.887796e-07 5.077918e-05
displacement 1.624454e+02 3.954551e+01 1.084759e-01
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 22
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 2.356795e-08 0.000000e+00
energy 5.077918e-03 8.425802e-06 5.077918e-05
displacement 1.624454e+02 4.887303e+01 1.124695e-01
Reforming stiffness matrix: reformation #12
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
12
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 23
stiffness matrix reformations = 12
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 5.637217e-08 0.000000e+00
energy 5.077918e-03 1.392026e-06 5.077918e-05
displacement 1.624454e+02 1.567361e+01 1.128722e-01
Reforming stiffness matrix: reformation #13
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
13
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 25
stiffness matrix reformations = 13
step from line search = 0.433291
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 5.106814e-08 0.000000e+00
energy 5.077918e-03 5.183976e-06 5.077918e-05
displacement 1.624454e+02 5.657137e+00 1.126882e-01
Reforming stiffness matrix: reformation #14
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
14
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 27
stiffness matrix reformations = 14
step from line search = 0.429930
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 3.464488e-08 0.000000e+00
energy 5.077918e-03 3.191497e-06 5.077918e-05
displacement 1.624454e+02 4.768064e+00 1.129668e-01
Reforming stiffness matrix: reformation #15
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
15
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 29
stiffness matrix reformations = 15
step from line search = 0.606065
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 3.435006e-08 0.000000e+00
energy 5.077918e-03 1.633347e-06 5.077918e-05
displacement 1.624454e+02 4.938509e+00 1.128262e-01
Reforming stiffness matrix: reformation #16
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
16
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 30
stiffness matrix reformations = 16
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 1.109171e-08 0.000000e+00
energy 5.077918e-03 8.798075e-07 5.077918e-05
displacement 1.624454e+02 3.207285e+00 1.126334e-01
Reforming stiffness matrix: reformation #17
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
17
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 31
stiffness matrix reformations = 17
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 3.204210e-10 0.000000e+00
energy 5.077918e-03 1.689394e-08 5.077918e-05
displacement 1.624454e+02 5.602361e-01 1.126721e-01
Reforming stiffness matrix: reformation #18
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
18
Nonlinear solution status: time= 0.0615091
stiffness updates = 0
right hand side evaluations = 32
stiffness matrix reformations = 18
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.836010e-06 3.097521e-11 0.000000e+00
energy 5.077918e-03 4.937131e-09 5.077918e-05
displacement 1.624454e+02 3.644612e-02 1.125736e-01
------- converged at time : 0.0615091
]0;(step 1/2: 6%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: decreasing time step, dt = 0.0224526
===== beginning time step 4 : 0.0839616 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.0839616
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.341261e-07 2.086111e-07 0.000000e+00
energy 4.352606e-05 5.253149e-06 4.352606e-07
displacement 3.725582e+00 3.725582e+00 3.725582e-06
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.0839616
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.341261e-07 1.377541e-08 0.000000e+00
energy 4.352606e-05 1.400960e-06 4.352606e-07
displacement 3.725582e+00 3.400558e-03 3.936580e-06
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.0839616
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.341261e-07 2.151636e-07 0.000000e+00
energy 4.352606e-05 1.541856e-06 4.352606e-07
displacement 3.725582e+00 1.036593e+02 1.413730e-04
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.0839616
stiffness updates = 0
right hand side evaluations = 6
stiffness matrix reformations = 4
step from line search = 0.415950
convergence norms : INITIAL CURRENT REQUIRED
residual 1.341261e-07 2.365809e-07 0.000000e+00
energy 4.352606e-05 1.111924e-04 4.352606e-07
displacement 3.725582e+00 1.350828e+02 5.504596e-04
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.0839616
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 5
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.365809e-07 3.252715e-08 0.000000e+00
energy 1.111924e-04 3.529655e-05 1.111924e-06
displacement 3.725582e+00 7.334071e+02 2.473933e-03
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.0839616
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 6
step from line search = 0.530585
convergence norms : INITIAL CURRENT REQUIRED
residual 2.365809e-07 1.165371e-08 0.000000e+00
energy 1.111924e-04 1.428554e-06 1.111924e-06
displacement 3.725582e+00 3.712971e+00 2.524508e-03
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.0839616
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 7
step from line search = 0.349574
convergence norms : INITIAL CURRENT REQUIRED
residual 2.365809e-07 9.774803e-09 0.000000e+00
energy 1.111924e-04 2.916055e-07 1.111924e-06
displacement 3.725582e+00 1.263383e+00 2.521943e-03
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.0839616
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.365809e-07 1.006261e-08 0.000000e+00
energy 1.111924e-04 4.557483e-07 1.111924e-06
displacement 3.725582e+00 4.808376e+00 2.505868e-03
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.0839616
stiffness updates = 0
right hand side evaluations = 13
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.365809e-07 4.036172e-11 0.000000e+00
energy 1.111924e-04 4.999430e-08 1.111924e-06
displacement 3.725582e+00 1.686390e-01 2.485114e-03
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.0839616
stiffness updates = 0
right hand side evaluations = 14
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.365809e-07 2.311336e-12 0.000000e+00
energy 1.111924e-04 3.685602e-10 1.111924e-06
displacement 3.725582e+00 2.223933e-02 2.486099e-03
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 0.0839616
stiffness updates = 0
right hand side evaluations = 15
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.365809e-07 5.189645e-20 0.000000e+00
energy 1.111924e-04 1.519504e-14 1.111924e-06
displacement 3.725582e+00 8.686178e-06 2.486138e-03
------- converged at time : 0.0839616
]0;(step 1/2: 8%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.035461
===== beginning time step 5 : 0.119423 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.119423
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.534233e-07 1.402752e-07 0.000000e+00
energy 5.156159e-06 3.185756e-06 5.156159e-08
displacement 1.181919e-01 1.181919e-01 1.181919e-07
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.119423
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.534233e-07 1.159372e-10 0.000000e+00
energy 5.156159e-06 6.509155e-08 5.156159e-08
displacement 1.181919e-01 8.551572e-04 1.027280e-07
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.119423
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.534233e-07 8.622109e-13 0.000000e+00
energy 5.156159e-06 1.376144e-10 5.156159e-08
displacement 1.181919e-01 2.858702e-05 1.001295e-07
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.119423
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.534233e-07 8.946075e-14 0.000000e+00
energy 5.156159e-06 1.008993e-11 5.156159e-08
displacement 1.181919e-01 2.607417e-07 9.997706e-08
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.119423
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.534233e-07 2.619037e-20 0.000000e+00
energy 5.156159e-06 5.804985e-16 5.156159e-08
displacement 1.181919e-01 5.917563e-09 9.996248e-08
------- converged at time : 0.119423
]0;(step 1/2: 12%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0483688
===== beginning time step 6 : 0.167791 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.167791
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.185853e-07 9.764365e-07 0.000000e+00
energy 1.218371e-03 3.119300e-05 1.218371e-05
displacement 1.123151e+02 1.123151e+02 1.123151e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.167791
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.185853e-07 3.776131e-07 0.000000e+00
energy 1.218371e-03 1.941898e-05 1.218371e-05
displacement 1.123151e+02 4.499682e-02 1.160821e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.167791
stiffness updates = 0
right hand side evaluations = 5
stiffness matrix reformations = 3
step from line search = 0.381291
convergence norms : INITIAL CURRENT REQUIRED
residual 2.185853e-07 5.537895e-06 0.000000e+00
energy 1.218371e-03 2.623150e-03 1.218371e-05
displacement 1.123151e+02 1.360990e+03 2.196124e-03
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.167791
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 4
step from line search = 0.431359
convergence norms : INITIAL CURRENT REQUIRED
residual 5.537895e-06 1.173171e-05 0.000000e+00
energy 2.623150e-03 2.836693e-03 2.623150e-05
displacement 1.123151e+02 1.505721e+04 2.860008e-02
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.167791
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 5
step from line search = 0.481667
convergence norms : INITIAL CURRENT REQUIRED
residual 1.173171e-05 4.867478e-06 0.000000e+00
energy 2.836693e-03 4.606798e-04 2.836693e-05
displacement 1.123151e+02 5.346534e+03 5.860581e-02
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.167791
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.173171e-05 1.881659e-07 0.000000e+00
energy 2.836693e-03 4.000155e-04 2.836693e-05
displacement 1.123151e+02 4.474233e+03 9.540352e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.167791
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.173171e-05 1.915226e-08 0.000000e+00
energy 2.836693e-03 5.333865e-06 2.836693e-05
displacement 1.123151e+02 1.784030e+01 9.347229e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.167791
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.173171e-05 2.757780e-10 0.000000e+00
energy 2.836693e-03 2.415605e-07 2.836693e-05
displacement 1.123151e+02 7.353441e-01 9.322012e-02
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.167791
stiffness updates = 0
right hand side evaluations = 13
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.173171e-05 6.643294e-11 0.000000e+00
energy 2.836693e-03 1.979858e-08 2.836693e-05
displacement 1.123151e+02 2.416906e-01 9.331362e-02
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.167791
stiffness updates = 0
right hand side evaluations = 14
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.173171e-05 1.266279e-12 0.000000e+00
energy 2.836693e-03 2.350457e-10 2.836693e-05
displacement 1.123151e+02 3.506832e-03 9.330915e-02
------- converged at time : 0.167791
]0;(step 1/2: 17%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0586951
===== beginning time step 7 : 0.226487 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.226487
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.119683e-07 5.709267e-07 0.000000e+00
energy 2.430569e-05 1.305119e-05 2.430569e-07
displacement 9.751271e-01 9.751271e-01 9.751271e-07
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.226487
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.119683e-07 3.335684e-09 0.000000e+00
energy 2.430569e-05 8.979231e-07 2.430569e-07
displacement 9.751271e-01 4.115905e-03 1.093514e-06
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.226487
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.119683e-07 4.340265e-09 0.000000e+00
energy 2.430569e-05 5.802320e-07 2.430569e-07
displacement 9.751271e-01 2.445250e-02 1.406079e-06
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.226487
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.119683e-07 1.946210e-08 0.000000e+00
energy 2.430569e-05 3.633071e-06 2.430569e-07
displacement 9.751271e-01 2.936996e+00 8.091559e-06
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.226487
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.119683e-07 1.636525e-08 0.000000e+00
energy 2.430569e-05 5.451837e-06 2.430569e-07
displacement 9.751271e-01 4.239532e+02 5.472605e-04
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.226487
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.119683e-07 3.212268e-09 0.000000e+00
energy 2.430569e-05 4.967339e-07 2.430569e-07
displacement 9.751271e-01 3.843862e+00 6.040027e-04
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.226487
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.119683e-07 1.196352e-12 0.000000e+00
energy 2.430569e-05 1.331105e-09 2.430569e-07
displacement 9.751271e-01 1.012968e-02 6.045736e-04
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.226487
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 5.119683e-07 8.280814e-19 0.000000e+00
energy 2.430569e-05 5.220764e-13 2.430569e-07
displacement 9.751271e-01 4.311309e-04 6.053949e-04
------- converged at time : 0.226487
]0;(step 1/2: 23%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.066956
===== beginning time step 8 : 0.293443 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.293443
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.643006e-07 1.376086e-06 0.000000e+00
energy 7.558383e-04 3.670681e-05 7.558383e-06
displacement 6.194696e+01 6.194696e+01 6.194696e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.293443
stiffness updates = 0
right hand side evaluations = 3
stiffness matrix reformations = 2
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.643006e-07 2.109350e-07 0.000000e+00
energy 7.558383e-04 1.532985e-05 7.558383e-06
displacement 6.194696e+01 3.164310e-02 6.445829e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.293443
stiffness updates = 0
right hand side evaluations = 5
stiffness matrix reformations = 3
step from line search = 0.362067
convergence norms : INITIAL CURRENT REQUIRED
residual 6.643006e-07 3.726251e-06 0.000000e+00
energy 7.558383e-04 1.665479e-03 7.558383e-06
displacement 6.194696e+01 6.970947e+02 1.132606e-03
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.293443
stiffness updates = 0
right hand side evaluations = 6
stiffness matrix reformations = 4
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.726251e-06 2.116301e-05 0.000000e+00
energy 1.665479e-03 1.040577e-02 1.665479e-05
displacement 6.194696e+01 2.833532e+04 4.048689e-02
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.293443
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 5
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.116301e-05 3.106113e-07 0.000000e+00
energy 1.040577e-02 1.546900e-04 1.040577e-04
displacement 6.194696e+01 3.816238e+02 3.303666e-02
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.293443
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.116301e-05 1.794175e-07 0.000000e+00
energy 1.040577e-02 5.505994e-05 1.040577e-04
displacement 6.194696e+01 8.454761e+01 3.644933e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.293443
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.116301e-05 6.419511e-09 0.000000e+00
energy 1.040577e-02 1.107709e-06 1.040577e-04
displacement 6.194696e+01 6.631430e+00 3.599583e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.293443
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.116301e-05 3.375076e-11 0.000000e+00
energy 1.040577e-02 4.877580e-09 1.040577e-04
displacement 6.194696e+01 1.661360e-01 3.613932e-02
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.293443
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.116301e-05 1.687931e-11 0.000000e+00
energy 1.040577e-02 6.478863e-10 1.040577e-04
displacement 6.194696e+01 5.032454e-02 3.615545e-02
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.293443
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.116301e-05 1.858118e-17 0.000000e+00
energy 1.040577e-02 2.611419e-12 1.040577e-04
displacement 6.194696e+01 2.957379e-04 3.615276e-02
------- converged at time : 0.293443
]0;(step 1/2: 29%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0735648
===== beginning time step 9 : 0.367007 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.367007
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.145874e-06 1.674674e-06 0.000000e+00
energy 3.563771e-04 4.162467e-05 3.563771e-06
displacement 2.624348e+01 2.624348e+01 2.624348e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.367007
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.145874e-06 8.661589e-08 0.000000e+00
energy 3.563771e-04 9.706195e-06 3.563771e-06
displacement 2.624348e+01 2.018141e-02 2.757480e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.367007
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.145874e-06 2.476654e-06 0.000000e+00
energy 3.563771e-04 2.103218e-04 3.563771e-06
displacement 2.624348e+01 5.385010e+02 7.778781e-04
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.367007
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.145874e-06 1.415817e-06 0.000000e+00
energy 3.563771e-04 3.598081e-04 3.563771e-06
displacement 2.624348e+01 5.274548e+03 9.972296e-03
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.367007
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.415817e-06 6.661665e-08 0.000000e+00
energy 3.598081e-04 9.607740e-06 3.598081e-06
displacement 2.624348e+01 1.684997e+01 1.052963e-02
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.367007
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 6
step from line search = 0.568524
convergence norms : INITIAL CURRENT REQUIRED
residual 1.415817e-06 2.514448e-08 0.000000e+00
energy 3.598081e-04 1.240959e-07 3.598081e-06
displacement 2.624348e+01 2.971555e+00 1.070401e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.367007
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 7
step from line search = 0.498437
convergence norms : INITIAL CURRENT REQUIRED
residual 1.415817e-06 1.372853e-08 0.000000e+00
energy 3.598081e-04 1.212519e-07 3.598081e-06
displacement 2.624348e+01 1.837980e+00 1.079776e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.367007
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.415817e-06 8.912339e-09 0.000000e+00
energy 3.598081e-04 8.704714e-07 3.598081e-06
displacement 2.624348e+01 2.006743e+00 1.077073e-02
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.367007
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.415817e-06 2.761408e-10 0.000000e+00
energy 3.598081e-04 7.762756e-08 3.598081e-06
displacement 2.624348e+01 5.853571e-01 1.082123e-02
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.367007
stiffness updates = 0
right hand side evaluations = 13
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.415817e-06 2.773009e-15 0.000000e+00
energy 3.598081e-04 1.966836e-11 3.598081e-06
displacement 2.624348e+01 1.157842e-02 1.080007e-02
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 0.367007
stiffness updates = 0
right hand side evaluations = 14
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.415817e-06 6.047590e-21 0.000000e+00
energy 3.598081e-04 1.682881e-15 3.598081e-06
displacement 2.624348e+01 1.178741e-05 1.080022e-02
------- converged at time : 0.367007
]0;(step 1/2: 37%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0779993
===== beginning time step 10 : 0.445007 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.445007
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.591789e-06 2.143555e-06 0.000000e+00
energy 3.145862e-04 5.191920e-05 3.145862e-06
displacement 2.154271e+01 2.154271e+01 2.154271e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.445007
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.591789e-06 7.354687e-08 0.000000e+00
energy 3.145862e-04 9.647136e-06 3.145862e-06
displacement 2.154271e+01 1.991034e-02 2.273429e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.445007
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.591789e-06 6.307475e-07 0.000000e+00
energy 3.145862e-04 1.272766e-04 3.145862e-06
displacement 2.154271e+01 3.421505e+02 5.200553e-04
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.445007
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.591789e-06 1.934936e-06 0.000000e+00
energy 3.145862e-04 9.399018e-05 3.145862e-06
displacement 2.154271e+01 1.285245e+03 3.412149e-03
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.445007
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.591789e-06 3.187613e-07 0.000000e+00
energy 3.145862e-04 9.800431e-05 3.145862e-06
displacement 2.154271e+01 7.573171e+02 7.343170e-03
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.445007
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.591789e-06 2.886382e-07 0.000000e+00
energy 3.145862e-04 2.894476e-05 3.145862e-06
displacement 2.154271e+01 2.431069e+01 7.484143e-03
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.445007
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.591789e-06 6.731554e-09 0.000000e+00
energy 3.145862e-04 4.452685e-06 3.145862e-06
displacement 2.154271e+01 4.864671e+00 7.474114e-03
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.445007
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.591789e-06 2.170224e-08 0.000000e+00
energy 3.145862e-04 6.102948e-08 3.145862e-06
displacement 2.154271e+01 6.196969e+00 7.567733e-03
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.445007
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.591789e-06 1.017898e-10 0.000000e+00
energy 3.145862e-04 7.053544e-08 3.145862e-06
displacement 2.154271e+01 2.459049e-01 7.511071e-03
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.445007
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.591789e-06 6.513461e-12 0.000000e+00
energy 3.145862e-04 2.640812e-09 3.145862e-06
displacement 2.154271e+01 5.998353e-02 7.499801e-03
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 0.445007
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.591789e-06 5.575871e-20 0.000000e+00
energy 3.145862e-04 1.490810e-14 3.145862e-06
displacement 2.154271e+01 4.348692e-06 7.499778e-03
------- converged at time : 0.445007
]0;(step 1/2: 45%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0816899
===== beginning time step 11 : 0.526697 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.526697
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.188237e-06 3.413314e-06 0.000000e+00
energy 9.436156e-04 8.904978e-05 9.436156e-06
displacement 6.827881e+01 6.827881e+01 6.827881e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.526697
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.188237e-06 2.410742e-07 0.000000e+00
energy 9.436156e-04 2.311104e-05 9.436156e-06
displacement 6.827881e+01 5.314388e-02 7.177152e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.526697
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.188237e-06 3.711967e-05 0.000000e+00
energy 9.436156e-04 1.885081e-03 9.436156e-06
displacement 6.827881e+01 1.786925e+03 2.478447e-03
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 6
stiffness matrix reformations = 4
step from line search = 0.426018
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 7.655689e-06 0.000000e+00
energy 1.885081e-03 1.022932e-03 1.885081e-05
displacement 6.827881e+01 1.263687e+03 7.248857e-03
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 5
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 1.248690e-06 0.000000e+00
energy 1.885081e-03 2.460534e-04 1.885081e-05
displacement 6.827881e+01 2.521303e+03 1.806148e-02
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 6
step from line search = 0.300667
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 6.315842e-07 0.000000e+00
energy 1.885081e-03 2.467434e-05 1.885081e-05
displacement 6.827881e+01 9.094814e+00 1.811777e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 2.394031e-07 0.000000e+00
energy 1.885081e-03 3.377446e-05 1.885081e-05
displacement 6.827881e+01 2.117316e+01 1.842755e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 1.546043e-09 0.000000e+00
energy 1.885081e-03 1.186029e-06 1.885081e-05
displacement 6.827881e+01 3.536482e+00 1.798642e-02
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 3.598158e-08 0.000000e+00
energy 1.885081e-03 2.139977e-06 1.885081e-05
displacement 6.827881e+01 1.157766e+01 1.786008e-02
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 13
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 5.519163e-09 0.000000e+00
energy 1.885081e-03 1.365184e-08 1.885081e-05
displacement 6.827881e+01 2.117546e+00 1.796107e-02
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 14
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 4.105585e-09 0.000000e+00
energy 1.885081e-03 7.587870e-07 1.885081e-05
displacement 6.827881e+01 2.358320e+00 1.801978e-02
Reforming stiffness matrix: reformation #12
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
12
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 15
stiffness matrix reformations = 12
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 4.378627e-12 0.000000e+00
energy 1.885081e-03 5.546205e-09 1.885081e-05
displacement 6.827881e+01 6.887414e-02 1.799500e-02
Reforming stiffness matrix: reformation #13
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
13
Nonlinear solution status: time= 0.526697
stiffness updates = 0
right hand side evaluations = 16
stiffness matrix reformations = 13
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.711967e-05 3.739086e-13 0.000000e+00
energy 1.885081e-03 4.362901e-11 1.885081e-05
displacement 6.827881e+01 2.576250e-03 1.799736e-02
------- converged at time : 0.526697
]0;(step 1/2: 53%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.083048
===== beginning time step 12 : 0.609745 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.609745
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.028207e-06 3.358304e-06 0.000000e+00
energy 1.391572e-04 8.120023e-05 1.391572e-06
displacement 4.768184e+00 4.768184e+00 4.768184e-06
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.609745
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.028207e-06 1.782412e-08 0.000000e+00
energy 1.391572e-04 5.187162e-06 1.391572e-06
displacement 4.768184e+00 1.948040e-02 5.326982e-06
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.609745
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.028207e-06 6.474137e-08 0.000000e+00
energy 1.391572e-04 1.351569e-05 1.391572e-06
displacement 4.768184e+00 1.308620e+01 3.257252e-05
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.609745
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.028207e-06 7.298861e-08 0.000000e+00
energy 1.391572e-04 2.638647e-05 1.391572e-06
displacement 4.768184e+00 2.619022e+01 1.139604e-04
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.609745
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.028207e-06 1.262765e-07 0.000000e+00
energy 1.391572e-04 7.067578e-05 1.391572e-06
displacement 4.768184e+00 2.882061e+02 7.546114e-04
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.609745
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.028207e-06 1.155405e-08 0.000000e+00
energy 1.391572e-04 8.892530e-07 1.391572e-06
displacement 4.768184e+00 3.194336e+01 1.092962e-03
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.609745
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.028207e-06 8.011662e-12 0.000000e+00
energy 1.391572e-04 2.310367e-08 1.391572e-06
displacement 4.768184e+00 7.925783e-02 1.097242e-03
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.609745
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.028207e-06 8.397966e-19 0.000000e+00
energy 1.391572e-04 3.542232e-13 1.391572e-06
displacement 4.768184e+00 2.071014e-05 1.097311e-03
------- converged at time : 0.609745
]0;(step 1/2: 61%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0864384
===== beginning time step 13 : 0.696183 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.696183
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.663311e-06 5.712315e-06 0.000000e+00
energy 1.489492e-03 1.427597e-04 1.489492e-05
displacement 1.022174e+02 1.022174e+02 1.022174e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.696183
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.663311e-06 3.726031e-07 0.000000e+00
energy 1.489492e-03 3.726453e-05 1.489492e-05
displacement 1.022174e+02 8.570306e-02 1.073673e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.696183
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.663311e-06 3.127768e-05 0.000000e+00
energy 1.489492e-03 3.562412e-03 1.489492e-05
displacement 1.022174e+02 1.010683e+03 1.672068e-03
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 6
stiffness matrix reformations = 4
step from line search = 0.342694
convergence norms : INITIAL CURRENT REQUIRED
residual 3.127768e-05 1.632694e-05 0.000000e+00
energy 3.562412e-03 2.362357e-03 3.562412e-05
displacement 1.022174e+02 7.316594e+02 4.536425e-03
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 5
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 3.127768e-05 1.545248e-05 0.000000e+00
energy 3.562412e-03 4.547041e-03 3.562412e-05
displacement 1.022174e+02 5.468562e+03 1.985758e-02
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.545248e-05 4.298552e-07 0.000000e+00
energy 4.547041e-03 1.476582e-05 4.547041e-05
displacement 1.022174e+02 9.524429e+00 1.994592e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 7
step from line search = 0.433031
convergence norms : INITIAL CURRENT REQUIRED
residual 1.545248e-05 2.747546e-07 0.000000e+00
energy 4.547041e-03 1.382275e-05 4.547041e-05
displacement 1.022174e+02 1.456949e+01 2.001844e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 8
step from line search = 0.508950
convergence norms : INITIAL CURRENT REQUIRED
residual 1.545248e-05 1.665367e-07 0.000000e+00
energy 4.547041e-03 1.994136e-06 4.547041e-05
displacement 1.022174e+02 9.565892e+00 2.034345e-02
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 13
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.545248e-05 1.128445e-07 0.000000e+00
energy 4.547041e-03 8.220531e-06 4.547041e-05
displacement 1.022174e+02 1.101622e+01 2.043254e-02
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 14
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.545248e-05 8.094464e-08 0.000000e+00
energy 4.547041e-03 1.175372e-06 4.547041e-05
displacement 1.022174e+02 8.637970e+00 2.034524e-02
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 15
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.545248e-05 8.166680e-08 0.000000e+00
energy 4.547041e-03 1.253046e-06 4.547041e-05
displacement 1.022174e+02 4.571820e+00 2.031349e-02
Reforming stiffness matrix: reformation #12
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
12
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 16
stiffness matrix reformations = 12
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.545248e-05 2.264454e-10 0.000000e+00
energy 4.547041e-03 1.599920e-07 4.547041e-05
displacement 1.022174e+02 2.141343e-01 2.027450e-02
Reforming stiffness matrix: reformation #13
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
13
Nonlinear solution status: time= 0.696183
stiffness updates = 0
right hand side evaluations = 17
stiffness matrix reformations = 13
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.545248e-05 8.299446e-15 0.000000e+00
energy 4.547041e-03 2.692284e-10 4.547041e-05
displacement 1.022174e+02 7.026549e-03 2.029101e-02
------- converged at time : 0.696183
]0;(step 1/2: 70%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0874443
===== beginning time step 14 : 0.783627 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.783627
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 1
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 4.932346e-06 5.063397e-06 0.000000e+00
energy 1.351325e-04 1.285408e-04 1.351325e-06
displacement 6.396130e-01 6.396130e-01 6.396130e-07
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.783627
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 2
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 4.932346e-06 3.539473e-09 0.000000e+00
energy 1.351325e-04 2.130252e-06 1.351325e-06
displacement 6.396130e-01 2.544203e-02 8.929201e-07
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.783627
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 3
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 4.932346e-06 3.575766e-09 0.000000e+00
energy 1.351325e-04 6.147363e-07 1.351325e-06
displacement 6.396130e-01 2.551824e-02 1.181990e-06
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.783627
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 4
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 4.932346e-06 8.804475e-09 0.000000e+00
energy 1.351325e-04 1.407417e-06 1.351325e-06
displacement 6.396130e-01 1.855703e+00 5.722793e-06
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.783627
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 5
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 4.932346e-06 1.191218e-08 0.000000e+00
energy 1.351325e-04 2.745149e-06 1.351325e-06
displacement 6.396130e-01 8.129266e+00 2.709861e-05
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.783627
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 4.932346e-06 2.737997e-09 0.000000e+00
energy 1.351325e-04 1.641104e-06 1.351325e-06
displacement 6.396130e-01 1.130858e+01 7.270075e-05
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.783627
stiffness updates = 0
right hand side evaluations = 13
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 4.932346e-06 5.753560e-10 0.000000e+00
energy 1.351325e-04 2.056485e-07 1.351325e-06
displacement 6.396130e-01 1.207335e+00 8.837900e-05
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.783627
stiffness updates = 0
right hand side evaluations = 14
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 4.932346e-06 1.187531e-13 0.000000e+00
energy 1.351325e-04 2.540695e-09 1.351325e-06
displacement 6.396130e-01 7.491655e-03 8.967525e-05
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.783627
stiffness updates = 0
right hand side evaluations = 15
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 4.932346e-06 1.337260e-20 0.000000e+00
energy 1.351325e-04 1.783472e-14 1.351325e-06
displacement 6.396130e-01 1.562499e-06 8.969288e-05
------- converged at time : 0.783627
]0;(step 1/2: 78%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0899554
===== beginning time step 15 : 0.873583 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.873583
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.603651e-06 8.398399e-06 0.000000e+00
energy 1.872019e-03 2.297560e-04 1.872019e-05
displacement 1.213110e+02 1.213110e+02 1.213110e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.873583
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.603651e-06 4.718323e-07 0.000000e+00
energy 1.872019e-03 5.262994e-05 1.872019e-05
displacement 1.213110e+02 1.445335e-01 1.289513e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.873583
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.603651e-06 8.561941e-06 0.000000e+00
energy 1.872019e-03 8.667906e-05 1.872019e-05
displacement 1.213110e+02 1.010336e+03 1.752689e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
*************************************************************************
* ERROR *
* *
* Negative jacobian was detected. *
*************************************************************************
------- failed to converge at time : 0.873583
Retrying time step. Retry attempt 1 of max 25
AUTO STEPPER: retry step, dt = 0.0864956
===== beginning time step 15 : 0.870123 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.870123
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.591411e-06 8.323286e-06 0.000000e+00
energy 1.808860e-03 2.273961e-04 1.808860e-05
displacement 1.168722e+02 1.168722e+02 1.168722e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.870123
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.591411e-06 4.546053e-07 0.000000e+00
energy 1.808860e-03 5.130720e-05 1.808860e-05
displacement 1.168722e+02 1.402065e-01 1.242593e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.870123
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.591411e-06 7.806390e-06 0.000000e+00
energy 1.808860e-03 1.091960e-04 1.808860e-05
displacement 1.168722e+02 9.729608e+02 1.688201e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
*************************************************************************
* ERROR *
* *
* Negative jacobian was detected. *
*************************************************************************
------- failed to converge at time : 0.870123
Retrying time step. Retry attempt 2 of max 25
AUTO STEPPER: retry step, dt = 0.0830358
===== beginning time step 15 : 0.866663 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.866663
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.579285e-06 8.249049e-06 0.000000e+00
energy 1.747145e-03 2.250591e-04 1.747145e-05
displacement 1.125351e+02 1.125351e+02 1.125351e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.866663
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.579285e-06 4.377742e-07 0.000000e+00
energy 1.747145e-03 5.000317e-05 1.747145e-05
displacement 1.125351e+02 1.360079e-01 1.196755e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.866663
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.579285e-06 7.123132e-06 0.000000e+00
energy 1.747145e-03 1.291281e-04 1.747145e-05
displacement 1.125351e+02 9.363653e+02 1.625090e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
*************************************************************************
* ERROR *
* *
* Negative jacobian was detected. *
*************************************************************************
------- failed to converge at time : 0.866663
Retrying time step. Retry attempt 3 of max 25
AUTO STEPPER: retry step, dt = 0.079576
===== beginning time step 15 : 0.863203 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.863203
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.567280e-06 8.175722e-06 0.000000e+00
energy 1.686895e-03 2.227470e-04 1.686895e-05
displacement 1.083012e+02 1.083012e+02 1.083012e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.863203
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.567280e-06 4.213443e-07 0.000000e+00
energy 1.686895e-03 4.871842e-05 1.686895e-05
displacement 1.083012e+02 1.319365e-01 1.152011e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.863203
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.567280e-06 6.505223e-06 0.000000e+00
energy 1.686895e-03 1.467523e-04 1.686895e-05
displacement 1.083012e+02 9.005653e+02 1.563381e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
*************************************************************************
* ERROR *
* *
* Negative jacobian was detected. *
*************************************************************************
------- failed to converge at time : 0.863203
Retrying time step. Retry attempt 4 of max 25
AUTO STEPPER: retry step, dt = 0.0761161
===== beginning time step 15 : 0.859743 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.859743
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.555402e-06 8.103343e-06 0.000000e+00
energy 1.628133e-03 2.204614e-04 1.628133e-05
displacement 1.041718e+02 1.041718e+02 1.041718e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.859743
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.555402e-06 4.053221e-07 0.000000e+00
energy 1.628133e-03 4.745363e-05 1.628133e-05
displacement 1.041718e+02 1.279910e-01 1.108378e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.859743
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.555402e-06 5.946326e-06 0.000000e+00
energy 1.628133e-03 1.621651e-04 1.628133e-05
displacement 1.041718e+02 8.655731e+02 1.503096e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
*************************************************************************
* ERROR *
* *
* Negative jacobian was detected. *
*************************************************************************
------- failed to converge at time : 0.859743
Retrying time step. Retry attempt 5 of max 25
AUTO STEPPER: retry step, dt = 0.0726563
===== beginning time step 15 : 0.856284 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.856284
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.543661e-06 8.031956e-06 0.000000e+00
energy 1.570886e-03 2.182045e-04 1.570886e-05
displacement 1.001490e+02 1.001490e+02 1.001490e-04
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.856284
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.543661e-06 3.897150e-07 0.000000e+00
energy 1.570886e-03 4.620961e-05 1.570886e-05
displacement 1.001490e+02 1.241703e-01 1.065875e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.856284
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.543661e-06 5.440512e-06 0.000000e+00
energy 1.570886e-03 1.755191e-04 1.570886e-05
displacement 1.001490e+02 8.314023e+02 1.444258e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
*************************************************************************
* ERROR *
* *
* Negative jacobian was detected. *
*************************************************************************
------- failed to converge at time : 0.856284
Retrying time step. Retry attempt 6 of max 25
AUTO STEPPER: retry step, dt = 0.0691965
===== beginning time step 15 : 0.852824 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.852824
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.532069e-06 7.961617e-06 0.000000e+00
energy 1.515189e-03 2.159786e-04 1.515189e-05
displacement 9.623483e+01 9.623483e+01 9.623483e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.852824
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.532069e-06 3.745322e-07 0.000000e+00
energy 1.515189e-03 4.498733e-05 1.515189e-05
displacement 9.623483e+01 1.204737e-01 1.024522e-04
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.852824
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.532069e-06 4.982441e-06 0.000000e+00
energy 1.515189e-03 1.869272e-04 1.515189e-05
displacement 9.623483e+01 7.980680e+02 1.386894e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
*************************************************************************
* ERROR *
* *
* Negative jacobian was detected. *
*************************************************************************
------- failed to converge at time : 0.852824
Retrying time step. Retry attempt 7 of max 25
AUTO STEPPER: retry step, dt = 0.0657367
===== beginning time step 15 : 0.849364 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.849364
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.520639e-06 7.892395e-06 0.000000e+00
energy 1.461083e-03 2.137866e-04 1.461083e-05
displacement 9.243201e+01 9.243201e+01 9.243201e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.849364
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.520639e-06 3.597828e-07 0.000000e+00
energy 1.461083e-03 4.378985e-05 1.461083e-05
displacement 9.243201e+01 1.169008e-01 9.843489e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.849364
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.520639e-06 4.434325e-06 0.000000e+00
energy 1.461083e-03 2.136972e-04 1.461083e-05
displacement 9.243201e+01 7.616559e+02 1.325815e-03
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
*************************************************************************
* ERROR *
* *
* Negative jacobian was detected. *
*************************************************************************
------- failed to converge at time : 0.849364
Retrying time step. Retry attempt 8 of max 25
AUTO STEPPER: retry step, dt = 0.0622768
===== beginning time step 15 : 0.845904 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.845904
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.509387e-06 7.824379e-06 0.000000e+00
energy 1.408617e-03 2.116329e-04 1.408617e-05
displacement 8.874383e+01 8.874383e+01 8.874383e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.845904
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.509387e-06 3.450593e-07 0.000000e+00
energy 1.408617e-03 4.248371e-05 1.408617e-05
displacement 8.874383e+01 1.104548e-01 9.445763e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.845904
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.509387e-06 1.147470e-05 0.000000e+00
energy 1.408617e-03 1.434735e-03 1.408617e-05
displacement 8.874383e+01 1.037078e+03 1.653076e-03
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 6
stiffness matrix reformations = 4
step from line search = 0.506803
convergence norms : INITIAL CURRENT REQUIRED
residual 1.147470e-05 6.800182e-06 0.000000e+00
energy 1.434735e-03 1.697653e-03 1.434735e-05
displacement 8.874383e+01 6.351975e+02 4.247327e-03
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.845904
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.800182e-06 1.742451e-05 0.000000e+00
energy 1.697653e-03 3.920135e-03 1.697653e-05
displacement 8.874383e+01 3.418796e+03 1.516452e-02
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 6.627503e-07 0.000000e+00
energy 3.920135e-03 7.485310e-05 3.920135e-05
displacement 8.874383e+01 5.984919e+01 1.523558e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 1.314157e-06 0.000000e+00
energy 3.920135e-03 1.207135e-05 3.920135e-05
displacement 8.874383e+01 3.081882e+01 1.481612e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 1.178824e-07 0.000000e+00
energy 3.920135e-03 3.100553e-06 3.920135e-05
displacement 8.874383e+01 1.272310e+01 1.503768e-02
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 9
step from line search = 0.098609
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 1.454973e-07 0.000000e+00
energy 3.920135e-03 6.549147e-06 3.920135e-05
displacement 8.874383e+01 2.172132e+00 1.502842e-02
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 14
stiffness matrix reformations = 10
step from line search = 0.289170
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 1.686420e-07 0.000000e+00
energy 3.920135e-03 1.227378e-05 3.920135e-05
displacement 8.874383e+01 3.585025e+00 1.493372e-02
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 15
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 3.556702e-07 0.000000e+00
energy 3.920135e-03 1.476280e-05 3.920135e-05
displacement 8.874383e+01 1.933098e+01 1.476913e-02
Reforming stiffness matrix: reformation #12
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
12
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 16
stiffness matrix reformations = 12
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 3.445841e-09 0.000000e+00
energy 3.920135e-03 1.778558e-06 3.920135e-05
displacement 8.874383e+01 2.139559e+00 1.472380e-02
Reforming stiffness matrix: reformation #13
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
13
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 17
stiffness matrix reformations = 13
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 7.014014e-09 0.000000e+00
energy 3.920135e-03 1.014950e-09 3.920135e-05
displacement 8.874383e+01 3.668748e+00 1.476797e-02
Reforming stiffness matrix: reformation #14
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
14
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 18
stiffness matrix reformations = 14
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 5.023604e-12 0.000000e+00
energy 3.920135e-03 2.841482e-09 3.920135e-05
displacement 8.874383e+01 4.665676e-02 1.478834e-02
Reforming stiffness matrix: reformation #15
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
15
Nonlinear solution status: time= 0.845904
stiffness updates = 0
right hand side evaluations = 19
stiffness matrix reformations = 15
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 1.742451e-05 2.138347e-12 0.000000e+00
energy 3.920135e-03 3.875235e-10 3.920135e-05
displacement 8.874383e+01 6.490233e-03 1.478922e-02
------- converged at time : 0.845904
]0;(step 1/2: 85%) tempModel.feb - FEBio 2.9.0.0
===== beginning time step 16 : 0.908181 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.908181
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.826730e-06 6.380709e-06 0.000000e+00
energy 2.126528e-04 1.727730e-04 2.126528e-06
displacement 2.210435e+00 2.210435e+00 2.210435e-06
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.908181
stiffness updates = 0
right hand side evaluations = 3
stiffness matrix reformations = 2
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.826730e-06 2.255148e-10 0.000000e+00
energy 2.126528e-04 1.434318e-07 2.126528e-06
displacement 2.210435e+00 3.063889e-02 1.843560e-06
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.908181
stiffness updates = 0
right hand side evaluations = 4
stiffness matrix reformations = 3
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.826730e-06 5.486948e-12 0.000000e+00
energy 2.126528e-04 8.621547e-10 2.126528e-06
displacement 2.210435e+00 3.932351e-05 1.845689e-06
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.908181
stiffness updates = 0
right hand side evaluations = 5
stiffness matrix reformations = 4
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.826730e-06 3.467092e-13 0.000000e+00
energy 2.126528e-04 4.218761e-12 2.126528e-06
displacement 2.210435e+00 8.879292e-08 1.845458e-06
------- converged at time : 0.908181
]0;(step 1/2: 91%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0698215
===== beginning time step 17 : 0.978002 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 0.978002
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.947588e-06 9.748602e-06 0.000000e+00
energy 1.565355e-03 2.675485e-04 1.565355e-05
displacement 9.093095e+01 9.093095e+01 9.093095e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 0.978002
stiffness updates = 0
right hand side evaluations = 3
stiffness matrix reformations = 2
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.947588e-06 2.704939e-07 0.000000e+00
energy 1.565355e-03 4.184498e-05 1.565355e-05
displacement 9.093095e+01 1.191337e-01 9.674233e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 0.978002
stiffness updates = 0
right hand side evaluations = 4
stiffness matrix reformations = 3
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 6.947588e-06 2.418663e-05 0.000000e+00
energy 1.565355e-03 2.498716e-03 1.565355e-05
displacement 9.093095e+01 4.786462e+02 9.276220e-04
*************************************************************************
* WARNING *
* *
* Problem is diverging. Stiffness matrix will now be reformed *
*************************************************************************
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 0.978002
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.418663e-05 3.616418e-06 0.000000e+00
energy 2.498716e-03 6.414054e-04 2.498716e-05
displacement 9.093095e+01 2.007908e+03 5.498886e-03
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 0.978002
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.418663e-05 4.831099e-06 0.000000e+00
energy 2.498716e-03 7.955602e-04 2.498716e-05
displacement 9.093095e+01 7.534814e+02 1.016767e-02
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 0.978002
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.418663e-05 2.909085e-07 0.000000e+00
energy 2.498716e-03 3.933778e-05 2.498716e-05
displacement 9.093095e+01 2.298597e+01 1.037617e-02
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 0.978002
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.418663e-05 2.057425e-09 0.000000e+00
energy 2.498716e-03 3.047161e-06 2.498716e-05
displacement 9.093095e+01 9.345077e-01 1.054217e-02
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 0.978002
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 8
step from line search = 0.295577
convergence norms : INITIAL CURRENT REQUIRED
residual 2.418663e-05 8.747347e-09 0.000000e+00
energy 2.498716e-03 1.325411e-06 2.498716e-05
displacement 9.093095e+01 2.769779e+00 1.052484e-02
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 0.978002
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.418663e-05 2.277641e-08 0.000000e+00
energy 2.498716e-03 1.832155e-06 2.498716e-05
displacement 9.093095e+01 7.418165e+00 1.044699e-02
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 0.978002
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.418663e-05 1.653743e-10 0.000000e+00
energy 2.498716e-03 4.368804e-09 2.498716e-05
displacement 9.093095e+01 4.254184e-01 1.049164e-02
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 0.978002
stiffness updates = 0
right hand side evaluations = 13
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.418663e-05 9.885905e-12 0.000000e+00
energy 2.498716e-03 3.666891e-09 2.498716e-05
displacement 9.093095e+01 5.391407e-02 1.048185e-02
Reforming stiffness matrix: reformation #12
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
12
Nonlinear solution status: time= 0.978002
stiffness updates = 0
right hand side evaluations = 14
stiffness matrix reformations = 12
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 2.418663e-05 6.568804e-20 0.000000e+00
energy 2.498716e-03 5.807036e-15 2.498716e-05
displacement 9.093095e+01 3.556318e-06 1.048164e-02
------- converged at time : 0.978002
]0;(step 1/2: 98%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 0.0733836
MUST POINT CONTROLLER: adjusting time step. dt = 0.0219976
===== beginning time step 18 : 1 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
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 8.615132e-06 7.577065e-06 0.000000e+00
energy 3.835741e-04 2.189394e-04 3.835741e-06
displacement 9.263867e+00 9.263867e+00 9.263867e-06
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
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 8.615132e-06 1.083060e-10 0.000000e+00
energy 3.835741e-04 9.453388e-08 3.835741e-06
displacement 9.263867e+00 4.065310e-02 8.357645e-06
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
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 8.615132e-06 3.523152e-12 0.000000e+00
energy 3.835741e-04 2.140281e-10 3.835741e-06
displacement 9.263867e+00 6.711000e-06 8.356653e-06
------- converged at time : 1
]0;(step 1/2: 100%) tempModel.feb - FEBio 2.9.0.0
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 .................... : 18
Total number of equilibrium iterations ............ : 214
Average number of equilibrium iterations .......... : 11.8889
Total number of right hand evaluations ............ : 289
Total number of stiffness reformations ............ : 223
Time in linear solver: 0:01:00
===== beginning time step 1 : 11 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 11
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.828326e-06 9.056247e-06 0.000000e+00
energy 5.008226e-04 2.617542e-04 5.008226e-06
displacement 1.671743e+01 1.671743e+01 1.671743e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 11
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.828326e-06 1.161914e-08 0.000000e+00
energy 5.008226e-04 3.548855e-06 5.008226e-06
displacement 1.671743e+01 7.107126e-02 1.869787e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 11
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.828326e-06 7.084082e-09 0.000000e+00
energy 5.008226e-04 1.026821e-06 5.008226e-06
displacement 1.671743e+01 2.505310e-02 1.986896e-05
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 11
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.828326e-06 7.548564e-09 0.000000e+00
energy 5.008226e-04 2.005615e-06 5.008226e-06
displacement 1.671743e+01 6.530418e-01 2.641856e-05
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 11
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.828326e-06 1.447898e-08 0.000000e+00
energy 5.008226e-04 3.286943e-06 5.008226e-06
displacement 1.671743e+01 1.242636e+01 7.139781e-05
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 11
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.828326e-06 1.971738e-09 0.000000e+00
energy 5.008226e-04 3.378581e-06 5.008226e-06
displacement 1.671743e+01 1.800957e+01 1.574857e-04
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 11
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.828326e-06 3.363532e-11 0.000000e+00
energy 5.008226e-04 1.349024e-08 5.008226e-06
displacement 1.671743e+01 1.023219e-01 1.644730e-04
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 11
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.828326e-06 1.635589e-17 0.000000e+00
energy 5.008226e-04 4.070578e-12 5.008226e-06
displacement 1.671743e+01 1.759604e-04 1.647366e-04
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 11
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.828326e-06 9.114350e-28 0.000000e+00
energy 5.008226e-04 7.858586e-20 5.008226e-06
displacement 1.671743e+01 1.434448e-09 1.647373e-04
------- converged at time : 11
]0;(step 2/2: 10%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 18
===== beginning time step 2 : 29 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 29
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.772527e-06 7.497223e-06 0.000000e+00
energy 4.606602e-04 2.186738e-04 4.606602e-06
displacement 1.390709e+01 1.390709e+01 1.390709e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 29
stiffness updates = 0
right hand side evaluations = 3
stiffness matrix reformations = 2
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 8.772527e-06 1.432059e-10 0.000000e+00
energy 4.606602e-04 1.102643e-07 4.606602e-06
displacement 1.390709e+01 4.092538e-02 1.278231e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 29
stiffness updates = 0
right hand side evaluations = 4
stiffness matrix reformations = 3
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 8.772527e-06 5.686262e-12 0.000000e+00
energy 4.606602e-04 1.553576e-09 4.606602e-06
displacement 1.390709e+01 4.618043e-05 1.279583e-05
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 29
stiffness updates = 0
right hand side evaluations = 5
stiffness matrix reformations = 4
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 8.772527e-06 1.812005e-18 0.000000e+00
energy 4.606602e-04 9.352310e-15 4.606602e-06
displacement 1.390709e+01 2.527292e-07 1.279681e-05
------- converged at time : 29
]0;(step 2/2: 28%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 24.4
===== beginning time step 3 : 53.4 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 53.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 7.745603e-06 8.965873e-06 0.000000e+00
energy 4.985069e-04 2.607728e-04 4.985069e-06
displacement 1.656476e+01 1.656476e+01 1.656476e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 53.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 7.745603e-06 7.042087e-09 0.000000e+00
energy 4.985069e-04 1.551039e-06 4.985069e-06
displacement 1.656476e+01 6.677272e-02 1.846182e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 53.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 7.745603e-06 2.107854e-09 0.000000e+00
energy 4.985069e-04 2.998167e-07 4.985069e-06
displacement 1.656476e+01 5.558880e-03 1.897087e-05
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 53.4
stiffness updates = 0
right hand side evaluations = 5
stiffness matrix reformations = 4
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.745603e-06 1.084774e-09 0.000000e+00
energy 4.985069e-04 1.562763e-07 4.985069e-06
displacement 1.656476e+01 5.804983e-03 1.946080e-05
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 53.4
stiffness updates = 0
right hand side evaluations = 6
stiffness matrix reformations = 5
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.745603e-06 6.364827e-10 0.000000e+00
energy 4.985069e-04 9.274552e-08 4.985069e-06
displacement 1.656476e+01 1.367750e-02 2.022035e-05
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 53.4
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.745603e-06 5.843343e-10 0.000000e+00
energy 4.985069e-04 4.169320e-08 4.985069e-06
displacement 1.656476e+01 5.332569e-03 2.073175e-05
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 53.4
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.745603e-06 2.232117e-10 0.000000e+00
energy 4.985069e-04 2.098170e-08 4.985069e-06
displacement 1.656476e+01 1.303778e-03 2.099989e-05
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 53.4
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.745603e-06 7.509793e-11 0.000000e+00
energy 4.985069e-04 1.036671e-08 4.985069e-06
displacement 1.656476e+01 4.332982e-04 2.115547e-05
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 53.4
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.745603e-06 3.037881e-11 0.000000e+00
energy 4.985069e-04 5.368524e-09 4.985069e-06
displacement 1.656476e+01 2.755492e-04 2.127417e-05
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 53.4
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.745603e-06 6.590410e-12 0.000000e+00
energy 4.985069e-04 6.641732e-10 4.985069e-06
displacement 1.656476e+01 4.814626e-05 2.132470e-05
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 53.4
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.745603e-06 1.232832e-12 0.000000e+00
energy 4.985069e-04 1.169316e-10 4.985069e-06
displacement 1.656476e+01 1.661907e-05 2.134905e-05
------- converged at time : 53.4
]0;(step 2/2: 52%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 28.6944
===== beginning time step 4 : 82.0944 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 82.0944
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.745225e-06 7.471889e-06 0.000000e+00
energy 4.596253e-04 2.181509e-04 4.596253e-06
displacement 1.387030e+01 1.387030e+01 1.387030e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 82.0944
stiffness updates = 0
right hand side evaluations = 3
stiffness matrix reformations = 2
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 8.745225e-06 1.407607e-10 0.000000e+00
energy 4.596253e-04 1.097185e-07 4.596253e-06
displacement 1.387030e+01 4.064184e-02 1.275255e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 82.0944
stiffness updates = 0
right hand side evaluations = 4
stiffness matrix reformations = 3
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 8.745225e-06 4.416299e-12 0.000000e+00
energy 4.596253e-04 1.212198e-09 4.596253e-06
displacement 1.387030e+01 4.468559e-05 1.276584e-05
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 82.0944
stiffness updates = 0
right hand side evaluations = 5
stiffness matrix reformations = 4
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 8.745225e-06 5.926058e-18 0.000000e+00
energy 4.596253e-04 1.484225e-13 4.596253e-06
displacement 1.387030e+01 1.592037e-07 1.276652e-05
------- converged at time : 82.0944
]0;(step 2/2: 81%) tempModel.feb - FEBio 2.9.0.0
AUTO STEPPER: increasing time step, dt = 32.9555
MUST POINT CONTROLLER: adjusting time step. dt = 18.9056
===== beginning time step 5 : 101 =====
Reforming stiffness matrix: reformation #1
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
1
Nonlinear solution status: time= 101
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.720190e-06 8.939397e-06 0.000000e+00
energy 4.975961e-04 2.601753e-04 4.975961e-06
displacement 1.653878e+01 1.653878e+01 1.653878e-05
Reforming stiffness matrix: reformation #2
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
2
Nonlinear solution status: time= 101
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.720190e-06 6.976094e-09 0.000000e+00
energy 4.975961e-04 1.546485e-06 4.975961e-06
displacement 1.653878e+01 6.670764e-02 1.843355e-05
Reforming stiffness matrix: reformation #3
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
3
Nonlinear solution status: time= 101
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.720190e-06 2.102849e-09 0.000000e+00
energy 4.975961e-04 2.956376e-07 4.975961e-06
displacement 1.653878e+01 5.481069e-03 1.893831e-05
Reforming stiffness matrix: reformation #4
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
4
Nonlinear solution status: time= 101
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.720190e-06 1.110626e-09 0.000000e+00
energy 4.975961e-04 1.473033e-07 4.975961e-06
displacement 1.653878e+01 7.540801e-03 1.947603e-05
Reforming stiffness matrix: reformation #5
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
5
Nonlinear solution status: time= 101
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.720190e-06 5.470915e-10 0.000000e+00
energy 4.975961e-04 9.227332e-08 4.975961e-06
displacement 1.653878e+01 1.264883e-02 2.021351e-05
Reforming stiffness matrix: reformation #6
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
6
Nonlinear solution status: time= 101
stiffness updates = 0
right hand side evaluations = 7
stiffness matrix reformations = 6
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.720190e-06 4.501554e-10 0.000000e+00
energy 4.975961e-04 4.186744e-08 4.975961e-06
displacement 1.653878e+01 5.149069e-03 2.071805e-05
Reforming stiffness matrix: reformation #7
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
7
Nonlinear solution status: time= 101
stiffness updates = 0
right hand side evaluations = 8
stiffness matrix reformations = 7
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.720190e-06 2.111125e-10 0.000000e+00
energy 4.975961e-04 1.660778e-08 4.975961e-06
displacement 1.653878e+01 1.238686e-03 2.097888e-05
Reforming stiffness matrix: reformation #8
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
8
Nonlinear solution status: time= 101
stiffness updates = 0
right hand side evaluations = 9
stiffness matrix reformations = 8
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.720190e-06 5.381675e-11 0.000000e+00
energy 4.975961e-04 9.685583e-09 4.975961e-06
displacement 1.653878e+01 4.174180e-04 2.112931e-05
Reforming stiffness matrix: reformation #9
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
9
Nonlinear solution status: time= 101
stiffness updates = 0
right hand side evaluations = 10
stiffness matrix reformations = 9
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.720190e-06 2.344683e-11 0.000000e+00
energy 4.975961e-04 2.573663e-09 4.975961e-06
displacement 1.653878e+01 2.213744e-04 2.123082e-05
Reforming stiffness matrix: reformation #10
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
10
Nonlinear solution status: time= 101
stiffness updates = 0
right hand side evaluations = 11
stiffness matrix reformations = 10
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.720190e-06 4.301338e-12 0.000000e+00
energy 4.975961e-04 5.169433e-10 4.975961e-06
displacement 1.653878e+01 2.605434e-05 2.126561e-05
Reforming stiffness matrix: reformation #11
===== reforming stiffness matrix:
Nr of equations ........................... : 11784
Nr of nonzeroes in stiffness matrix ....... : 906480
11
Nonlinear solution status: time= 101
stiffness updates = 0
right hand side evaluations = 12
stiffness matrix reformations = 11
step from line search = 1.000000
convergence norms : INITIAL CURRENT REQUIRED
residual 7.720190e-06 1.939583e-12 0.000000e+00
energy 4.975961e-04 1.061025e-10 4.975961e-06
displacement 1.653878e+01 1.255913e-05 2.128429e-05
------- converged at time : 101
]0;(step 2/2: 100%) tempModel.feb - FEBio 2.9.0.0
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 .................... : 5
Total number of equilibrium iterations ............ : 39
Average number of equilibrium iterations .......... : 7.8
Total number of right hand evaluations ............ : 44
Total number of stiffness reformations ............ : 39
Time in linear solver: 0:01:09
Elapsed time : 0:02:30
N O R M A L T E R M I N A T I O N
Waiting for log file...
Proceeding to check log file...04-Jun-2019 13:03:24
------- converged at time : 0.01
------- converged at time : 0.0369231
------- converged at time : 0.0615091
------- converged at time : 0.0839616
------- converged at time : 0.119423
------- converged at time : 0.167791
------- converged at time : 0.226487
------- converged at time : 0.293443
------- converged at time : 0.367007
------- converged at time : 0.445007
------- converged at time : 0.526697
------- converged at time : 0.609745
------- converged at time : 0.696183
------- converged at time : 0.783627
------- converged at time : 0.845904
------- converged at time : 0.908181
------- converged at time : 0.978002
------- converged at time : 1
------- converged at time : 11
------- converged at time : 29
------- converged at time : 53.4
------- converged at time : 82.0944
------- converged at time : 101
--- Done --- 04-Jun-2019 13:03:25
Import FEBio results
if 1%runFlag==1 %i.e. a succesful run
% Importing nodal displacements from a log file [time_mat, N_disp_mat,~]=importFEBio_logfile(fullfile(savePath,febioLogFileName_disp)); %Nodal displacements time_mat=[0; time_mat(:)]; %Time N_disp_mat=N_disp_mat(:,2:end,:); sizImport=size(N_disp_mat); sizImport(3)=sizImport(3)+1; N_disp_mat_n=zeros(sizImport); N_disp_mat_n(:,:,2:end)=N_disp_mat; N_disp_mat=N_disp_mat_n; DN=N_disp_mat(:,:,end); DN_magnitude=sqrt(sum(DN.^2,2)); V_def=V+DN; V_DEF=N_disp_mat+repmat(V,[1 1 size(N_disp_mat,3)]); X_DEF=V_DEF(:,1,:); Y_DEF=V_DEF(:,2,:); Z_DEF=V_DEF(:,3,:); % [CF]=vertexToFaceMeasure(Fb_all,DN_magnitude);
Importing element strain energies from a log file
[~,E_energy,~]=importFEBio_logfile(fullfile(savePath,febioLogFileName_strainEnergy)); %Element stresses %Remove nodal index column E_energy=E_energy(:,2:end,:); %Add initial state i.e. zero displacement sizImport=size(E_energy); sizImport(3)=sizImport(3)+1; E_energy_mat_n=zeros(sizImport); E_energy_mat_n(:,:,2:end)=E_energy; E_energy=E_energy_mat_n;
indBlob=unique(Fb_blob(:));
t=linspace(0,2*pi,250)';
V_def_blob=V(indBlob,:)+N_disp_mat(indBlob,:,end);
[~,indMin]=min(V_def_blob(:,1));
[~,indMax]=max(V_def_blob(:,1));
xEnd=V_def_blob(indMin,1);
xStart=V_def_blob(indMax,1);
rEnd=sphereRadius+(xEnd.*((sphereRadius-tubeRadius(2))/tubeLength));
rStart=sphereRadius+(xStart.*((sphereRadius-tubeRadius(2))/tubeLength));
xMid=mean([xStart xEnd]);%sum([rStart rEnd].*[xStart xEnd])./sum([rStart rEnd]);
rMid=sphereRadius+(xMid.*((sphereRadius-tubeRadius(2))/tubeLength));
V_plot_xEnd=[xEnd*ones(size(t)) rEnd*cos(t) rEnd*sin(t)];
V_plot_xMid=[xMid*ones(size(t)) rMid*cos(t) rMid*sin(t)];
V_plot_xStart=[xStart*ones(size(t)) rStart*cos(t) rStart*sin(t)];
Plotting the simulated results using anim8 to visualize and animate deformations
% Create basic view and store graphics handle to initiate animation hf=cFigure; hold on; ht=gtitle(['Radial stretch: ',num2str(rMid/sphereRadius)]); hp1=gpatch(Fb_blob,V_def,DN_magnitude,'none',1); %Add graphics object to animate hp2=plotV(V_plot_xEnd ,'r-','LineWidth',3); hp3=plotV(V_plot_xMid ,'r-','LineWidth',3); hp4=plotV(V_plot_xStart,'r-','LineWidth',3); gpatch(F_tube,V_def,'kw','none',0.25); %Add graphics object to animate axisGeom(gca,fontSize); colormap(cMap); colorbar; caxis([0 max(DN_magnitude(:))]); caxis manual; axis([min(X_DEF(:)) max(X_DEF(:)) min(Y_DEF(:)) max(Y_DEF(:)) min(Z_DEF(:)) max(Z_DEF(:))]); camlight headlight; lighting gouraud; view(0,0); % view(-30,30); zoom(1.5); axis off; drawnow; LMid=1; % Set up animation features animStruct.Time=time_mat; %The time vector for qt=1:1:size(N_disp_mat,3) %Loop over time increments DN=N_disp_mat(:,:,qt); %Current displacement DN_magnitude=sqrt(sum(DN.^2,2)); V_def=V+DN; %Current nodal coordinates V_def_blob=V_def(indBlob,:); [~,indMin]=min(V_def_blob(:,1)); [~,indMax]=max(V_def_blob(:,1)); xEnd=V_def_blob(indMin,1); xStart=V_def_blob(indMax,1); rEnd=sphereRadius+(xEnd.*((sphereRadius-tubeRadius(2))/tubeLength)); rStart=sphereRadius+(xStart.*((sphereRadius-tubeRadius(2))/tubeLength)); xMid=mean([xStart xEnd]);%sum([rStart rEnd].*[xStart xEnd])./sum([rStart rEnd]); rMid=sphereRadius+(xMid.*((sphereRadius-tubeRadius(2))/tubeLength)); LMid=min(LMid,rMid/sphereRadius); V_plot_xEnd=[xEnd*ones(size(t)) rEnd*cos(t) rEnd*sin(t)]; V_plot_xMid=[xMid*ones(size(t)) rMid*cos(t) rMid*sin(t)]; V_plot_xStart=[xStart*ones(size(t)) rStart*cos(t) rStart*sin(t)]; %Set entries in animation structure animStruct.Handles{qt}=[hp1 hp1 hp2 hp2 hp2 hp3 hp3 hp3 hp4 hp4 hp4 ht]; %Handles of objects to animate animStruct.Props{qt}={'Vertices','CData','XData','YData','ZData','XData','YData','ZData','XData','YData','ZData','String'}; %Properties of objects to animate animStruct.Set{qt}={V_def,DN_magnitude,... V_plot_xEnd(:,1),V_plot_xEnd(:,2),V_plot_xEnd(:,3),... V_plot_xMid(:,1),V_plot_xMid(:,2),V_plot_xMid(:,3),... V_plot_xStart(:,1),V_plot_xStart(:,2),V_plot_xStart(:,3),... ['Radial stretch: ',num2str(rMid/sphereRadius)]}; %Property values for to set in order to animate end anim8(hf,animStruct); %Initiate animation feature drawnow;
end
GIBBON www.gibboncode.org
Kevin Mattheus Moerman, [email protected]
GIBBON footer text
License: https://github.com/gibbonCode/GIBBON/blob/master/LICENSE
GIBBON: The Geometry and Image-based Bioengineering add-On. A toolbox for image segmentation, image-based modeling, meshing, and finite element analysis.
Copyright (C) 2019 Kevin Mattheus Moerman
This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.