DEMO_febio_0072_pneunet_actuator_set_disp_01
Below is a demonstration for:
- Building geometry for a simple pneunet actuator
- Defining the boundary conditions
- Coding the febio structure
- Running the model
- Importing and visualizing the displacement and stress results
Contents
- Keywords
- Plot settings
- Control parameters
- Removed unused nodes and clean up index matrices
- Defining the boundary conditions
- Find top surface elements
- Split element regions (central plate and sides);
- Defining the FEBio input structure
- Quick viewing of the FEBio input file structure
- Exporting the FEBio input file
- Running the FEBio analysis
- Import FEBio results
Keywords
- febio_spec version 3.0
- febio, FEBio
- pressure loading
- hexahedral elements, hex8
- pneunet actuator
- soft robotic
- static, solid
- hyperelastic, Ogden
- displacement logfile
- stress logfile
clear; close all; clc;
Plot settings
fontSize=20; faceAlpha1=0.8; markerSize=40; markerSize2=20; lineWidth=3;
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_stress=[febioFebFileNamePart,'_stress_out.txt']; %Log file name for exporting stress febioLogFileName_force=[febioFebFileNamePart,'_force_out.txt']; %Log file name for exporting force %Load appliedPressure1=0.15; appliedPressure2=appliedPressure1/30; %Define applied displacement perturbation prescribedDisplacement_X=2; %Material parameter set c1=1; %Shear-modulus-like parameter m1=2; %Material parameter setting degree of non-linearity k_factor=100; %Bulk modulus factor k=c1*k_factor; %Bulk modulus c2=c1*2; %Shear-modulus-like parameter m2=2; %Material parameter setting degree of non-linearity k2=c2*k_factor; %Bulk modulus % FEA control settings numTimeSteps=25; %Number of time steps desired opt_iter=25; %Optimum number of iterations max_refs=opt_iter*2; %Max reforms max_ups=0; %Set to zero to use full-Newton iterations max_retries=5; %Maximum number of retires dtmin=(1/numTimeSteps)/100; %Minimum time step size dtmax=(1/numTimeSteps)*4; %Maximum time step size runMode='internal';%'internal';
actuatorLengthDesired=100; pointSpacing=2; periodSize=10; numPeriods=ceil(actuatorLengthDesired/periodSize); heightFactor=1; %To scale height wrt period size widthFactor=1; %To scale width wrt period size numElementsPeriod=ceil(periodSize./pointSpacing); if numElementsPeriod<4 numElementsPeriod=4; end numElementsLength=((numPeriods-1)*numElementsPeriod)+(numElementsPeriod-1); modelLength=numElementsLength.*pointSpacing; modelHeight_X=periodSize*heightFactor; modelWidth_Y=periodSize*widthFactor; numElementsHeight_X=ceil(modelHeight_X./pointSpacing); if numElementsHeight_X<6 numElementsHeight_X=6; end numElementsWidth_Y=ceil(modelWidth_Y./pointSpacing); if numElementsWidth_Y<1 numElementsWidth_Y=1; end boxDim=[modelHeight_X modelWidth_Y modelLength]; boxEl=[numElementsHeight_X numElementsWidth_Y numElementsLength]; [meshStruct]=hexMeshBox(boxDim,boxEl); E_bar=meshStruct.E; V_bar=meshStruct.V; F_bar=meshStruct.F; Fb_bar=meshStruct.Fb; Cb_bar=meshStruct.faceBoundaryMarker; VE_bar=patchCentre(E_bar,V_bar); CZ=VE_bar(:,3); CZ=CZ-min(CZ); CZ=CZ./max(CZ); CZ=round((CZ.*(numElementsLength-1)))+1; CW=VE_bar(:,1); CW=CW-min(CW); CW=CW./max(CW); CW=round((CW.*(numElementsHeight_X-1)))+1; CD=rem(CZ,numElementsPeriod); logicKeep1=~(CD==0 & CW>3); E1=E_bar(logicKeep1,:); F1=element2patch(E1); [indBoundary1]=tesBoundary(F1,V_bar); logicKeep2=any(ismember(E1,F1(indBoundary1,:)),2); F2=element2patch(E1(logicKeep2,:)); [indBoundary2]=tesBoundary(F2,V_bar); Fb=F2(indBoundary2,:); Cb=7*ones(size(Fb,1),1); for q=1:1:6 F_Cb1=Fb_bar(Cb_bar==q,:); logicNow=all(ismember(Fb,F_Cb1),2); Cb(logicNow)=q; end Cb(~any(ismember(Fb,F1(indBoundary1,:)),2))=0;
Removed unused nodes and clean up index matrices
[E,V,indFix2]=patchCleanUnused(E1(logicKeep2,:),V_bar); V(:,1)=V(:,1)-min(V(:,1)); Fb=indFix2(Fb); F=indFix2(F2); V2=V; V2(:,1)=-V2(:,1); V=[V;V2]; Fb=[Fb;fliplr(Fb)+size(V2,1)]; F=[F;fliplr(F)+size(V2,1)]; cShift=max(Cb(:))+1; Cb=[Cb;Cb+cShift]; indicesInnerSurface=unique(Fb(Cb==1,:)); logicElementsInner=any(ismember(E,indicesInnerSurface),2); E1=E; E2=E(:,[5 6 7 8 1 2 3 4]); E2=E2(~logicElementsInner,:); E=[E1;E2+size(V2,1)]; C=[ones(size(E1,1),1); 2*ones(size(E2,1),1);]; [F,V,~,ind2]=mergeVertices(F,V); Fb=ind2(Fb); E=ind2(E); indicesInnerSurface=ind2(indicesInnerSurface); actuatorLength=max(V(:,3))-min(V(:,3));
Defining the boundary conditions
The visualization of the model boundary shows colors for each side of the disc. These labels can be used to define boundary conditions.
%Define supported node sets bcSupportList=unique(Fb(ismember(Cb,[5 5+cShift]),:)); %Node set part of selected face bcPrescribeList=unique(Fb(ismember(Cb,[6 6+cShift]),:)); %Node set part of selected face %Get pressure faces F_pressure1=Fb(Cb==0,:); F_pressure2=Fb(Cb==cShift,:);
ind2=unique(E(C==2,:)); d=max(V(ind2,1)); V(ind2,1)=V(ind2,1)-d; V(:,1)=V(:,1)+d/2; % [E,V,indFix]=patchCleanUnused(E,V); % bcSupportList=indFix(bcSupportList); % bcPrescribeList=indFix(bcPrescribeList); % F_pressure1=indFix(F_pressure1); % F_pressure2=indFix(F_pressure2); [F,CF]=element2patch(E,C,'hex8'); [F,V,~,indFix]=mergeVertices(F,V); E=indFix(E); bcSupportList=indFix(bcSupportList); bcPrescribeList=indFix(bcPrescribeList); F_pressure1=indFix(F_pressure1); F_pressure2=indFix(F_pressure2); indb=tesBoundary(F,V); Fb=F(indb,:);
cFigure; hold on; gpatch(F,V,CF,'k',0.5); % plotV(V(indicesInnerSurface,:),'k.','MarkerSize',markerSize); axisGeom; colormap(turbo(250)); icolorbar; camlight headlight; gdrawnow;

Find top surface elements
logicTopSurface=all(ismember(Fb,bcPrescribeList),2); F_top=Fb(logicTopSurface,:); center_of_mass=mean(V(unique(F_top(:)),:),1);
Visualizing boundary conditions. Markers plotted on the semi-transparent model denote the nodes in the various boundary condition lists.
hf=cFigure; title('Boundary conditions','FontSize',fontSize); xlabel('X','FontSize',fontSize); ylabel('Y','FontSize',fontSize); zlabel('Z','FontSize',fontSize); hold on; gpatch(Fb,V,'w','none',0.5); hl(1)=plotV(V(bcSupportList,:),'k.','MarkerSize',markerSize); hl(2)=gpatch(F_pressure1,V,'r','k',1); hl(3)=gpatch(F_pressure2,V,'b','k',1); hl(4)=gpatch(F_top,V,'g','k',1); patchNormPlot(F_pressure1,V); legend(hl,{'BC full support','Pressure surface 1','Pressure surface 2','BC force'}); axisGeom(gca,fontSize); camlight headlight; gdrawnow;

Split element regions (central plate and sides);
XE=patchCentre(E,V(:,1)); logicElementsInner=abs(XE)<eps; E1=E(~logicElementsInner,:); %Other elements E2=E(logicElementsInner,:); %Inner element layer E=[E1;E2]; [F1]=element2patch(E1); [F2]=element2patch(E2);
cFigure; hold on; gpatch(F1,V,'bw','k',0.5); gpatch(F2,V,'rw','k',0.5); axisGeom; colormap(turbo(250)); icolorbar; camlight headlight; gdrawnow;

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='3.0'; %Module section febio_spec.Module.ATTR.type='solid'; %Create control structure for use by all steps stepStruct.Control.analysis='STATIC'; stepStruct.Control.time_steps=numTimeSteps; stepStruct.Control.step_size=1/numTimeSteps; stepStruct.Control.solver.max_refs=max_refs; stepStruct.Control.solver.max_ups=max_ups; stepStruct.Control.time_stepper.dtmin=dtmin; stepStruct.Control.time_stepper.dtmax=dtmax; stepStruct.Control.time_stepper.max_retries=max_retries; stepStruct.Control.time_stepper.opt_iter=opt_iter; %Add template based default settings to proposed control section [stepStruct.Control]=structComplete(stepStruct.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.step{1}.Control=stepStruct.Control; febio_spec.Step.step{1}.ATTR.id=1; febio_spec.Step.step{2}.Control=stepStruct.Control; febio_spec.Step.step{2}.ATTR.id=2; febio_spec.Step.step{3}.Control=stepStruct.Control; febio_spec.Step.step{3}.ATTR.id=3; febio_spec.Step.step{4}.Control=stepStruct.Control; febio_spec.Step.step{4}.ATTR.id=4; %Material section materialName1='Material1'; febio_spec.Material.material{1}.ATTR.name=materialName1; febio_spec.Material.material{1}.ATTR.type='Ogden'; 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}.k=k; materialName2='Material2'; febio_spec.Material.material{2}.ATTR.name=materialName2; febio_spec.Material.material{2}.ATTR.type='Ogden'; febio_spec.Material.material{2}.ATTR.id=2; febio_spec.Material.material{2}.c1=c2; febio_spec.Material.material{2}.m1=m2; febio_spec.Material.material{2}.c2=c2; febio_spec.Material.material{2}.m2=-m2; febio_spec.Material.material{2}.k=k2; materialName3='Material3'; febio_spec.Material.material{3}.ATTR.name=materialName3; febio_spec.Material.material{3}.ATTR.type='rigid body'; febio_spec.Material.material{3}.ATTR.id=3; febio_spec.Material.material{3}.density=1; febio_spec.Material.material{3}.center_of_mass=center_of_mass; %Mesh section % -> Nodes febio_spec.Mesh.Nodes{1}.ATTR.name='nodeSet_all'; %The node set name febio_spec.Mesh.Nodes{1}.node.ATTR.id=(1:size(V,1))'; %The node id's febio_spec.Mesh.Nodes{1}.node.VAL=V; %The nodel coordinates % -> Elements partName1='Part1'; febio_spec.Mesh.Elements{1}.ATTR.name=partName1; %Name of this part febio_spec.Mesh.Elements{1}.ATTR.type='hex8'; %Element type febio_spec.Mesh.Elements{1}.elem.ATTR.id=(1:1:size(E1,1))'; %Element id's febio_spec.Mesh.Elements{1}.elem.VAL=E1; %The element matrix partName2='Part2'; febio_spec.Mesh.Elements{2}.ATTR.name=partName2; %Name of this part febio_spec.Mesh.Elements{2}.ATTR.type='hex8'; %Element type febio_spec.Mesh.Elements{2}.elem.ATTR.id=size(E1,1)+(1:1:size(E2,1))'; %Element id's febio_spec.Mesh.Elements{2}.elem.VAL=E2; %The element matrix partName3='Part3'; febio_spec.Mesh.Elements{3}.ATTR.name=partName3; %Name of this part febio_spec.Mesh.Elements{3}.ATTR.type='quad4'; %Element type febio_spec.Mesh.Elements{3}.elem.ATTR.id=size(E1,1)+size(E2,1)+(1:1:size(F_top,1))'; %Element id's febio_spec.Mesh.Elements{3}.elem.VAL=F_top; %The element matrix % -> Surfaces surfaceName1='LoadedSurface1'; febio_spec.Mesh.Surface{1}.ATTR.name=surfaceName1; febio_spec.Mesh.Surface{1}.quad4.ATTR.id=(1:1:size(F_pressure1,1))'; febio_spec.Mesh.Surface{1}.quad4.VAL=F_pressure1; surfaceName2='LoadedSurface2'; febio_spec.Mesh.Surface{2}.ATTR.name=surfaceName2; febio_spec.Mesh.Surface{2}.quad4.ATTR.id=(1:1:size(F_pressure2,1))'; febio_spec.Mesh.Surface{2}.quad4.VAL=F_pressure2; % -> NodeSets nodeSetName1='bcSupportList'; febio_spec.Mesh.NodeSet{1}.ATTR.name=nodeSetName1; febio_spec.Mesh.NodeSet{1}.node.ATTR.id=bcSupportList(:); %MeshDomains section febio_spec.MeshDomains.SolidDomain{1}.ATTR.name=partName1; febio_spec.MeshDomains.SolidDomain{1}.ATTR.mat=materialName1; febio_spec.MeshDomains.SolidDomain{2}.ATTR.name=partName2; febio_spec.MeshDomains.SolidDomain{2}.ATTR.mat=materialName2; febio_spec.MeshDomains.ShellDomain.ATTR.name=partName3; febio_spec.MeshDomains.ShellDomain.ATTR.mat=materialName3; %Boundary condition section % -> Fix boundary conditions febio_spec.Boundary.bc{1}.ATTR.type='fix'; febio_spec.Boundary.bc{1}.ATTR.node_set=nodeSetName1; febio_spec.Boundary.bc{1}.dofs='x,y,z'; %Loads section % -> Surface load febio_spec.Loads.surface_load{1}.ATTR.type='pressure'; febio_spec.Loads.surface_load{1}.ATTR.surface=surfaceName1; febio_spec.Loads.surface_load{1}.pressure.ATTR.lc=1; febio_spec.Loads.surface_load{1}.pressure.VAL=appliedPressure1; febio_spec.Loads.surface_load{1}.symmetric_stiffness=1; febio_spec.Loads.surface_load{2}.ATTR.type='pressure'; febio_spec.Loads.surface_load{2}.ATTR.surface=surfaceName2; febio_spec.Loads.surface_load{2}.pressure.ATTR.lc=1; febio_spec.Loads.surface_load{2}.pressure.VAL=appliedPressure2; febio_spec.Loads.surface_load{2}.symmetric_stiffness=1; %Rigid section % ->Rigid body fix boundary conditions febio_spec.Rigid.rigid_constraint{1}.ATTR.name='RigidFix_1'; febio_spec.Rigid.rigid_constraint{1}.ATTR.type='fix'; febio_spec.Rigid.rigid_constraint{1}.rb=3; febio_spec.Rigid.rigid_constraint{1}.dofs='Ry'; % ->Rigid body prescribe boundary conditions febio_spec.Step.step{2}.Rigid.rigid_constraint{1}.ATTR.name='RigidPrescribe'; febio_spec.Step.step{2}.Rigid.rigid_constraint{1}.ATTR.type='prescribe'; febio_spec.Step.step{2}.Rigid.rigid_constraint{1}.rb=3; febio_spec.Step.step{2}.Rigid.rigid_constraint{1}.dof='Rx'; febio_spec.Step.step{2}.Rigid.rigid_constraint{1}.value.ATTR.lc=2; febio_spec.Step.step{2}.Rigid.rigid_constraint{1}.value.VAL=prescribedDisplacement_X; febio_spec.Step.step{2}.Rigid.rigid_constraint{1}.relative=1; febio_spec.Step.step{3}.Rigid.rigid_constraint{1}.ATTR.name='RigidPrescribe'; febio_spec.Step.step{3}.Rigid.rigid_constraint{1}.ATTR.type='prescribe'; febio_spec.Step.step{3}.Rigid.rigid_constraint{1}.rb=3; febio_spec.Step.step{3}.Rigid.rigid_constraint{1}.dof='Rx'; febio_spec.Step.step{3}.Rigid.rigid_constraint{1}.value.ATTR.lc=3; febio_spec.Step.step{3}.Rigid.rigid_constraint{1}.value.VAL=-prescribedDisplacement_X; febio_spec.Step.step{3}.Rigid.rigid_constraint{1}.relative=1; febio_spec.Step.step{4}.Rigid.rigid_constraint{1}.ATTR.name='RigidPrescribe'; febio_spec.Step.step{4}.Rigid.rigid_constraint{1}.ATTR.type='prescribe'; febio_spec.Step.step{4}.Rigid.rigid_constraint{1}.rb=3; febio_spec.Step.step{4}.Rigid.rigid_constraint{1}.dof='Rx'; febio_spec.Step.step{4}.Rigid.rigid_constraint{1}.value.ATTR.lc=4; febio_spec.Step.step{4}.Rigid.rigid_constraint{1}.value.VAL=-prescribedDisplacement_X; febio_spec.Step.step{4}.Rigid.rigid_constraint{1}.relative=1; %LoadData section % -> load_controller febio_spec.LoadData.load_controller{1}.ATTR.id=1; febio_spec.LoadData.load_controller{1}.ATTR.type='loadcurve'; febio_spec.LoadData.load_controller{1}.interpolate='LINEAR'; febio_spec.LoadData.load_controller{1}.points.point.VAL=[0 0; 1 1; 2 1; 3 1; 4 1]; febio_spec.LoadData.load_controller{2}.ATTR.id=2; febio_spec.LoadData.load_controller{2}.ATTR.type='loadcurve'; febio_spec.LoadData.load_controller{2}.interpolate='LINEAR'; febio_spec.LoadData.load_controller{2}.points.point.VAL=[0 0; 1 0; 2 1;]; febio_spec.LoadData.load_controller{3}.ATTR.id=3; febio_spec.LoadData.load_controller{3}.ATTR.type='loadcurve'; febio_spec.LoadData.load_controller{3}.interpolate='LINEAR'; febio_spec.LoadData.load_controller{3}.points.point.VAL=[0 0; 1 0; 2 0; 3 1]; febio_spec.LoadData.load_controller{4}.ATTR.id=4; febio_spec.LoadData.load_controller{4}.ATTR.type='loadcurve'; febio_spec.LoadData.load_controller{4}.interpolate='LINEAR'; febio_spec.LoadData.load_controller{4}.points.point.VAL=[0 0; 1 0; 2 0; 3 0; 4 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_stress; febio_spec.Output.logfile.element_data{1}.ATTR.data='s1'; febio_spec.Output.logfile.element_data{1}.ATTR.delim=','; febio_spec.Output.logfile.element_data{1}.VAL=1:size(E,1); febio_spec.Output.logfile.rigid_body_data{1}.ATTR.file=febioLogFileName_force; febio_spec.Output.logfile.rigid_body_data{1}.ATTR.data='Fx'; febio_spec.Output.logfile.rigid_body_data{1}.ATTR.delim=',';
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 % febView(febioFebFileName);
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.runMode=runMode; [runFlag]=runMonitorFEBio(febioAnalysis);%START FEBio NOW!!!!!!!!
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% --------> RUNNING/MONITORING FEBIO JOB <-------- 23-Apr-2021 19:43:27 FEBio path: /home/kevin/FEBioStudio/bin/febio3 # Attempt removal of existing log files 23-Apr-2021 19:43:27 * Removal succesful 23-Apr-2021 19:43:27 # Attempt removal of existing .xplt files 23-Apr-2021 19:43:27 * Removal succesful 23-Apr-2021 19:43:27 # Starting FEBio... 23-Apr-2021 19:43:27 Max. total analysis time is: Inf s =========================================================================== ________ _________ _______ __ _________ | |\ | |\ | \\ | |\ / \\ | ____|| | ____|| | __ || |__|| | ___ || | |\___\| | |\___\| | |\_| || \_\| | // \ || | ||__ | ||__ | ||_| || | |\ | || | || | |\ | |\ | \\ | || | || | || | ___|| | ___|| | ___ || | || | || | || | |\__\| | |\__\| | |\__| || | || | || | || | || | ||___ | ||__| || | || | \\__/ || | || | |\ | || | || | || |___|| |________|| |_________// |__|| \_________// 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 version 3.2.0 FEBio is a registered trademark. copyright (c) 2006-2020 - All rights reserved =========================================================================== Default linear solver: pardiso Success loading plugin libfebiochem_lnx64.so (version 0.0.0) Success loading plugin libfebioheat_lnx64.so (version 0.0.0) Reading file /mnt/data/MATLAB/GIBBON/data/temp/tempModel.feb ...SUCCESS! ************************************************************************* * WARNING * * * * 20 isolated vertices removed. * * * ************************************************************************* Setting parameter "pressure" to : 0 Setting parameter "pressure" to : 0 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 ]0;(step 1/4: 0%) tempModel.feb - FEBio 3.2.0 ===== beginning time step 1 : 0.04 ===== Setting parameter "pressure" to : 0.04 Setting parameter "pressure" to : 0.04 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 ===== reforming stiffness matrix: Nr of equations ........................... : 10013 Nr of nonzeroes in stiffness matrix ....... : 277827 1 Nonlinear solution status: time= 0.04 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 1 step from line search = 0.048997 convergence norms : INITIAL CURRENT REQUIRED residual 2.295302e-01 2.698756e+00 0.000000e+00 energy 2.652857e-01 7.313284e-03 2.652857e-03 displacement 8.006763e+04 1.922212e+02 1.922212e-04 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.04 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.295302e-01 1.562803e+00 0.000000e+00 energy 2.652857e-01 6.598725e-04 2.652857e-03 displacement 8.006763e+04 2.104440e+02 7.998178e-04 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.04 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.295302e-01 5.186576e-01 0.000000e+00 energy 2.652857e-01 1.254143e-03 2.652857e-03 displacement 8.006763e+04 6.681779e+01 1.323833e-03 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.04 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.295302e-01 4.180192e-01 0.000000e+00 energy 2.652857e-01 5.401625e-04 2.652857e-03 displacement 8.006763e+04 9.128212e+01 2.109357e-03 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.04 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.295302e-01 1.012538e-01 0.000000e+00 energy 2.652857e-01 4.753420e-04 2.652857e-03 displacement 8.006763e+04 3.335606e+01 2.671342e-03 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.04 stiffness updates = 0 right hand side evaluations = 10 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.295302e-01 9.929942e-02 0.000000e+00 energy 2.652857e-01 1.047963e-04 2.652857e-03 displacement 8.006763e+04 3.725865e+01 3.338417e-03 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.04 stiffness updates = 0 right hand side evaluations = 11 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.295302e-01 4.818248e-03 0.000000e+00 energy 2.652857e-01 6.232108e-05 2.652857e-03 displacement 8.006763e+04 7.683721e+00 3.665767e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.04 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.295302e-01 1.622374e-03 0.000000e+00 energy 2.652857e-01 5.187091e-06 2.652857e-03 displacement 8.006763e+04 4.557345e+00 3.928391e-03 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.04 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.295302e-01 1.198253e-06 0.000000e+00 energy 2.652857e-01 1.318987e-07 2.652857e-03 displacement 8.006763e+04 1.252880e-01 3.972812e-03 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.04 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.295302e-01 1.548937e-10 0.000000e+00 energy 2.652857e-01 1.037610e-10 2.652857e-03 displacement 8.006763e+04 1.280259e-03 3.977317e-03 convergence summary number of iterations : 10 number of reformations : 10 ------- converged at time : 0.04 Data Record #1 =========================================================================== Step = 1 Time = 0.04 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 1 Time = 0.04 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 1 Time = 0.04 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 4%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.064 ===== beginning time step 2 : 0.104 ===== Setting parameter "pressure" to : 0.104 Setting parameter "pressure" to : 0.104 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 1 step from line search = 0.500000 convergence norms : INITIAL CURRENT REQUIRED residual 5.907132e-01 6.231083e+04 0.000000e+00 energy 1.806716e-01 1.342208e+02 1.806716e-03 displacement 1.221495e+05 3.053737e+04 3.053737e-02 ************************************************************************* * WARNING * * * * Problem is diverging. Stiffness matrix will now be reformed * * * ************************************************************************* Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 1.198601e+01 0.000000e+00 energy 1.342208e+02 1.079822e+00 1.342208e+00 displacement 1.221495e+05 6.030032e+01 2.980058e-02 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 3 step from line search = 0.487059 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 7.021480e+00 0.000000e+00 energy 1.342208e+02 2.277605e-02 1.342208e+00 displacement 1.221495e+05 4.063719e+01 2.939586e-02 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 4 step from line search = 0.149482 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 8.849540e+00 0.000000e+00 energy 1.342208e+02 2.372289e-02 1.342208e+00 displacement 1.221495e+05 2.082227e+02 2.468011e-02 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 11 stiffness matrix reformations = 5 step from line search = 0.553787 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 2.108300e+01 0.000000e+00 energy 1.342208e+02 1.946482e-02 1.342208e+00 displacement 1.221495e+05 3.915100e+02 1.887282e-02 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 12 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 1.129822e+01 0.000000e+00 energy 1.342208e+02 2.317369e-02 1.342208e+00 displacement 1.221495e+05 3.002108e+02 1.442722e-02 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 13 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 2.696373e+01 0.000000e+00 energy 1.342208e+02 9.719095e-03 1.342208e+00 displacement 1.221495e+05 6.017286e+02 9.148635e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 14 stiffness matrix reformations = 8 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 2.766985e+00 0.000000e+00 energy 1.342208e+02 1.999120e-02 1.342208e+00 displacement 1.221495e+05 1.294972e+02 7.104326e-03 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 16 stiffness matrix reformations = 9 step from line search = 0.343674 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 6.269141e+00 0.000000e+00 energy 1.342208e+02 9.063896e-03 1.342208e+00 displacement 1.221495e+05 1.718725e+02 5.073204e-03 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 10 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 2.075519e+01 0.000000e+00 energy 1.342208e+02 7.609878e-03 1.342208e+00 displacement 1.221495e+05 4.633148e+02 2.479398e-03 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 18 stiffness matrix reformations = 11 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 1.341287e+00 0.000000e+00 energy 1.342208e+02 1.363437e-02 1.342208e+00 displacement 1.221495e+05 1.214931e+02 1.509604e-03 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 12 step from line search = 0.341714 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 3.894210e+00 0.000000e+00 energy 1.342208e+02 6.132659e-03 1.342208e+00 displacement 1.221495e+05 1.305400e+02 7.607773e-04 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 1.279628e+01 0.000000e+00 energy 1.342208e+02 4.699797e-03 1.342208e+00 displacement 1.221495e+05 3.826722e+02 8.433184e-05 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 22 stiffness matrix reformations = 14 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 7.932325e-01 0.000000e+00 energy 1.342208e+02 8.566616e-03 1.342208e+00 displacement 1.221495e+05 9.650132e+01 3.297273e-05 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 24 stiffness matrix reformations = 15 step from line search = 0.364374 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 2.408696e+00 0.000000e+00 energy 1.342208e+02 3.715915e-03 1.342208e+00 displacement 1.221495e+05 1.086059e+02 1.895922e-04 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 25 stiffness matrix reformations = 16 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 6.109649e+00 0.000000e+00 energy 1.342208e+02 7.755157e-04 1.342208e+00 displacement 1.221495e+05 2.721800e+02 8.813096e-04 Reforming stiffness matrix: reformation #17 17 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 26 stiffness matrix reformations = 17 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 5.434971e-01 0.000000e+00 energy 1.342208e+02 4.804868e-03 1.342208e+00 displacement 1.221495e+05 8.165008e+01 1.490553e-03 Reforming stiffness matrix: reformation #18 18 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 28 stiffness matrix reformations = 18 step from line search = 0.469227 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 1.544789e+00 0.000000e+00 energy 1.342208e+02 1.798150e-03 1.342208e+00 displacement 1.221495e+05 9.466198e+01 2.329182e-03 Reforming stiffness matrix: reformation #19 19 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 29 stiffness matrix reformations = 19 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 1.607571e+00 0.000000e+00 energy 1.342208e+02 1.658395e-03 1.342208e+00 displacement 1.221495e+05 1.418035e+02 3.613004e-03 Reforming stiffness matrix: reformation #20 20 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 30 stiffness matrix reformations = 20 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 5.487953e-01 0.000000e+00 energy 1.342208e+02 1.901135e-03 1.342208e+00 displacement 1.221495e+05 8.288979e+01 4.785577e-03 Reforming stiffness matrix: reformation #21 21 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 31 stiffness matrix reformations = 21 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 1.016976e+00 0.000000e+00 energy 1.342208e+02 2.162908e-04 1.342208e+00 displacement 1.221495e+05 1.129123e+02 6.363642e-03 Reforming stiffness matrix: reformation #22 22 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 32 stiffness matrix reformations = 22 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 8.615811e-02 0.000000e+00 energy 1.342208e+02 7.896480e-04 1.342208e+00 displacement 1.221495e+05 3.311078e+01 7.312212e-03 Reforming stiffness matrix: reformation #23 23 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 33 stiffness matrix reformations = 23 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 5.458314e-01 0.000000e+00 energy 1.342208e+02 5.022580e-04 1.342208e+00 displacement 1.221495e+05 8.255329e+01 8.944830e-03 Reforming stiffness matrix: reformation #24 24 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 34 stiffness matrix reformations = 24 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 2.441015e-03 0.000000e+00 energy 1.342208e+02 1.153288e-04 1.342208e+00 displacement 1.221495e+05 5.670228e+00 9.399866e-03 Reforming stiffness matrix: reformation #25 25 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 35 stiffness matrix reformations = 25 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 3.154793e-02 0.000000e+00 energy 1.342208e+02 4.639273e-05 1.342208e+00 displacement 1.221495e+05 1.985244e+01 1.028189e-02 Reforming stiffness matrix: reformation #26 26 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 36 stiffness matrix reformations = 26 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 5.097707e-07 0.000000e+00 energy 1.342208e+02 4.298392e-07 1.342208e+00 displacement 1.221495e+05 9.341430e-02 1.034383e-02 Reforming stiffness matrix: reformation #27 27 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 37 stiffness matrix reformations = 27 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 3.675021e-08 0.000000e+00 energy 1.342208e+02 5.968382e-09 1.342208e+00 displacement 1.221495e+05 2.083017e-02 1.037316e-02 Reforming stiffness matrix: reformation #28 28 Nonlinear solution status: time= 0.104 stiffness updates = 0 right hand side evaluations = 38 stiffness matrix reformations = 28 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.231083e+04 1.121553e-11 0.000000e+00 energy 1.342208e+02 1.998680e-11 1.342208e+00 displacement 1.221495e+05 3.636979e-05 1.037193e-02 convergence summary number of iterations : 28 number of reformations : 28 ------- converged at time : 0.104 Data Record #1 =========================================================================== Step = 2 Time = 0.104 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 2 Time = 0.104 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 2 Time = 0.104 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 10%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: decreasing time step, dt = 0.0604964 ===== beginning time step 3 : 0.164496 ===== Setting parameter "pressure" to : 0.164496 Setting parameter "pressure" to : 0.164496 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 1 step from line search = 0.500000 convergence norms : INITIAL CURRENT REQUIRED residual 5.324311e-01 8.123214e+04 0.000000e+00 energy 1.349515e-01 1.708780e+02 1.349515e-03 displacement 1.375729e+05 3.439323e+04 3.439323e-02 ************************************************************************* * WARNING * * * * Problem is diverging. Stiffness matrix will now be reformed * * * ************************************************************************* Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 2.779009e+01 0.000000e+00 energy 1.708780e+02 1.922981e+00 1.708780e+00 displacement 1.375729e+05 8.625508e+01 3.363714e-02 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 2.172856e+02 0.000000e+00 energy 1.708780e+02 5.084882e-01 1.708780e+00 displacement 1.375729e+05 1.141664e+03 2.402263e-02 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 2.682291e+00 0.000000e+00 energy 1.708780e+02 3.474478e-02 1.708780e+00 displacement 1.375729e+05 1.600447e+01 2.303427e-02 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 5 step from line search = 0.192880 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 8.007808e+00 0.000000e+00 energy 1.708780e+02 1.833996e-02 1.708780e+00 displacement 1.375729e+05 1.025685e+02 2.006727e-02 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 11 stiffness matrix reformations = 6 step from line search = 0.481606 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 2.075231e+01 0.000000e+00 energy 1.708780e+02 1.518406e-02 1.708780e+00 displacement 1.375729e+05 1.263640e+02 1.701953e-02 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 12 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 2.487740e+01 0.000000e+00 energy 1.708780e+02 1.394520e-02 1.708780e+00 displacement 1.375729e+05 2.470291e+02 1.317419e-02 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 13 stiffness matrix reformations = 8 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.016811e+01 0.000000e+00 energy 1.708780e+02 2.171726e-02 1.708780e+00 displacement 1.375729e+05 2.265163e+02 9.949716e-03 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 14 stiffness matrix reformations = 9 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 2.892511e+01 0.000000e+00 energy 1.708780e+02 6.093128e-03 1.708780e+00 displacement 1.375729e+05 3.963995e+02 6.379122e-03 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 15 stiffness matrix reformations = 10 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 2.365897e+00 0.000000e+00 energy 1.708780e+02 1.802437e-02 1.708780e+00 displacement 1.375729e+05 1.309009e+02 4.685931e-03 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 11 step from line search = 0.368589 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 6.186173e+00 0.000000e+00 energy 1.708780e+02 7.922529e-03 1.708780e+00 displacement 1.375729e+05 1.439374e+02 3.192473e-03 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 18 stiffness matrix reformations = 12 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.645447e+01 0.000000e+00 energy 1.708780e+02 2.571503e-03 1.708780e+00 displacement 1.375729e+05 3.704600e+02 1.396583e-03 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 19 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.520071e+00 0.000000e+00 energy 1.708780e+02 1.190643e-02 1.708780e+00 displacement 1.375729e+05 1.233330e+02 6.978142e-04 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 14 step from line search = 0.409727 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 4.053198e+00 0.000000e+00 energy 1.708780e+02 4.984183e-03 1.708780e+00 displacement 1.375729e+05 1.321886e+02 2.344825e-04 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 22 stiffness matrix reformations = 15 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 7.514962e+00 0.000000e+00 energy 1.708780e+02 1.548315e-03 1.708780e+00 displacement 1.375729e+05 2.818841e+02 3.236170e-05 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 23 stiffness matrix reformations = 16 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.259081e+00 0.000000e+00 energy 1.708780e+02 7.097031e-03 1.708780e+00 displacement 1.375729e+05 1.201827e+02 2.046691e-04 Reforming stiffness matrix: reformation #17 17 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 25 stiffness matrix reformations = 17 step from line search = 0.569771 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 2.978119e+00 0.000000e+00 energy 1.708780e+02 2.499095e-03 1.708780e+00 displacement 1.375729e+05 1.339638e+02 6.474586e-04 Reforming stiffness matrix: reformation #18 18 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 26 stiffness matrix reformations = 18 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.953299e+00 0.000000e+00 energy 1.708780e+02 3.571496e-03 1.708780e+00 displacement 1.375729e+05 1.527263e+02 1.415786e-03 Reforming stiffness matrix: reformation #19 19 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 27 stiffness matrix reformations = 19 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.788821e+00 0.000000e+00 energy 1.708780e+02 2.188917e-03 1.708780e+00 displacement 1.375729e+05 1.474767e+02 2.468222e-03 Reforming stiffness matrix: reformation #20 20 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 28 stiffness matrix reformations = 20 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 8.115546e-01 0.000000e+00 energy 1.708780e+02 2.162980e-03 1.708780e+00 displacement 1.375729e+05 1.002992e+02 3.557861e-03 Reforming stiffness matrix: reformation #21 21 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 29 stiffness matrix reformations = 21 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.043860e+00 0.000000e+00 energy 1.708780e+02 7.130614e-04 1.708780e+00 displacement 1.375729e+05 1.138978e+02 4.939642e-03 Reforming stiffness matrix: reformation #22 22 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 30 stiffness matrix reformations = 22 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.834291e-01 0.000000e+00 energy 1.708780e+02 1.032522e-03 1.708780e+00 displacement 1.375729e+05 4.832916e+01 5.962060e-03 Reforming stiffness matrix: reformation #23 23 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 31 stiffness matrix reformations = 23 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 5.546453e-01 0.000000e+00 energy 1.708780e+02 1.941228e-04 1.708780e+00 displacement 1.375729e+05 8.316376e+01 7.449708e-03 Reforming stiffness matrix: reformation #24 24 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 32 stiffness matrix reformations = 24 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.096337e-02 0.000000e+00 energy 1.708780e+02 2.324696e-04 1.708780e+00 displacement 1.375729e+05 1.207641e+01 8.060253e-03 Reforming stiffness matrix: reformation #25 25 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 33 stiffness matrix reformations = 25 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 8.578816e-02 0.000000e+00 energy 1.708780e+02 1.015600e-04 1.708780e+00 displacement 1.375729e+05 3.268593e+01 9.117257e-03 Reforming stiffness matrix: reformation #26 26 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 34 stiffness matrix reformations = 26 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 2.170131e-05 0.000000e+00 energy 1.708780e+02 4.517919e-06 1.708780e+00 displacement 1.375729e+05 5.906090e-01 9.264320e-03 Reforming stiffness matrix: reformation #27 27 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 35 stiffness matrix reformations = 27 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.049618e-05 0.000000e+00 energy 1.708780e+02 1.486675e-07 1.708780e+00 displacement 1.375729e+05 3.663803e-01 9.380985e-03 Reforming stiffness matrix: reformation #28 28 Nonlinear solution status: time= 0.164496 stiffness updates = 0 right hand side evaluations = 36 stiffness matrix reformations = 28 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.123214e+04 1.098117e-09 0.000000e+00 energy 1.708780e+02 9.499124e-10 1.708780e+00 displacement 1.375729e+05 1.131319e-03 9.374497e-03 convergence summary number of iterations : 28 number of reformations : 28 ------- converged at time : 0.164496 Data Record #1 =========================================================================== Step = 3 Time = 0.164496351 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 3 Time = 0.164496351 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 3 Time = 0.164496351 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 16%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: decreasing time step, dt = 0.0571857 ===== beginning time step 4 : 0.221682 ===== Setting parameter "pressure" to : 0.221682 Setting parameter "pressure" to : 0.221682 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 1 step from line search = 0.500000 convergence norms : INITIAL CURRENT REQUIRED residual 4.798471e-01 1.207051e+05 0.000000e+00 energy 8.689072e-02 2.458420e+02 8.689072e-04 displacement 1.648955e+05 4.122388e+04 4.122388e-02 ************************************************************************* * WARNING * * * * Problem is diverging. Stiffness matrix will now be reformed * * * ************************************************************************* Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 1.316244e+02 0.000000e+00 energy 2.458420e+02 5.959434e+00 2.458420e+00 displacement 1.648955e+05 1.692249e+02 4.038715e-02 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 5.851480e+01 0.000000e+00 energy 2.458420e+02 3.789729e-01 2.458420e+00 displacement 1.648955e+05 3.464884e+02 4.745303e-02 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 9.100997e+00 0.000000e+00 energy 2.458420e+02 3.262378e-02 2.458420e+00 displacement 1.648955e+05 4.440656e+02 3.875139e-02 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 2.027563e+01 0.000000e+00 energy 2.458420e+02 2.750678e-02 2.458420e+00 displacement 1.648955e+05 7.456757e+02 2.878530e-02 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 9.417092e+00 0.000000e+00 energy 2.458420e+02 2.061952e-02 2.458420e+00 displacement 1.648955e+05 2.599339e+02 2.358565e-02 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 11 stiffness matrix reformations = 7 step from line search = 0.544297 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 8.451319e+00 0.000000e+00 energy 2.458420e+02 1.202778e-02 2.458420e+00 displacement 1.648955e+05 3.825610e+02 1.798680e-02 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.221682 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.207051e+05 1.122303e+01 0.000000e+00 energy 2.458420e+02 1.368970e-02 2.458420e+00 displacement 1.648955e+05 4.072026e+02 1.298597e-02 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.221682 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.207051e+05 6.166175e+00 0.000000e+00 energy 2.458420e+02 1.587965e-02 2.458420e+00 displacement 1.648955e+05 4.228708e+02 8.735886e-03 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.221682 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.207051e+05 1.409427e+01 0.000000e+00 energy 2.458420e+02 3.822683e-03 2.458420e+00 displacement 1.648955e+05 3.870110e+02 5.448940e-03 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.221682 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.207051e+05 2.539081e+00 0.000000e+00 energy 2.458420e+02 1.449470e-02 2.458420e+00 displacement 1.648955e+05 2.783652e+02 3.276697e-03 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 12 step from line search = 0.584078 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 6.836284e+00 0.000000e+00 energy 2.458420e+02 5.635469e-03 2.458420e+00 displacement 1.648955e+05 1.850572e+02 1.907680e-03 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 18 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 4.400514e+00 0.000000e+00 energy 2.458420e+02 8.380694e-03 2.458420e+00 displacement 1.648955e+05 2.937314e+02 7.146415e-04 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 19 stiffness matrix reformations = 14 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 5.545142e+00 0.000000e+00 energy 2.458420e+02 5.338002e-03 2.458420e+00 displacement 1.648955e+05 2.534923e+02 1.285708e-04 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 15 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 2.376617e+00 0.000000e+00 energy 2.458420e+02 6.580792e-03 2.458420e+00 displacement 1.648955e+05 2.089557e+02 3.799285e-05 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 16 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 5.276621e+00 0.000000e+00 energy 2.458420e+02 8.832466e-04 2.458420e+00 displacement 1.648955e+05 2.404310e+02 4.014690e-04 Reforming stiffness matrix: reformation #17 17 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 22 stiffness matrix reformations = 17 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 6.817512e-01 0.000000e+00 energy 2.458420e+02 4.757734e-03 2.458420e+00 displacement 1.648955e+05 1.098658e+02 9.213475e-04 Reforming stiffness matrix: reformation #18 18 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 24 stiffness matrix reformations = 18 step from line search = 0.572704 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 1.836446e+00 0.000000e+00 energy 2.458420e+02 1.580680e-03 2.458420e+00 displacement 1.648955e+05 1.066091e+02 1.646185e-03 Reforming stiffness matrix: reformation #19 19 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 25 stiffness matrix reformations = 19 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 1.027560e+00 0.000000e+00 energy 2.458420e+02 2.256126e-03 2.458420e+00 displacement 1.648955e+05 1.198337e+02 2.648278e-03 Reforming stiffness matrix: reformation #20 20 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 26 stiffness matrix reformations = 20 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 9.767035e-01 0.000000e+00 energy 2.458420e+02 1.175933e-03 2.458420e+00 displacement 1.648955e+05 1.109853e+02 3.838291e-03 Reforming stiffness matrix: reformation #21 21 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 27 stiffness matrix reformations = 21 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 3.196014e-01 0.000000e+00 energy 2.458420e+02 1.124357e-03 2.458420e+00 displacement 1.648955e+05 6.554048e+01 4.903582e-03 Reforming stiffness matrix: reformation #22 22 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 28 stiffness matrix reformations = 22 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 4.253665e-01 0.000000e+00 energy 2.458420e+02 2.585781e-04 2.458420e+00 displacement 1.648955e+05 7.312437e+01 6.170808e-03 Reforming stiffness matrix: reformation #23 23 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 29 stiffness matrix reformations = 23 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 3.829065e-02 0.000000e+00 energy 2.458420e+02 3.306068e-04 2.458420e+00 displacement 1.648955e+05 2.267424e+01 6.939802e-03 Reforming stiffness matrix: reformation #24 24 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 30 stiffness matrix reformations = 24 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 7.827187e-02 0.000000e+00 energy 2.458420e+02 8.316508e-06 2.458420e+00 displacement 1.648955e+05 3.135514e+01 7.901987e-03 Reforming stiffness matrix: reformation #25 25 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 31 stiffness matrix reformations = 25 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 4.012335e-04 0.000000e+00 energy 2.458420e+02 1.702159e-05 2.458420e+00 displacement 1.648955e+05 2.424645e+00 8.180709e-03 Reforming stiffness matrix: reformation #26 26 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 32 stiffness matrix reformations = 26 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 1.147089e-04 0.000000e+00 energy 2.458420e+02 2.945895e-07 2.458420e+00 displacement 1.648955e+05 1.219455e+00 8.381300e-03 Reforming stiffness matrix: reformation #27 27 Nonlinear solution status: time= 0.221682 stiffness updates = 0 right hand side evaluations = 33 stiffness matrix reformations = 27 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.207051e+05 3.064950e-09 0.000000e+00 energy 2.458420e+02 2.026700e-09 2.458420e+00 displacement 1.648955e+05 1.152524e-03 8.375144e-03 convergence summary number of iterations : 27 number of reformations : 27 ------- converged at time : 0.221682 Data Record #1 =========================================================================== Step = 4 Time = 0.221682065 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 4 Time = 0.221682065 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 4 Time = 0.221682065 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 22%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: decreasing time step, dt = 0.0550421 ===== beginning time step 5 : 0.276724 ===== Setting parameter "pressure" to : 0.276724 Setting parameter "pressure" to : 0.276724 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 1 step from line search = 0.500000 convergence norms : INITIAL CURRENT REQUIRED residual 4.482838e-01 1.480659e+05 0.000000e+00 energy 6.025074e-02 2.951531e+02 6.025074e-04 displacement 1.802037e+05 4.505092e+04 4.505092e-02 ************************************************************************* * WARNING * * * * Problem is diverging. Stiffness matrix will now be reformed * * * ************************************************************************* Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 4.273453e+02 0.000000e+00 energy 2.951531e+02 1.435545e+01 2.951531e+00 displacement 1.802037e+05 2.730580e+02 4.417323e-02 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 3.068171e+01 0.000000e+00 energy 2.951531e+02 3.567316e-02 2.951531e+00 displacement 1.802037e+05 1.123636e+02 4.522453e-02 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 4 step from line search = 0.590358 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 6.031386e+00 0.000000e+00 energy 2.951531e+02 1.639917e-02 2.951531e+00 displacement 1.802037e+05 6.534369e+01 4.213400e-02 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 10 stiffness matrix reformations = 5 step from line search = 0.057916 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 6.063418e+00 0.000000e+00 energy 2.951531e+02 1.354853e-02 2.951531e+00 displacement 1.802037e+05 3.274524e+01 3.983350e-02 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 12 stiffness matrix reformations = 6 step from line search = 0.217182 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 9.644660e+00 0.000000e+00 energy 2.951531e+02 1.909116e-02 2.951531e+00 displacement 1.802037e+05 1.287982e+02 3.544336e-02 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 14 stiffness matrix reformations = 7 step from line search = 0.625396 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 1.406329e+01 0.000000e+00 energy 2.951531e+02 1.515095e-02 2.951531e+00 displacement 1.802037e+05 3.783791e+02 2.850880e-02 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 15 stiffness matrix reformations = 8 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 7.528904e+00 0.000000e+00 energy 2.951531e+02 2.541045e-02 2.951531e+00 displacement 1.802037e+05 3.771953e+02 2.233557e-02 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 16 stiffness matrix reformations = 9 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 1.893701e+01 0.000000e+00 energy 2.951531e+02 2.383072e-03 2.951531e+00 displacement 1.802037e+05 5.698594e+02 1.577924e-02 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 10 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 2.943919e+00 0.000000e+00 energy 2.951531e+02 2.280001e-02 2.951531e+00 displacement 1.802037e+05 2.331550e+02 1.218028e-02 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 19 stiffness matrix reformations = 11 step from line search = 0.478205 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 7.938929e+00 0.000000e+00 energy 2.951531e+02 9.334487e-03 2.951531e+00 displacement 1.802037e+05 2.289197e+02 9.074522e-03 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 12 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 1.015823e+01 0.000000e+00 energy 2.951531e+02 7.988144e-03 2.951531e+00 displacement 1.802037e+05 4.157502e+02 5.612459e-03 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 4.344029e+00 0.000000e+00 energy 2.951531e+02 1.416175e-02 2.951531e+00 displacement 1.802037e+05 2.455061e+02 3.515161e-03 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 22 stiffness matrix reformations = 14 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 1.530311e+01 0.000000e+00 energy 2.951531e+02 7.853974e-03 2.951531e+00 displacement 1.802037e+05 4.850801e+02 1.397997e-03 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 23 stiffness matrix reformations = 15 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 8.532596e-01 0.000000e+00 energy 2.951531e+02 1.057822e-02 2.951531e+00 displacement 1.802037e+05 1.053794e+02 7.412056e-04 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 25 stiffness matrix reformations = 16 step from line search = 0.325421 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 2.673857e+00 0.000000e+00 energy 2.951531e+02 4.751199e-03 2.951531e+00 displacement 1.802037e+05 1.175560e+02 2.768092e-04 Reforming stiffness matrix: reformation #17 17 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 26 stiffness matrix reformations = 17 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 9.715380e+00 0.000000e+00 energy 2.951531e+02 4.970124e-03 2.951531e+00 displacement 1.802037e+05 3.670379e+02 3.066364e-05 Reforming stiffness matrix: reformation #18 18 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 27 stiffness matrix reformations = 18 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 4.513690e-01 0.000000e+00 energy 2.951531e+02 6.065529e-03 2.951531e+00 displacement 1.802037e+05 7.703670e+01 1.623910e-04 Reforming stiffness matrix: reformation #19 19 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 29 stiffness matrix reformations = 19 step from line search = 0.344805 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 1.526873e+00 0.000000e+00 energy 2.951531e+02 2.623103e-03 2.951531e+00 displacement 1.802037e+05 8.972933e+01 4.774508e-04 Reforming stiffness matrix: reformation #20 20 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 30 stiffness matrix reformations = 20 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 4.236401e+00 0.000000e+00 energy 2.951531e+02 1.107602e-03 2.951531e+00 displacement 1.802037e+05 2.373932e+02 1.373174e-03 Reforming stiffness matrix: reformation #21 21 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 31 stiffness matrix reformations = 21 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 2.653330e-01 0.000000e+00 energy 2.951531e+02 2.954774e-03 2.951531e+00 displacement 1.802037e+05 5.899390e+01 1.996578e-03 Reforming stiffness matrix: reformation #22 22 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 33 stiffness matrix reformations = 22 step from line search = 0.456877 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 8.484132e-01 0.000000e+00 energy 2.951531e+02 1.045527e-03 2.951531e+00 displacement 1.802037e+05 7.165358e+01 2.820238e-03 Reforming stiffness matrix: reformation #23 23 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 34 stiffness matrix reformations = 23 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 7.921304e-01 0.000000e+00 energy 2.951531e+02 9.041335e-04 2.951531e+00 displacement 1.802037e+05 1.022261e+02 3.991786e-03 Reforming stiffness matrix: reformation #24 24 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 35 stiffness matrix reformations = 24 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 2.175476e-01 0.000000e+00 energy 2.951531e+02 8.781468e-04 2.951531e+00 displacement 1.802037e+05 5.295490e+01 4.961230e-03 Reforming stiffness matrix: reformation #25 25 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 36 stiffness matrix reformations = 25 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 3.034890e-01 0.000000e+00 energy 2.951531e+02 1.413173e-04 2.951531e+00 displacement 1.802037e+05 6.224774e+01 6.131861e-03 Reforming stiffness matrix: reformation #26 26 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 37 stiffness matrix reformations = 26 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 1.687172e-02 0.000000e+00 energy 2.951531e+02 1.933224e-04 2.951531e+00 displacement 1.802037e+05 1.511235e+01 6.754376e-03 Reforming stiffness matrix: reformation #27 27 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 38 stiffness matrix reformations = 27 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 2.841555e-02 0.000000e+00 energy 2.951531e+02 3.952218e-06 2.951531e+00 displacement 1.802037e+05 1.894620e+01 7.487190e-03 Reforming stiffness matrix: reformation #28 28 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 39 stiffness matrix reformations = 28 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 2.829913e-05 0.000000e+00 energy 2.951531e+02 2.734830e-06 2.951531e+00 displacement 1.802037e+05 6.926432e-01 7.631662e-03 Reforming stiffness matrix: reformation #29 29 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 40 stiffness matrix reformations = 29 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 9.512882e-07 0.000000e+00 energy 2.951531e+02 6.057919e-08 2.951531e+00 displacement 1.802037e+05 1.082629e-01 7.689160e-03 Reforming stiffness matrix: reformation #30 30 Nonlinear solution status: time= 0.276724 stiffness updates = 0 right hand side evaluations = 41 stiffness matrix reformations = 30 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.480659e+05 2.420406e-09 0.000000e+00 energy 2.951531e+02 1.521471e-09 2.951531e+00 displacement 1.802037e+05 9.571598e-04 7.683760e-03 convergence summary number of iterations : 30 number of reformations : 30 ------- converged at time : 0.276724 Data Record #1 =========================================================================== Step = 5 Time = 0.276724145 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 5 Time = 0.276724145 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 5 Time = 0.276724145 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 28%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: decreasing time step, dt = 0.0502812 ===== beginning time step 6 : 0.327005 ===== Setting parameter "pressure" to : 0.327005 Setting parameter "pressure" to : 0.327005 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.327005 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 1 step from line search = 0.500000 convergence norms : INITIAL CURRENT REQUIRED residual 3.772053e-01 7.667380e+05 0.000000e+00 energy 7.258176e-02 1.318535e+03 7.258176e-04 displacement 3.858004e+05 9.645009e+04 9.645009e-02 ************************************************************************* * WARNING * * * * Problem is diverging. Stiffness matrix will now be reformed * * * ************************************************************************* Reforming stiffness matrix: reformation #2 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected. * * * ************************************************************************* ------- failed to converge at time : 0.327005 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.041901 ===== beginning time step 6 : 0.318625 ===== Setting parameter "pressure" to : 0.318625 Setting parameter "pressure" to : 0.318625 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.318625 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 1 step from line search = 0.500000 convergence norms : INITIAL CURRENT REQUIRED residual 2.619485e-01 3.221634e+07 0.000000e+00 energy 8.204391e-01 3.418596e+04 8.204391e-03 displacement 2.186975e+06 5.467438e+05 5.467438e-01 ************************************************************************* * WARNING * * * * Problem is diverging. Stiffness matrix will now be reformed * * * ************************************************************************* Reforming stiffness matrix: reformation #2 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected. * * * ************************************************************************* ------- failed to converge at time : 0.318625 Retrying time step. Retry attempt 2 of max 5 AUTO STEPPER: retry step, dt = 0.0335208 ===== beginning time step 6 : 0.310245 ===== Setting parameter "pressure" to : 0.310245 Setting parameter "pressure" to : 0.310245 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.310245 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 1 step from line search = 0.055645 convergence norms : INITIAL CURRENT REQUIRED residual 1.676475e-01 2.640781e+00 0.000000e+00 energy 1.978035e-01 5.269290e-03 1.978035e-03 displacement 6.018977e+04 1.863724e+02 1.863724e-04 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.310245 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.676475e-01 7.237364e-01 0.000000e+00 energy 1.978035e-01 1.131259e-03 1.978035e-03 displacement 6.018977e+04 1.551686e+02 6.765041e-04 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.310245 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.676475e-01 7.166895e-01 0.000000e+00 energy 1.978035e-01 3.413976e-04 1.978035e-03 displacement 6.018977e+04 7.977213e+01 1.216009e-03 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.310245 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.676475e-01 7.119652e-02 0.000000e+00 energy 1.978035e-01 5.556901e-04 1.978035e-03 displacement 6.018977e+04 3.996341e+01 1.696341e-03 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.310245 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.676475e-01 2.210272e-01 0.000000e+00 energy 1.978035e-01 9.354546e-05 1.978035e-03 displacement 6.018977e+04 5.071425e+01 2.331829e-03 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.310245 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.676475e-01 1.440604e-03 0.000000e+00 energy 1.978035e-01 5.537107e-05 1.978035e-03 displacement 6.018977e+04 5.086029e+00 2.554417e-03 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.310245 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.676475e-01 3.058290e-03 0.000000e+00 energy 1.978035e-01 4.973656e-06 1.978035e-03 displacement 6.018977e+04 6.212007e+00 2.812077e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.310245 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.676475e-01 1.154366e-08 0.000000e+00 energy 1.978035e-01 3.085244e-09 1.978035e-03 displacement 6.018977e+04 3.936720e-04 2.810313e-03 convergence summary number of iterations : 8 number of reformations : 8 ------- converged at time : 0.310245 Data Record #1 =========================================================================== Step = 6 Time = 0.310244922 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 6 Time = 0.310244922 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 6 Time = 0.310244922 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 31%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.0588166 ===== beginning time step 7 : 0.369062 ===== Setting parameter "pressure" to : 0.369062 Setting parameter "pressure" to : 0.369062 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 1 step from line search = 0.040264 convergence norms : INITIAL CURRENT REQUIRED residual 5.188179e-01 2.155757e+00 0.000000e+00 energy 1.668360e-01 3.313725e-03 1.668360e-03 displacement 8.560362e+04 1.387811e+02 1.387811e-04 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.188179e-01 1.335534e+02 0.000000e+00 energy 1.668360e-01 2.826754e-01 1.668360e-03 displacement 8.560362e+04 2.069360e+03 1.145100e-03 ************************************************************************* * WARNING * * * * Problem is diverging. Stiffness matrix will now be reformed * * * ************************************************************************* Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 1.362291e-02 0.000000e+00 energy 2.826754e-01 1.372103e-03 2.826754e-03 displacement 8.560362e+04 1.591578e+00 1.207400e-03 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 4 step from line search = 0.098024 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 1.408491e-01 0.000000e+00 energy 2.826754e-01 2.093552e-03 2.826754e-03 displacement 8.560362e+04 3.491825e+01 1.647244e-03 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 11 stiffness matrix reformations = 5 step from line search = 0.270141 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 6.744706e-01 0.000000e+00 energy 2.826754e-01 1.707983e-03 2.826754e-03 displacement 8.560362e+04 5.772980e+01 2.314988e-03 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 12 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 3.874374e+00 0.000000e+00 energy 2.826754e-01 3.367448e-03 2.826754e-03 displacement 8.560362e+04 2.065651e+02 3.893122e-03 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 13 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 4.482806e-02 0.000000e+00 energy 2.826754e-01 1.261913e-03 2.826754e-03 displacement 8.560362e+04 2.415997e+01 4.527955e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 15 stiffness matrix reformations = 8 step from line search = 0.373127 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 2.660130e-01 0.000000e+00 energy 2.826754e-01 5.040562e-04 2.826754e-03 displacement 8.560362e+04 4.162450e+01 5.434571e-03 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 16 stiffness matrix reformations = 9 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 3.583795e-01 0.000000e+00 energy 2.826754e-01 1.617295e-04 2.826754e-03 displacement 8.560362e+04 6.756685e+01 6.710315e-03 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 10 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 2.148071e-02 0.000000e+00 energy 2.826754e-01 2.333142e-04 2.826754e-03 displacement 8.560362e+04 1.704060e+01 7.401907e-03 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 18 stiffness matrix reformations = 11 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 4.185014e-02 0.000000e+00 energy 2.826754e-01 1.513668e-05 2.826754e-03 displacement 8.560362e+04 2.292461e+01 8.246655e-03 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 19 stiffness matrix reformations = 12 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 3.553742e-05 0.000000e+00 energy 2.826754e-01 3.840113e-06 2.826754e-03 displacement 8.560362e+04 8.176248e-01 8.411414e-03 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 2.655284e-06 0.000000e+00 energy 2.826754e-01 1.362536e-07 2.826754e-03 displacement 8.560362e+04 1.807862e-01 8.489441e-03 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.369062 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 14 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.335534e+02 1.372971e-08 0.000000e+00 energy 2.826754e-01 6.110374e-09 2.826754e-03 displacement 8.560362e+04 2.858825e-03 8.479644e-03 convergence summary number of iterations : 14 number of reformations : 14 ------- converged at time : 0.369062 Data Record #1 =========================================================================== Step = 7 Time = 0.369061543 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 7 Time = 0.369061543 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 7 Time = 0.369061543 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 37%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.0790533 ===== beginning time step 8 : 0.448115 ===== Setting parameter "pressure" to : 0.448115 Setting parameter "pressure" to : 0.448115 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 1 step from line search = 0.064147 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 2.068917e+00 0.000000e+00 energy 4.972056e-01 2.770191e-02 4.972056e-03 displacement 2.612144e+04 1.074845e+02 1.074845e-04 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 2 step from line search = 0.303978 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 2.502275e+01 0.000000e+00 energy 4.972056e-01 9.275139e-02 4.972056e-03 displacement 2.612144e+04 7.232271e+02 2.773106e-04 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 1.067496e+00 0.000000e+00 energy 4.972056e-01 9.569836e-03 4.972056e-03 displacement 2.612144e+04 1.236450e+02 7.697236e-04 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 4 step from line search = 0.295210 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 2.349303e+00 0.000000e+00 energy 4.972056e-01 3.979274e-03 4.972056e-03 displacement 2.612144e+04 1.004776e+02 1.414236e-03 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 8.695335e+00 0.000000e+00 energy 4.972056e-01 5.954660e-03 4.972056e-03 displacement 2.612144e+04 3.819505e+02 3.252626e-03 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 10 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 2.338375e-01 0.000000e+00 energy 4.972056e-01 4.235171e-03 4.972056e-03 displacement 2.612144e+04 6.387879e+01 4.223944e-03 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 12 stiffness matrix reformations = 7 step from line search = 0.316280 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 9.363983e-01 0.000000e+00 energy 4.972056e-01 1.866026e-03 4.972056e-03 displacement 2.612144e+04 7.232732e+01 5.395369e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 13 stiffness matrix reformations = 8 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 3.072384e+00 0.000000e+00 energy 4.972056e-01 1.463507e-03 4.972056e-03 displacement 2.612144e+04 2.112085e+02 7.732641e-03 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 14 stiffness matrix reformations = 9 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 1.039371e-01 0.000000e+00 energy 4.972056e-01 1.669977e-03 4.972056e-03 displacement 2.612144e+04 3.876084e+01 8.862234e-03 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 16 stiffness matrix reformations = 10 step from line search = 0.447751 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 4.135317e-01 0.000000e+00 energy 4.972056e-01 5.437294e-04 4.972056e-03 displacement 2.612144e+04 5.163653e+01 1.026188e-02 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 11 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 3.008192e-01 0.000000e+00 energy 4.972056e-01 4.362093e-04 4.972056e-03 displacement 2.612144e+04 6.469005e+01 1.195114e-02 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 18 stiffness matrix reformations = 12 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 5.739417e-02 0.000000e+00 energy 4.972056e-01 2.887427e-04 4.972056e-03 displacement 2.612144e+04 2.770217e+01 1.312619e-02 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 19 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 3.406938e-02 0.000000e+00 energy 4.972056e-01 4.825395e-05 4.972056e-03 displacement 2.612144e+04 2.117689e+01 1.419888e-02 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 14 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 4.675196e-04 0.000000e+00 energy 4.972056e-01 1.029514e-05 4.972056e-03 displacement 2.612144e+04 2.700589e+00 1.459229e-02 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 15 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 9.418255e-06 0.000000e+00 energy 4.972056e-01 1.373742e-07 4.972056e-03 displacement 2.612144e+04 3.641127e-01 1.473812e-02 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 0.448115 stiffness updates = 0 right hand side evaluations = 22 stiffness matrix reformations = 16 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.459875e-01 2.955427e-08 0.000000e+00 energy 4.972056e-01 1.139610e-08 4.972056e-03 displacement 2.612144e+04 3.274416e-03 1.472436e-02 convergence summary number of iterations : 16 number of reformations : 16 ------- converged at time : 0.448115 Data Record #1 =========================================================================== Step = 8 Time = 0.44811484 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 8 Time = 0.44811484 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 8 Time = 0.44811484 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 45%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.0952426 ===== beginning time step 9 : 0.543357 ===== Setting parameter "pressure" to : 0.543357 Setting parameter "pressure" to : 0.543357 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 1 step from line search = 0.118104 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 6.231512e+00 0.000000e+00 energy 7.977507e-01 7.241910e-02 7.977507e-03 displacement 1.662545e+04 2.319005e+02 2.319005e-04 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 1.568538e+02 0.000000e+00 energy 7.977507e-01 3.433476e-01 7.977507e-03 displacement 1.662545e+04 2.650646e+03 1.352701e-03 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 6.977920e-02 0.000000e+00 energy 7.977507e-01 5.317204e-03 7.977507e-03 displacement 1.662545e+04 6.256307e+00 1.485168e-03 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 4 step from line search = 0.080334 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 3.937702e-01 0.000000e+00 energy 7.977507e-01 4.969161e-03 7.977507e-03 displacement 1.662545e+04 4.835218e+01 2.052938e-03 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 5 step from line search = 0.217162 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 1.609909e+00 0.000000e+00 energy 7.977507e-01 4.664153e-03 7.977507e-03 displacement 1.662545e+04 8.259880e+01 2.939685e-03 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 11 stiffness matrix reformations = 6 step from line search = 0.537998 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 3.734768e+00 0.000000e+00 energy 7.977507e-01 3.279599e-03 7.977507e-03 displacement 1.662545e+04 1.406147e+02 4.343341e-03 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 12 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 2.828694e+00 0.000000e+00 energy 7.977507e-01 4.197755e-03 7.977507e-03 displacement 1.662545e+04 1.800776e+02 6.270981e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 13 stiffness matrix reformations = 8 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 1.990516e+00 0.000000e+00 energy 7.977507e-01 3.296601e-03 7.977507e-03 displacement 1.662545e+04 1.554929e+02 8.386755e-03 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 14 stiffness matrix reformations = 9 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 1.459127e+00 0.000000e+00 energy 7.977507e-01 2.240336e-03 7.977507e-03 displacement 1.662545e+04 1.351998e+02 1.063931e-02 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 15 stiffness matrix reformations = 10 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 8.174078e-01 0.000000e+00 energy 7.977507e-01 1.697595e-03 7.977507e-03 displacement 1.662545e+04 1.015067e+02 1.280931e-02 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 16 stiffness matrix reformations = 11 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 5.862798e-01 0.000000e+00 energy 7.977507e-01 8.909753e-04 7.977507e-03 displacement 1.662545e+04 8.638792e+01 1.499079e-02 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 12 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 1.901277e-01 0.000000e+00 energy 7.977507e-01 6.285702e-04 7.977507e-03 displacement 1.662545e+04 4.986764e+01 1.676341e-02 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 18 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 1.278718e-01 0.000000e+00 energy 7.977507e-01 1.826646e-04 7.977507e-03 displacement 1.662545e+04 4.046712e+01 1.844531e-02 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 19 stiffness matrix reformations = 14 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 9.435748e-03 0.000000e+00 energy 7.977507e-01 8.257831e-05 7.977507e-03 displacement 1.662545e+04 1.161313e+01 1.937968e-02 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 15 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 1.940556e-03 0.000000e+00 energy 7.977507e-01 3.447355e-06 7.977507e-03 displacement 1.662545e+04 5.127989e+00 2.001336e-02 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 16 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 4.929075e-08 0.000000e+00 energy 7.977507e-01 6.400881e-08 7.977507e-03 displacement 1.662545e+04 3.172596e-02 2.006365e-02 Reforming stiffness matrix: reformation #17 17 Nonlinear solution status: time= 0.543357 stiffness updates = 0 right hand side evaluations = 22 stiffness matrix reformations = 17 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.390910e+00 1.696443e-09 0.000000e+00 energy 7.977507e-01 8.397954e-11 7.977507e-03 displacement 1.662545e+04 6.243648e-04 2.007072e-02 convergence summary number of iterations : 17 number of reformations : 17 ------- converged at time : 0.543357 Data Record #1 =========================================================================== Step = 9 Time = 0.543357478 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 9 Time = 0.543357478 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 9 Time = 0.543357478 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 54%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.108194 ===== beginning time step 10 : 0.651552 ===== Setting parameter "pressure" to : 0.651552 Setting parameter "pressure" to : 0.651552 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 1 step from line search = 0.167155 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 1.342237e+01 0.000000e+00 energy 1.086715e+00 1.240556e-01 1.086715e-02 displacement 1.270095e+04 3.548744e+02 3.548744e-04 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 5.355046e+01 0.000000e+00 energy 1.086715e+00 1.021332e-01 1.086715e-02 displacement 1.270095e+04 1.695902e+03 5.675129e-04 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 1.002993e+00 0.000000e+00 energy 1.086715e+00 1.343341e-02 1.086715e-02 displacement 1.270095e+04 4.788023e+01 8.543225e-04 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 4 step from line search = 0.262470 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 3.113955e+00 0.000000e+00 energy 1.086715e+00 6.822554e-03 1.086715e-02 displacement 1.270095e+04 1.366864e+02 1.643760e-03 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 2.246021e+01 0.000000e+00 energy 1.086715e+00 2.223809e-02 1.086715e-02 displacement 1.270095e+04 5.590171e+02 4.067765e-03 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 2.652969e-01 0.000000e+00 energy 1.086715e+00 7.365273e-03 1.086715e-02 displacement 1.270095e+04 5.592465e+01 5.062259e-03 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 11 stiffness matrix reformations = 7 step from line search = 0.212061 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 1.240783e+00 0.000000e+00 energy 1.086715e+00 3.934979e-03 1.086715e-02 displacement 1.270095e+04 7.760721e+01 6.379851e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 13 stiffness matrix reformations = 8 step from line search = 0.532646 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 2.951920e+00 0.000000e+00 energy 1.086715e+00 2.675850e-03 1.086715e-02 displacement 1.270095e+04 1.365493e+02 8.367405e-03 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 14 stiffness matrix reformations = 9 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 2.114648e+00 0.000000e+00 energy 1.086715e+00 3.372095e-03 1.086715e-02 displacement 1.270095e+04 1.697127e+02 1.090414e-02 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 15 stiffness matrix reformations = 10 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 1.429248e+00 0.000000e+00 energy 1.086715e+00 2.499722e-03 1.086715e-02 displacement 1.270095e+04 1.345800e+02 1.344710e-02 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 16 stiffness matrix reformations = 11 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 9.425885e-01 0.000000e+00 energy 1.086715e+00 1.604204e-03 1.086715e-02 displacement 1.270095e+04 1.117521e+02 1.599831e-02 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 12 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 4.779206e-01 0.000000e+00 energy 1.086715e+00 1.072828e-03 1.086715e-02 displacement 1.270095e+04 7.830698e+01 1.830478e-02 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 18 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 2.652061e-01 0.000000e+00 energy 1.086715e+00 5.084133e-04 1.086715e-02 displacement 1.270095e+04 5.890987e+01 2.043179e-02 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 19 stiffness matrix reformations = 14 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 6.524229e-02 0.000000e+00 energy 1.086715e+00 2.526927e-04 1.086715e-02 displacement 1.270095e+04 2.959194e+01 2.201041e-02 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 15 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 1.785696e-02 0.000000e+00 energy 1.086715e+00 5.157089e-05 1.086715e-02 displacement 1.270095e+04 1.551608e+01 2.319036e-02 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 16 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 2.647184e-04 0.000000e+00 energy 1.086715e+00 4.895939e-06 1.086715e-02 displacement 1.270095e+04 2.103321e+00 2.363285e-02 Reforming stiffness matrix: reformation #17 17 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 22 stiffness matrix reformations = 17 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 3.737505e-07 0.000000e+00 energy 1.086715e+00 3.885521e-08 1.086715e-02 displacement 1.270095e+04 6.244630e-02 2.370960e-02 Reforming stiffness matrix: reformation #18 18 Nonlinear solution status: time= 0.651552 stiffness updates = 0 right hand side evaluations = 23 stiffness matrix reformations = 18 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.824005e+00 6.432104e-09 0.000000e+00 energy 1.086715e+00 2.221703e-09 1.086715e-02 displacement 1.270095e+04 1.499998e-04 2.370603e-02 convergence summary number of iterations : 18 number of reformations : 18 ------- converged at time : 0.651552 Data Record #1 =========================================================================== Step = 10 Time = 0.651551588 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 10 Time = 0.651551588 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 10 Time = 0.651551588 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 65%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.117442 ===== beginning time step 11 : 0.768994 ===== Setting parameter "pressure" to : 0.768994 Setting parameter "pressure" to : 0.768994 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 1 step from line search = 0.208196 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 2.242172e+01 0.000000e+00 energy 1.333962e+00 1.707785e-01 1.333962e-02 displacement 1.052806e+04 4.563460e+02 4.563460e-04 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 3.132591e+01 0.000000e+00 energy 1.333962e+00 4.932922e-02 1.333962e-02 displacement 1.052806e+04 1.346844e+03 3.337730e-04 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 6.937501e+00 0.000000e+00 energy 1.333962e+00 1.261707e-02 1.333962e-02 displacement 1.052806e+04 1.392412e+02 7.310924e-04 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 1.452990e+01 0.000000e+00 energy 1.333962e+00 9.803829e-03 1.333962e-02 displacement 1.052806e+04 7.012481e+02 2.804713e-03 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.768994 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.190411e+00 2.056246e+00 0.000000e+00 energy 1.333962e+00 1.116958e-02 1.333962e-02 displacement 1.052806e+04 8.193817e+01 3.782437e-03 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 6 step from line search = 0.453755 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 4.043949e+00 0.000000e+00 energy 1.333962e+00 4.528767e-03 1.333962e-02 displacement 1.052806e+04 1.804384e+02 5.596844e-03 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 10 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 5.285255e+00 0.000000e+00 energy 1.333962e+00 3.413596e-03 1.333962e-02 displacement 1.052806e+04 2.477284e+02 8.168940e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 11 stiffness matrix reformations = 8 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 1.427669e+00 0.000000e+00 energy 1.333962e+00 5.681730e-03 1.333962e-02 displacement 1.052806e+04 1.404090e+02 1.043309e-02 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 12 stiffness matrix reformations = 9 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 5.839223e+00 0.000000e+00 energy 1.333962e+00 3.374728e-03 1.333962e-02 displacement 1.052806e+04 2.566362e+02 1.393618e-02 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 13 stiffness matrix reformations = 10 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 1.359805e-01 0.000000e+00 energy 1.333962e+00 2.666518e-03 1.333962e-02 displacement 1.052806e+04 4.734644e+01 1.559955e-02 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 15 stiffness matrix reformations = 11 step from line search = 0.371910 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 6.087826e-01 0.000000e+00 energy 1.333962e+00 9.711361e-04 1.333962e-02 displacement 1.052806e+04 5.672110e+01 1.752561e-02 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 16 stiffness matrix reformations = 12 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 8.577332e-01 0.000000e+00 energy 1.333962e+00 2.884295e-04 1.333962e-02 displacement 1.052806e+04 1.067245e+02 2.035392e-02 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 6.110281e-02 0.000000e+00 energy 1.333962e+00 6.064457e-04 1.333962e-02 displacement 1.052806e+04 2.960106e+01 2.192896e-02 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 18 stiffness matrix reformations = 14 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 2.613147e-01 0.000000e+00 energy 1.333962e+00 2.603667e-04 1.333962e-02 displacement 1.052806e+04 5.590507e+01 2.418867e-02 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 19 stiffness matrix reformations = 15 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 1.162700e-04 0.000000e+00 energy 1.333962e+00 2.264056e-05 1.333962e-02 displacement 1.052806e+04 1.878036e+00 2.461569e-02 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 16 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 5.069946e-04 0.000000e+00 energy 1.333962e+00 5.317144e-06 1.333962e-02 displacement 1.052806e+04 2.503073e+00 2.511259e-02 Reforming stiffness matrix: reformation #17 17 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 17 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 9.112555e-06 0.000000e+00 energy 1.333962e+00 8.606956e-07 1.333962e-02 displacement 1.052806e+04 1.636786e-01 2.498585e-02 Reforming stiffness matrix: reformation #18 18 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 22 stiffness matrix reformations = 18 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 9.791833e-07 0.000000e+00 energy 1.333962e+00 1.829960e-07 1.333962e-02 displacement 1.052806e+04 4.213391e-02 2.505033e-02 Reforming stiffness matrix: reformation #19 19 Nonlinear solution status: time= 0.768994 stiffness updates = 0 right hand side evaluations = 23 stiffness matrix reformations = 19 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.190411e+00 1.761706e-07 0.000000e+00 energy 1.333962e+00 4.199515e-08 1.333962e-02 displacement 1.052806e+04 9.256478e-03 2.502011e-02 convergence summary number of iterations : 19 number of reformations : 19 ------- converged at time : 0.768994 Data Record #1 =========================================================================== Step = 11 Time = 0.768993636 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 11 Time = 0.768993636 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 11 Time = 0.768993636 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 77%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.123701 ===== beginning time step 12 : 0.892695 ===== Setting parameter "pressure" to : 0.892695 Setting parameter "pressure" to : 0.892695 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 1 step from line search = 0.243920 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 3.238200e+01 0.000000e+00 energy 1.537509e+00 2.093864e-01 1.537509e-02 displacement 9.031469e+03 5.373452e+02 5.373452e-04 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 2.200540e+01 0.000000e+00 energy 1.537509e+00 2.746265e-02 1.537509e-02 displacement 9.031469e+03 1.157064e+03 2.446710e-04 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 2.570848e+01 0.000000e+00 energy 1.537509e+00 1.449558e-02 1.537509e-02 displacement 9.031469e+03 2.759621e+02 7.388736e-04 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 2.312679e+00 0.000000e+00 energy 1.537509e+00 1.160505e-02 1.537509e-02 displacement 9.031469e+03 3.712193e+02 2.122900e-03 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 5 step from line search = 0.373700 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 5.296833e+00 0.000000e+00 energy 1.537509e+00 6.598057e-03 1.537509e-02 displacement 9.031469e+03 9.014000e+01 3.025448e-03 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 9.085888e+00 0.000000e+00 energy 1.537509e+00 1.595632e-03 1.537509e-02 displacement 9.031469e+03 3.258936e+02 5.285278e-03 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 10 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 1.352331e+00 0.000000e+00 energy 1.537509e+00 8.679258e-03 1.537509e-02 displacement 9.031469e+03 1.691869e+02 7.330926e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 12 stiffness matrix reformations = 8 step from line search = 0.543115 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 3.375854e+00 0.000000e+00 energy 1.537509e+00 2.866172e-03 1.537509e-02 displacement 9.031469e+03 1.336143e+02 9.418590e-03 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 0.892695 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.482976e+00 2.044387e+00 0.000000e+00 energy 1.537509e+00 3.843432e-03 1.537509e-02 displacement 9.031469e+03 1.796212e+02 1.218042e-02 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 0.892695 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.482976e+00 1.672716e+00 0.000000e+00 energy 1.537509e+00 2.362046e-03 1.537509e-02 displacement 9.031469e+03 1.439526e+02 1.495315e-02 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 0.892695 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.482976e+00 7.155105e-01 0.000000e+00 energy 1.537509e+00 1.854472e-03 1.537509e-02 displacement 9.031469e+03 1.026728e+02 1.751973e-02 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 16 stiffness matrix reformations = 12 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 6.282703e-01 0.000000e+00 energy 1.537509e+00 7.055925e-04 1.537509e-02 displacement 9.031469e+03 8.736679e+01 2.006646e-02 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 17 stiffness matrix reformations = 13 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 9.938465e-02 0.000000e+00 energy 1.537509e+00 5.556899e-04 1.537509e-02 displacement 9.031469e+03 3.847531e+01 2.185381e-02 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 18 stiffness matrix reformations = 14 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 1.030611e-01 0.000000e+00 energy 1.537509e+00 3.725630e-05 1.537509e-02 displacement 9.031469e+03 3.552013e+01 2.364206e-02 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 19 stiffness matrix reformations = 15 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 4.860539e-04 0.000000e+00 energy 1.537509e+00 2.551247e-05 1.537509e-02 displacement 9.031469e+03 3.274742e+00 2.420006e-02 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 16 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 1.187891e-04 0.000000e+00 energy 1.537509e+00 1.368184e-06 1.537509e-02 displacement 9.031469e+03 1.256948e+00 2.454872e-02 Reforming stiffness matrix: reformation #17 17 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 21 stiffness matrix reformations = 17 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 3.701989e-06 0.000000e+00 energy 1.537509e+00 6.522035e-07 1.537509e-02 displacement 9.031469e+03 7.011435e-02 2.446704e-02 Reforming stiffness matrix: reformation #18 18 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 22 stiffness matrix reformations = 18 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 1.332825e-06 0.000000e+00 energy 1.537509e+00 2.434649e-07 1.537509e-02 displacement 9.031469e+03 4.620406e-02 2.453382e-02 Reforming stiffness matrix: reformation #19 19 Nonlinear solution status: time= 0.892695 stiffness updates = 0 right hand side evaluations = 23 stiffness matrix reformations = 19 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.482976e+00 4.038641e-07 0.000000e+00 energy 1.537509e+00 9.061761e-08 1.537509e-02 displacement 9.031469e+03 1.513937e-02 2.449562e-02 convergence summary number of iterations : 19 number of reformations : 19 ------- converged at time : 0.892695 Data Record #1 =========================================================================== Step = 12 Time = 0.89269505 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 12 Time = 0.89269505 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 12 Time = 0.89269505 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 89%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.12904 MUST POINT CONTROLLER: adjusting time step. dt = 0.107305 ===== beginning time step 13 : 1 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 0 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 1 step from line search = 0.206832 convergence norms : INITIAL CURRENT REQUIRED residual 1.912925e+00 2.004027e+01 0.000000e+00 energy 1.172354e+00 1.489090e-01 1.172354e-02 displacement 9.776428e+03 4.182301e+02 4.182301e-04 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.912925e+00 2.741698e+01 0.000000e+00 energy 1.172354e+00 4.067789e-02 1.172354e-02 displacement 9.776428e+03 1.260415e+03 3.218974e-04 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.912925e+00 5.284825e+00 0.000000e+00 energy 1.172354e+00 1.167741e-02 1.172354e-02 displacement 9.776428e+03 1.171493e+02 6.664263e-04 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.912925e+00 1.425518e+01 0.000000e+00 energy 1.172354e+00 1.069149e-02 1.172354e-02 displacement 9.776428e+03 6.231761e+02 2.517387e-03 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.912925e+00 9.931017e-01 0.000000e+00 energy 1.172354e+00 8.710093e-03 1.172354e-02 displacement 9.776428e+03 6.007670e+01 3.306654e-03 Reforming stiffness matrix: reformation #6 6 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 6 step from line search = 0.373875 convergence norms : INITIAL CURRENT REQUIRED residual 1.912925e+00 2.570862e+00 0.000000e+00 energy 1.172354e+00 3.486886e-03 1.172354e-02 displacement 9.776428e+03 1.235333e+02 4.690372e-03 Reforming stiffness matrix: reformation #7 7 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 10 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.912925e+00 4.751565e+00 0.000000e+00 energy 1.172354e+00 3.715718e-04 1.172354e-02 displacement 9.776428e+03 2.437117e+02 7.045205e-03 Reforming stiffness matrix: reformation #8 8 Nonlinear solution status: time= 1 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.912925e+00 4.571057e-01 0.000000e+00 energy 1.172354e+00 3.776770e-03 1.172354e-02 displacement 9.776428e+03 7.513391e+01 8.561179e-03 Reforming stiffness matrix: reformation #9 9 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 13 stiffness matrix reformations = 9 step from line search = 0.536651 convergence norms : INITIAL CURRENT REQUIRED residual 1.912925e+00 1.222922e+00 0.000000e+00 energy 1.172354e+00 1.004883e-03 1.172354e-02 displacement 9.776428e+03 8.301118e+01 1.031570e-02 Reforming stiffness matrix: reformation #10 10 Nonlinear solution status: time= 1 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.912925e+00 4.928355e-01 0.000000e+00 energy 1.172354e+00 1.307345e-03 1.172354e-02 displacement 9.776428e+03 8.197105e+01 1.222499e-02 Reforming stiffness matrix: reformation #11 11 Nonlinear solution status: time= 1 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.912925e+00 3.179401e-01 0.000000e+00 energy 1.172354e+00 4.886606e-04 1.172354e-02 displacement 9.776428e+03 6.208281e+01 1.401852e-02 Reforming stiffness matrix: reformation #12 12 Nonlinear solution status: time= 1 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.912925e+00 4.051973e-02 0.000000e+00 energy 1.172354e+00 2.546739e-04 1.172354e-02 displacement 9.776428e+03 2.435568e+01 1.520580e-02 Reforming stiffness matrix: reformation #13 13 Nonlinear solution status: time= 1 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.912925e+00 1.384898e-02 0.000000e+00 energy 1.172354e+00 2.292116e-05 1.172354e-02 displacement 9.776428e+03 1.353395e+01 1.612213e-02 Reforming stiffness matrix: reformation #14 14 Nonlinear solution status: time= 1 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.912925e+00 2.362979e-06 0.000000e+00 energy 1.172354e+00 2.406753e-06 1.172354e-02 displacement 9.776428e+03 2.764919e-01 1.625558e-02 Reforming stiffness matrix: reformation #15 15 Nonlinear solution status: time= 1 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.912925e+00 2.590307e-06 0.000000e+00 energy 1.172354e+00 1.638271e-07 1.172354e-02 displacement 9.776428e+03 1.126046e-01 1.634095e-02 Reforming stiffness matrix: reformation #16 16 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 20 stiffness matrix reformations = 16 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.912925e+00 5.717863e-07 0.000000e+00 energy 1.172354e+00 1.460608e-07 1.172354e-02 displacement 9.776428e+03 1.389014e-02 1.631122e-02 convergence summary number of iterations : 16 number of reformations : 16 ------- converged at time : 1 Data Record #1 =========================================================================== Step = 13 Time = 1 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 13 Time = 1 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 13 Time = 1 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 1/4: 100%) tempModel.feb - FEBio 3.2.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 .................... : 13 Total number of equilibrium iterations ............ : 250 Average number of equilibrium iterations .......... : 19.2308 Total number of right hand evaluations ............ : 328 Total number of stiffness reformations ............ : 250 L I N E A R S O L V E R S T A T S Total calls to linear solver ........ : 254 Avg iterations per solve ............ : 1 ===== beginning time step 1 : 1.04 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 0.08 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 ===== reforming stiffness matrix: Nr of equations ........................... : 10013 Nr of nonzeroes in stiffness matrix ....... : 277827 1 Nonlinear solution status: time= 1.04 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.650221e+05 2.383652e-03 0.000000e+00 energy 3.954474e+01 2.697325e-03 3.954474e-01 displacement 7.817614e+00 7.817614e+00 7.817614e-06 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1.04 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.650221e+05 3.084955e-06 0.000000e+00 energy 3.954474e+01 6.200294e-07 3.954474e-01 displacement 7.817614e+00 1.977225e-02 8.266001e-06 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1.04 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.650221e+05 5.930563e-08 0.000000e+00 energy 3.954474e+01 1.626846e-08 3.954474e-01 displacement 7.817614e+00 5.710344e-04 8.198825e-06 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1.04 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.650221e+05 1.860370e-09 0.000000e+00 energy 3.954474e+01 4.652857e-10 3.954474e-01 displacement 7.817614e+00 1.550323e-05 8.211826e-06 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 1.04 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.650221e+05 5.794170e-11 0.000000e+00 energy 3.954474e+01 1.353197e-11 3.954474e-01 displacement 7.817614e+00 4.840579e-07 8.209539e-06 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 1.04 Data Record #1 =========================================================================== Step = 1 Time = 1.04 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 1 Time = 1.04 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 1 Time = 1.04 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 4%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.064 ===== beginning time step 2 : 1.104 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 0.208 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 1.104 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.228896e+05 1.542497e-02 0.000000e+00 energy 1.009968e+02 1.120613e-02 1.009968e+00 displacement 1.989533e+01 1.989533e+01 1.989533e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1.104 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.228896e+05 8.437036e-06 0.000000e+00 energy 1.009968e+02 1.168140e-06 1.009968e+00 displacement 1.989533e+01 4.781786e-02 2.098877e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1.104 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.228896e+05 8.368905e-08 0.000000e+00 energy 1.009968e+02 2.704225e-08 1.009968e+00 displacement 1.989533e+01 8.482596e-04 2.087092e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1.104 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.228896e+05 2.855099e-09 0.000000e+00 energy 1.009968e+02 7.452581e-10 1.009968e+00 displacement 1.989533e+01 2.505002e-05 2.089662e-05 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 1.104 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.228896e+05 8.985176e-11 0.000000e+00 energy 1.009968e+02 2.140924e-11 1.009968e+00 displacement 1.989533e+01 7.647479e-07 2.089209e-05 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 1.104 Data Record #1 =========================================================================== Step = 2 Time = 1.104 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 2 Time = 1.104 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 2 Time = 1.104 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 10%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.0832 ===== beginning time step 3 : 1.1872 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 0.3744 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 1.1872 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.158409e+05 4.351970e-02 0.000000e+00 energy 1.700354e+02 2.449758e-02 1.700354e+00 displacement 3.350601e+01 3.350601e+01 3.350601e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1.1872 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.158409e+05 1.450183e-05 0.000000e+00 energy 1.700354e+02 1.214395e-06 1.700354e+00 displacement 3.350601e+01 7.390587e-02 3.525170e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1.1872 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.158409e+05 5.547077e-08 0.000000e+00 energy 1.700354e+02 2.541938e-08 1.700354e+00 displacement 3.350601e+01 6.473771e-04 3.514914e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1.1872 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.158409e+05 2.254855e-09 0.000000e+00 energy 1.700354e+02 6.448566e-10 1.700354e+00 displacement 3.350601e+01 2.181180e-05 3.517895e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 1.1872 Data Record #1 =========================================================================== Step = 3 Time = 1.1872 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 3 Time = 1.1872 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 3 Time = 1.1872 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 19%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.09856 ===== beginning time step 4 : 1.28576 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 0.57152 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 1.28576 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.006596e+06 8.512041e-02 0.000000e+00 energy 2.374293e+02 4.051999e-02 2.374293e+00 displacement 4.680094e+01 4.680094e+01 4.680094e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1.28576 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.006596e+06 2.036348e-05 0.000000e+00 energy 2.374293e+02 7.941826e-07 2.374293e+00 displacement 4.680094e+01 9.531021e-02 4.912449e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1.28576 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.006596e+06 1.707899e-08 0.000000e+00 energy 2.374293e+02 1.843421e-08 2.374293e+00 displacement 4.680094e+01 3.351708e-04 4.910330e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1.28576 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.006596e+06 1.072784e-09 0.000000e+00 energy 2.374293e+02 3.840406e-10 2.374293e+00 displacement 4.680094e+01 1.280230e-05 4.912784e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 1.28576 Data Record #1 =========================================================================== Step = 4 Time = 1.28576 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 4 Time = 1.28576 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 4 Time = 1.28576 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 29%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.110848 ===== beginning time step 5 : 1.39661 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 0.793216 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 1.39661 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.276186e+06 1.352351e-01 0.000000e+00 energy 2.985490e+02 5.724949e-02 2.985490e+00 displacement 5.888809e+01 5.888809e+01 5.888809e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1.39661 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.276186e+06 2.552389e-05 0.000000e+00 energy 2.985490e+02 8.628850e-08 2.985490e+00 displacement 5.888809e+01 1.116626e-01 6.169343e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1.39661 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.276186e+06 1.679806e-09 0.000000e+00 energy 2.985490e+02 1.286684e-08 2.985490e+00 displacement 5.888809e+01 2.164856e-04 6.179766e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1.39661 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.276186e+06 2.148955e-10 0.000000e+00 energy 2.985490e+02 1.710713e-10 2.985490e+00 displacement 5.888809e+01 4.911587e-06 6.181019e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 1.39661 Data Record #1 =========================================================================== Step = 5 Time = 1.396608 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 5 Time = 1.396608 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 5 Time = 1.396608 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 40%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.120678 ===== beginning time step 6 : 1.51729 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 1.03457 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 1.51729 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.516325e+06 1.886806e-01 0.000000e+00 energy 3.514856e+02 7.331109e-02 3.514856e+00 displacement 6.938785e+01 6.938785e+01 6.938785e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1.51729 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.516325e+06 2.989146e-05 0.000000e+00 energy 3.514856e+02 7.427626e-07 3.514856e+00 displacement 6.938785e+01 1.238455e-01 7.258473e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1.51729 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.516325e+06 1.888319e-08 0.000000e+00 energy 3.514856e+02 1.155789e-08 3.514856e+00 displacement 6.938785e+01 3.883186e-04 7.283496e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1.51729 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.516325e+06 2.753058e-11 0.000000e+00 energy 3.514856e+02 9.241807e-11 3.514856e+00 displacement 6.938785e+01 1.192745e-06 7.283187e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 1.51729 Data Record #1 =========================================================================== Step = 6 Time = 1.5172864 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 6 Time = 1.5172864 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 6 Time = 1.5172864 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 52%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.128543 ===== beginning time step 7 : 1.64583 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 1.29166 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 1.64583 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.724775e+06 2.412205e-01 0.000000e+00 energy 3.958702e+02 8.787178e-02 3.958702e+00 displacement 7.822910e+01 7.822910e+01 7.822910e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1.64583 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.724775e+06 3.350643e-05 0.000000e+00 energy 3.958702e+02 1.579474e-06 3.958702e+00 displacement 7.822910e+01 1.328065e-01 8.174088e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1.64583 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.724775e+06 6.380236e-08 0.000000e+00 energy 3.958702e+02 1.448000e-08 3.958702e+00 displacement 7.822910e+01 8.103078e-04 8.213969e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1.64583 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.724775e+06 4.711702e-10 0.000000e+00 energy 3.958702e+02 1.471082e-10 3.958702e+00 displacement 7.822910e+01 1.842372e-06 8.211983e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 1.64583 Data Record #1 =========================================================================== Step = 7 Time = 1.64582912 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 7 Time = 1.64582912 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 7 Time = 1.64582912 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 65%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.134834 ===== beginning time step 8 : 1.78066 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 1.56133 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 1.78066 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.902562e+06 2.900016e-01 0.000000e+00 energy 4.321504e+02 1.005351e-01 4.321504e+00 displacement 8.550060e+01 8.550060e+01 8.550060e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1.78066 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.902562e+06 3.645692e-05 0.000000e+00 energy 4.321504e+02 2.359464e-06 4.321504e+00 displacement 8.550060e+01 1.393574e-01 8.926513e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1.78066 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.902562e+06 1.264257e-07 0.000000e+00 energy 4.321504e+02 2.038539e-08 4.321504e+00 displacement 8.550060e+01 1.391269e-03 8.980364e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1.78066 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.902562e+06 1.350676e-09 0.000000e+00 energy 4.321504e+02 2.971589e-10 4.321504e+00 displacement 8.550060e+01 5.779685e-06 8.976748e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 1.78066 Data Record #1 =========================================================================== Step = 8 Time = 1.7806633 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 8 Time = 1.7806633 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 8 Time = 1.7806633 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 78%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.139867 ===== beginning time step 9 : 1.92053 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 1.84106 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 1.92053 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.052317e+06 3.334301e-01 0.000000e+00 energy 4.611584e+02 1.112001e-01 4.611584e+00 displacement 9.136608e+01 9.136608e+01 9.136608e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 1.92053 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.052317e+06 3.883819e-05 0.000000e+00 energy 4.611584e+02 3.052593e-06 4.611584e+00 displacement 9.136608e+01 1.441260e-01 9.533390e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 1.92053 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.052317e+06 1.968664e-07 0.000000e+00 energy 4.611584e+02 2.784741e-08 4.611584e+00 displacement 9.136608e+01 2.040913e-03 9.599711e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 1.92053 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.052317e+06 2.456011e-09 0.000000e+00 energy 4.611584e+02 4.992304e-10 4.611584e+00 displacement 9.136608e+01 1.166973e-05 9.594607e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 1.92053 Data Record #1 =========================================================================== Step = 9 Time = 1.92053064 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 9 Time = 1.92053064 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 9 Time = 1.92053064 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 92%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.143894 MUST POINT CONTROLLER: adjusting time step. dt = 0.0794694 ===== beginning time step 10 : 2 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -0 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 2 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.641035e+05 3.477714e-02 0.000000e+00 energy 1.475826e+02 2.019942e-02 1.475826e+00 displacement 2.928770e+01 2.928770e+01 2.928770e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.641035e+05 1.255610e-05 0.000000e+00 energy 1.475826e+02 9.957315e-07 1.475826e+00 displacement 2.928770e+01 6.570871e-02 3.089985e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.641035e+05 5.543679e-08 0.000000e+00 energy 1.475826e+02 2.045577e-08 1.475826e+00 displacement 2.928770e+01 5.620318e-04 3.079052e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.641035e+05 2.023682e-09 0.000000e+00 energy 1.475826e+02 5.151612e-10 1.475826e+00 displacement 2.928770e+01 1.875942e-05 3.081805e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 2 Data Record #1 =========================================================================== Step = 10 Time = 2 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 10 Time = 2 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 10 Time = 2 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 2/4: 100%) tempModel.feb - FEBio 3.2.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 .................... : 10 Total number of equilibrium iterations ............ : 42 Average number of equilibrium iterations .......... : 4.2 Total number of right hand evaluations ............ : 52 Total number of stiffness reformations ............ : 42 L I N E A R S O L V E R S T A T S Total calls to linear solver ........ : 42 Avg iterations per solve ............ : 1 ===== beginning time step 1 : 2.04 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -0.08 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 ===== reforming stiffness matrix: Nr of equations ........................... : 10013 Nr of nonzeroes in stiffness matrix ....... : 277827 1 Nonlinear solution status: time= 2.04 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.684605e+05 2.319094e-03 0.000000e+00 energy 3.720447e+01 2.515286e-03 3.720447e-01 displacement 7.391237e+00 7.391237e+00 7.391237e-06 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2.04 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.684605e+05 2.788612e-06 0.000000e+00 energy 3.720447e+01 4.860362e-07 3.720447e-01 displacement 7.391237e+00 1.932350e-02 7.853448e-06 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2.04 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.684605e+05 5.099893e-08 0.000000e+00 energy 3.720447e+01 1.258767e-08 3.720447e-01 displacement 7.391237e+00 4.617729e-04 7.790113e-06 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2.04 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.684605e+05 1.495991e-09 0.000000e+00 energy 3.720447e+01 3.462545e-10 3.720447e-01 displacement 7.391237e+00 1.230468e-05 7.801950e-06 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 2.04 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.684605e+05 4.393173e-11 0.000000e+00 energy 3.720447e+01 9.643643e-12 3.720447e-01 displacement 7.391237e+00 3.620968e-07 7.799931e-06 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 2.04 Data Record #1 =========================================================================== Step = 1 Time = 2.04 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 1 Time = 2.04 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 1 Time = 2.04 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 4%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.064 ===== beginning time step 2 : 2.104 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -0.208 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 2.104 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.309918e+05 1.468558e-02 0.000000e+00 energy 9.548254e+01 1.037115e-02 9.548254e-01 displacement 1.895478e+01 1.895478e+01 1.895478e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2.104 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.309918e+05 7.204117e-06 0.000000e+00 energy 9.548254e+01 8.183350e-07 9.548254e-01 displacement 1.895478e+01 4.639651e-02 2.010527e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2.104 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.309918e+05 7.139130e-08 0.000000e+00 energy 9.548254e+01 2.002383e-08 9.548254e-01 displacement 1.895478e+01 6.657282e-04 1.999430e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2.104 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.309918e+05 2.248674e-09 0.000000e+00 energy 9.548254e+01 5.406592e-10 9.548254e-01 displacement 1.895478e+01 1.923188e-05 2.001757e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 2.104 Data Record #1 =========================================================================== Step = 2 Time = 2.104 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 2 Time = 2.104 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 2 Time = 2.104 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 10%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.0832 ===== beginning time step 3 : 2.1872 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -0.3744 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 2.1872 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.276286e+05 4.167580e-02 0.000000e+00 energy 1.620122e+02 2.291744e-02 1.620122e+00 displacement 3.213229e+01 3.213229e+01 3.213229e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2.1872 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.276286e+05 1.208736e-05 0.000000e+00 energy 1.620122e+02 6.825786e-07 1.620122e+00 displacement 3.213229e+01 7.307097e-02 3.401391e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2.1872 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.276286e+05 4.766464e-08 0.000000e+00 energy 1.620122e+02 1.798158e-08 1.620122e+00 displacement 3.213229e+01 4.968223e-04 3.391704e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2.1872 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.276286e+05 1.765954e-09 0.000000e+00 energy 1.620122e+02 4.610760e-10 1.620122e+00 displacement 3.213229e+01 1.636785e-05 3.394398e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 2.1872 Data Record #1 =========================================================================== Step = 3 Time = 2.1872 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 3 Time = 2.1872 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 3 Time = 2.1872 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 19%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.09856 ===== beginning time step 4 : 2.28576 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -0.57152 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 2.28576 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.019663e+06 8.200603e-02 0.000000e+00 energy 2.285349e+02 3.830954e-02 2.285349e+00 displacement 4.528559e+01 4.528559e+01 4.528559e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2.28576 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.019663e+06 1.663142e-05 0.000000e+00 energy 2.285349e+02 1.201683e-07 2.285349e+00 displacement 4.528559e+01 9.586153e-02 4.784305e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2.28576 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.019663e+06 1.448015e-08 0.000000e+00 energy 2.285349e+02 1.202240e-08 2.285349e+00 displacement 4.528559e+01 2.519512e-04 4.782373e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2.28576 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.019663e+06 8.167341e-10 0.000000e+00 energy 2.285349e+02 2.642419e-10 2.285349e+00 displacement 4.528559e+01 9.124837e-06 4.784569e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 2.28576 Data Record #1 =========================================================================== Step = 4 Time = 2.28576 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 4 Time = 2.28576 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 4 Time = 2.28576 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 29%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.110848 ===== beginning time step 5 : 2.39661 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -0.793216 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 2.39661 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.287507e+06 1.314768e-01 0.000000e+00 energy 2.908431e+02 5.484628e-02 2.908431e+00 displacement 5.757554e+01 5.757554e+01 5.757554e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2.39661 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.287507e+06 2.054948e-05 0.000000e+00 energy 2.908431e+02 7.027944e-07 2.908431e+00 displacement 5.757554e+01 1.143012e-01 6.071940e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2.39661 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.287507e+06 1.400979e-09 0.000000e+00 energy 2.908431e+02 7.731719e-09 2.908431e+00 displacement 5.757554e+01 1.844227e-04 6.082244e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2.39661 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.287507e+06 1.370191e-10 0.000000e+00 energy 2.908431e+02 1.085982e-10 2.908431e+00 displacement 5.757554e+01 3.013480e-06 6.083300e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 2.39661 Data Record #1 =========================================================================== Step = 5 Time = 2.396608 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 5 Time = 2.396608 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 5 Time = 2.396608 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 40%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.120678 ===== beginning time step 6 : 2.51729 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -1.03457 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 2.51729 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.522800e+06 1.853806e-01 0.000000e+00 energy 3.470778e+02 7.127533e-02 3.470778e+00 displacement 6.863729e+01 6.863729e+01 6.863729e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2.51729 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.522800e+06 2.382309e-05 0.000000e+00 energy 3.470778e+02 1.632622e-06 3.470778e+00 displacement 6.863729e+01 1.289131e-01 7.227339e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2.51729 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.522800e+06 1.895230e-08 0.000000e+00 energy 3.470778e+02 7.609781e-09 3.470778e+00 displacement 6.863729e+01 3.812354e-04 7.252252e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2.51729 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.522800e+06 4.418011e-11 0.000000e+00 energy 3.470778e+02 6.672989e-11 3.470778e+00 displacement 6.863729e+01 5.690673e-07 7.251801e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 2.51729 Data Record #1 =========================================================================== Step = 6 Time = 2.5172864 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 6 Time = 2.5172864 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 6 Time = 2.5172864 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 52%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.128543 ===== beginning time step 7 : 2.64583 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -1.29166 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 2.64583 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.723547e+06 2.397874e-01 0.000000e+00 energy 3.967063e+02 8.680877e-02 3.967063e+00 displacement 7.837187e+01 7.837187e+01 7.837187e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2.64583 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.723547e+06 2.653538e-05 0.000000e+00 energy 3.967063e+02 2.563714e-06 3.967063e+00 displacement 7.837187e+01 1.404499e-01 8.241360e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2.64583 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.723547e+06 6.479842e-08 0.000000e+00 energy 3.967063e+02 1.178792e-08 3.967063e+00 displacement 7.837187e+01 8.143724e-04 8.281584e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2.64583 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.723547e+06 5.398173e-10 0.000000e+00 energy 3.967063e+02 1.444066e-10 3.967063e+00 displacement 7.837187e+01 2.023856e-06 8.279481e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 2.64583 Data Record #1 =========================================================================== Step = 7 Time = 2.64582912 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 7 Time = 2.64582912 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 7 Time = 2.64582912 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 65%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.134834 ===== beginning time step 8 : 2.78066 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -1.56133 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 2.78066 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.891161e+06 2.919335e-01 0.000000e+00 energy 4.399114e+02 1.010249e-01 4.399114e+00 displacement 8.682323e+01 8.682323e+01 8.682323e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2.78066 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.891161e+06 2.878685e-05 0.000000e+00 energy 4.399114e+02 3.434979e-06 4.399114e+00 displacement 8.682323e+01 1.496081e-01 9.119556e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2.78066 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.891161e+06 1.313536e-07 0.000000e+00 energy 4.399114e+02 1.929537e-08 4.399114e+00 displacement 8.682323e+01 1.409804e-03 9.174683e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2.78066 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.891161e+06 1.482799e-09 0.000000e+00 energy 4.399114e+02 3.158137e-10 4.399114e+00 displacement 8.682323e+01 6.524323e-06 9.170936e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 2.78066 Data Record #1 =========================================================================== Step = 8 Time = 2.7806633 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 8 Time = 2.7806633 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 8 Time = 2.7806633 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 78%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.139867 ===== beginning time step 9 : 2.92053 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -1.84106 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 2.92053 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.028705e+06 3.401548e-01 0.000000e+00 energy 4.772304e+02 1.137576e-01 4.772304e+00 displacement 9.410557e+01 9.410557e+01 9.410557e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 2.92053 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.028705e+06 3.066670e-05 0.000000e+00 energy 4.772304e+02 4.216928e-06 4.772304e+00 displacement 9.410557e+01 1.569550e-01 9.874544e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 2.92053 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.028705e+06 2.103262e-07 0.000000e+00 energy 4.772304e+02 2.890969e-08 4.772304e+00 displacement 9.410557e+01 2.090353e-03 9.943530e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 2.92053 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.028705e+06 2.702091e-09 0.000000e+00 energy 4.772304e+02 5.473506e-10 4.772304e+00 displacement 9.410557e+01 1.295100e-05 9.938239e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 2.92053 Data Record #1 =========================================================================== Step = 9 Time = 2.92053064 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 9 Time = 2.92053064 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 9 Time = 2.92053064 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 92%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.143894 MUST POINT CONTROLLER: adjusting time step. dt = 0.0794694 ===== beginning time step 10 : 3 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -0 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 3 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.526824e+05 3.613772e-02 0.000000e+00 energy 1.553560e+02 2.105407e-02 1.553560e+00 displacement 3.061337e+01 3.061337e+01 3.061337e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.526824e+05 1.151543e-05 0.000000e+00 energy 1.553560e+02 9.603184e-07 1.553560e+00 displacement 3.061337e+01 7.066726e-02 3.235432e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.526824e+05 6.363784e-08 0.000000e+00 energy 1.553560e+02 2.459974e-08 1.553560e+00 displacement 3.061337e+01 7.029077e-04 3.224625e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.526824e+05 2.371082e-09 0.000000e+00 energy 1.553560e+02 6.625030e-10 1.553560e+00 displacement 3.061337e+01 2.180164e-05 3.227541e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 3 Data Record #1 =========================================================================== Step = 10 Time = 3 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 10 Time = 3 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 10 Time = 3 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 3/4: 100%) tempModel.feb - FEBio 3.2.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 .................... : 10 Total number of equilibrium iterations ............ : 41 Average number of equilibrium iterations .......... : 4.1 Total number of right hand evaluations ............ : 51 Total number of stiffness reformations ............ : 41 L I N E A R S O L V E R S T A T S Total calls to linear solver ........ : 41 Avg iterations per solve ............ : 1 ===== beginning time step 1 : 3.04 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -0.08 Reforming stiffness matrix: reformation #1 ===== reforming stiffness matrix: Nr of equations ........................... : 10013 Nr of nonzeroes in stiffness matrix ....... : 277827 1 Nonlinear solution status: time= 3.04 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.650206e+05 2.446402e-03 0.000000e+00 energy 3.954559e+01 2.696328e-03 3.954559e-01 displacement 7.786018e+00 7.786018e+00 7.786018e-06 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3.04 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.650206e+05 2.967249e-06 0.000000e+00 energy 3.954559e+01 5.931839e-07 3.954559e-01 displacement 7.786018e+00 2.057694e-02 8.253753e-06 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3.04 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.650206e+05 6.062206e-08 0.000000e+00 energy 3.954559e+01 1.621501e-08 3.954559e-01 displacement 7.786018e+00 5.753941e-04 8.187024e-06 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3.04 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.650206e+05 1.892024e-09 0.000000e+00 energy 3.954559e+01 4.702264e-10 3.954559e-01 displacement 7.786018e+00 1.547367e-05 8.199980e-06 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 3.04 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.650206e+05 1.414870e-10 0.000000e+00 energy 3.954559e+01 2.016441e-11 3.954559e-01 displacement 7.786018e+00 5.117189e-07 8.197700e-06 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 3.04 Data Record #1 =========================================================================== Step = 1 Time = 3.04 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 1 Time = 3.04 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 1 Time = 3.04 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 4%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.064 ===== beginning time step 2 : 3.104 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -0.208 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 3.104 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.220080e+05 1.544526e-02 0.000000e+00 energy 1.014770e+02 1.108500e-02 1.014770e+00 displacement 1.998278e+01 1.998278e+01 1.998278e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3.104 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.220080e+05 7.610130e-06 0.000000e+00 energy 1.014770e+02 1.028911e-06 1.014770e+00 displacement 1.998278e+01 4.909452e-02 2.114456e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3.104 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.220080e+05 8.550960e-08 0.000000e+00 energy 1.014770e+02 2.642920e-08 1.014770e+00 displacement 1.998278e+01 8.521508e-04 2.102881e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3.104 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.220080e+05 5.039956e-09 0.000000e+00 energy 1.014770e+02 9.189084e-10 1.014770e+00 displacement 1.998278e+01 2.505073e-05 2.105432e-05 Reforming stiffness matrix: reformation #5 5 Nonlinear solution status: time= 3.104 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 5 step from line search = 0.559367 convergence norms : INITIAL CURRENT REQUIRED residual 4.220080e+05 1.165562e-09 0.000000e+00 energy 1.014770e+02 8.829884e-11 1.014770e+00 displacement 1.998278e+01 9.146967e-07 2.105179e-05 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 3.104 Data Record #1 =========================================================================== Step = 2 Time = 3.104 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 2 Time = 3.104 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 2 Time = 3.104 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 10%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.0832 ===== beginning time step 3 : 3.1872 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -0.3744 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 3.1872 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.119677e+05 4.383026e-02 0.000000e+00 energy 1.721453e+02 2.448439e-02 1.721453e+00 displacement 3.388744e+01 3.388744e+01 3.388744e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3.1872 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.119677e+05 1.270969e-05 0.000000e+00 energy 1.721453e+02 9.471866e-07 1.721453e+00 displacement 3.388744e+01 7.707032e-02 3.578427e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3.1872 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.119677e+05 5.911142e-08 0.000000e+00 energy 1.721453e+02 2.514418e-08 1.721453e+00 displacement 3.388744e+01 6.811226e-04 3.568570e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3.1872 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.119677e+05 8.437335e-09 0.000000e+00 energy 1.721453e+02 1.148908e-09 1.721453e+00 displacement 3.388744e+01 2.314133e-05 3.571541e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 3.1872 Data Record #1 =========================================================================== Step = 3 Time = 3.1872 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 3 Time = 3.1872 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 3 Time = 3.1872 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 19%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.09856 ===== beginning time step 4 : 3.28576 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -0.57152 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 3.28576 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.968134e+05 8.634848e-02 0.000000e+00 energy 2.427582e+02 4.094891e-02 2.427582e+00 displacement 4.777153e+01 4.777153e+01 4.777153e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3.28576 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.968134e+05 1.746099e-05 0.000000e+00 energy 2.427582e+02 3.887125e-07 2.427582e+00 displacement 4.777153e+01 1.009229e-01 5.034572e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3.28576 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.968134e+05 2.076159e-08 0.000000e+00 energy 2.427582e+02 1.874684e-08 2.427582e+00 displacement 4.777153e+01 3.973788e-04 5.033196e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3.28576 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.968134e+05 9.459555e-09 0.000000e+00 energy 2.427582e+02 1.074090e-09 2.427582e+00 displacement 4.777153e+01 1.508765e-05 5.035660e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 3.28576 Data Record #1 =========================================================================== Step = 4 Time = 3.28576 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 4 Time = 3.28576 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 4 Time = 3.28576 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 29%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.110848 ===== beginning time step 5 : 3.39661 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -0.793216 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 3.39661 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.257297e+06 1.385129e-01 0.000000e+00 energy 3.088387e+02 5.862366e-02 3.088387e+00 displacement 6.076251e+01 6.076251e+01 6.076251e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3.39661 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.257297e+06 2.149734e-05 0.000000e+00 energy 3.088387e+02 4.607204e-07 3.088387e+00 displacement 6.076251e+01 1.201114e-01 6.392202e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3.39661 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.257297e+06 2.546962e-09 0.000000e+00 energy 3.088387e+02 1.371453e-08 3.088387e+00 displacement 6.076251e+01 2.949147e-04 6.404035e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3.39661 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.257297e+06 5.094526e-10 0.000000e+00 energy 3.088387e+02 2.335838e-10 3.088387e+00 displacement 6.076251e+01 5.137969e-06 6.405307e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 3.39661 Data Record #1 =========================================================================== Step = 5 Time = 3.396608 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 5 Time = 3.396608 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 5 Time = 3.396608 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 40%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.120678 ===== beginning time step 6 : 3.51729 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -1.03457 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 3.51729 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.485253e+06 1.954073e-01 0.000000e+00 energy 3.684118e+02 7.618323e-02 3.684118e+00 displacement 7.247224e+01 7.247224e+01 7.247224e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3.51729 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.485253e+06 2.480156e-05 0.000000e+00 energy 3.684118e+02 1.430009e-06 3.684118e+00 displacement 7.247224e+01 1.352832e-01 7.612124e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3.51729 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.485253e+06 2.157880e-08 0.000000e+00 energy 3.684118e+02 1.348188e-08 3.684118e+00 displacement 7.247224e+01 4.788893e-04 7.639646e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3.51729 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.485253e+06 8.556930e-11 0.000000e+00 energy 3.684118e+02 1.354412e-10 3.684118e+00 displacement 7.247224e+01 1.308015e-06 7.639335e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 3.51729 Data Record #1 =========================================================================== Step = 6 Time = 3.5172864 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 6 Time = 3.5172864 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 6 Time = 3.5172864 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 52%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.128543 ===== beginning time step 7 : 3.64583 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -1.29166 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 3.64583 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.678796e+06 2.529207e-01 0.000000e+00 energy 4.209160e+02 9.278659e-02 4.209160e+00 displacement 8.279556e+01 8.279556e+01 8.279556e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3.64583 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.678796e+06 2.759190e-05 0.000000e+00 energy 4.209160e+02 2.403345e-06 4.209160e+00 displacement 8.279556e+01 1.472618e-01 8.684623e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3.64583 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.678796e+06 7.295917e-08 0.000000e+00 energy 4.209160e+02 1.811759e-08 4.209160e+00 displacement 8.279556e+01 9.228324e-04 8.728551e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3.64583 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.678796e+06 5.413758e-10 0.000000e+00 energy 4.209160e+02 1.928307e-10 4.209160e+00 displacement 8.279556e+01 1.887859e-06 8.726503e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 3.64583 Data Record #1 =========================================================================== Step = 7 Time = 3.64582912 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 7 Time = 3.64582912 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 7 Time = 3.64582912 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 65%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.134834 ===== beginning time step 8 : 3.78066 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -1.56133 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 3.78066 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.839405e+06 3.081387e-01 0.000000e+00 energy 4.665520e+02 1.079833e-01 4.665520e+00 displacement 9.177752e+01 9.177752e+01 9.177752e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3.78066 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.839405e+06 2.992564e-05 0.000000e+00 energy 4.665520e+02 3.313188e-06 4.665520e+00 displacement 9.177752e+01 1.567894e-01 9.615400e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3.78066 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.839405e+06 1.480032e-07 0.000000e+00 energy 4.665520e+02 2.656062e-08 4.665520e+00 displacement 9.177752e+01 1.547433e-03 9.675279e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3.78066 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.839405e+06 1.546226e-09 0.000000e+00 energy 4.665520e+02 3.643795e-10 4.665520e+00 displacement 9.177752e+01 5.919582e-06 9.671506e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 3.78066 Data Record #1 =========================================================================== Step = 8 Time = 3.7806633 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 8 Time = 3.7806633 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 8 Time = 3.7806633 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 78%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.139867 ===== beginning time step 9 : 3.92053 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -1.84106 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 3.92053 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.970176e+06 3.593054e-01 0.000000e+00 energy 5.058979e+02 1.215960e-01 5.058979e+00 displacement 9.953709e+01 9.953709e+01 9.953709e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 3.92053 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.970176e+06 3.189473e-05 0.000000e+00 energy 5.058979e+02 4.127050e-06 5.058979e+00 displacement 9.953709e+01 1.644644e-01 1.041757e-04 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 3.92053 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.970176e+06 2.373076e-07 0.000000e+00 energy 5.058979e+02 3.744461e-08 5.058979e+00 displacement 9.953709e+01 2.268867e-03 1.049228e-04 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 3.92053 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.970176e+06 2.919844e-09 0.000000e+00 energy 5.058979e+02 6.126135e-10 5.058979e+00 displacement 9.953709e+01 1.217142e-05 1.048689e-04 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 3.92053 Data Record #1 =========================================================================== Step = 9 Time = 3.92053064 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 9 Time = 3.92053064 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 9 Time = 3.92053064 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 92%) tempModel.feb - FEBio 3.2.0 AUTO STEPPER: increasing time step, dt = 0.143894 MUST POINT CONTROLLER: adjusting time step. dt = 0.0794694 ===== beginning time step 10 : 4 ===== Setting parameter "pressure" to : 1 Setting parameter "pressure" to : 1 Setting parameter "value" to : 2 Setting parameter "value" to : -2 Setting parameter "value" to : -2 Reforming stiffness matrix: reformation #1 1 Nonlinear solution status: time= 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 6.328395e+05 3.819623e-02 0.000000e+00 energy 1.646099e+02 2.251018e-02 1.646099e+00 displacement 3.240127e+01 3.240127e+01 3.240127e-05 Reforming stiffness matrix: reformation #2 2 Nonlinear solution status: time= 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 6.328395e+05 1.220027e-05 0.000000e+00 energy 1.646099e+02 1.302766e-06 1.646099e+00 displacement 3.240127e+01 7.509938e-02 3.415082e-05 Reforming stiffness matrix: reformation #3 3 Nonlinear solution status: time= 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 6.328395e+05 7.849413e-08 0.000000e+00 energy 1.646099e+02 3.480584e-08 1.646099e+00 displacement 3.240127e+01 9.644009e-04 3.404133e-05 Reforming stiffness matrix: reformation #4 4 Nonlinear solution status: time= 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 6.328395e+05 4.057502e-09 0.000000e+00 energy 1.646099e+02 1.061991e-09 1.646099e+00 displacement 3.240127e+01 2.956982e-05 3.407365e-05 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 4 Data Record #1 =========================================================================== Step = 10 Time = 4 Data = ux;uy;uz File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt Data Record #2 =========================================================================== Step = 10 Time = 4 Data = s1 File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt Data Record #3 =========================================================================== Step = 10 Time = 4 Data = Fx File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_force_out.txt ]0;(step 4/4: 100%) tempModel.feb - FEBio 3.2.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 .................... : 10 Total number of equilibrium iterations ............ : 42 Average number of equilibrium iterations .......... : 4.2 Total number of right hand evaluations ............ : 53 Total number of stiffness reformations ............ : 42 L I N E A R S O L V E R S T A T S Total calls to linear solver ........ : 42 Avg iterations per solve ............ : 1 N O N L I N E A R I T E R A T I O N S U M M A R Y Number of time steps completed .................... : 43 Total number of equilibrium iterations ............ : 375 Total number of right hand evaluations ............ : 484 Total number of stiffness reformations ............ : 375 Time in linear solver: 0:00:26 ]0;(step 4/4: 100%) tempModel.feb - FEBio 3.2.0 Elapsed time : 0:00:37 N O R M A L T E R M I N A T I O N * Log file found. 23-Apr-2021 19:44:04 # Parsing log file... 23-Apr-2021 19:44:04 number of iterations : 10 23-Apr-2021 19:44:05 number of reformations : 10 23-Apr-2021 19:44:05 ------- converged at time : 0.04 23-Apr-2021 19:44:05 number of iterations : 28 23-Apr-2021 19:44:05 number of reformations : 28 23-Apr-2021 19:44:05 ------- converged at time : 0.104 23-Apr-2021 19:44:05 number of iterations : 28 23-Apr-2021 19:44:05 number of reformations : 28 23-Apr-2021 19:44:05 ------- converged at time : 0.164496 23-Apr-2021 19:44:05 number of iterations : 27 23-Apr-2021 19:44:05 number of reformations : 27 23-Apr-2021 19:44:05 ------- converged at time : 0.221682 23-Apr-2021 19:44:05 number of iterations : 30 23-Apr-2021 19:44:05 number of reformations : 30 23-Apr-2021 19:44:05 ------- converged at time : 0.276724 23-Apr-2021 19:44:05 number of iterations : 8 23-Apr-2021 19:44:05 number of reformations : 8 23-Apr-2021 19:44:05 ------- converged at time : 0.310245 23-Apr-2021 19:44:05 number of iterations : 14 23-Apr-2021 19:44:05 number of reformations : 14 23-Apr-2021 19:44:05 ------- converged at time : 0.369062 23-Apr-2021 19:44:05 number of iterations : 16 23-Apr-2021 19:44:05 number of reformations : 16 23-Apr-2021 19:44:05 ------- converged at time : 0.448115 23-Apr-2021 19:44:05 number of iterations : 17 23-Apr-2021 19:44:05 number of reformations : 17 23-Apr-2021 19:44:05 ------- converged at time : 0.543357 23-Apr-2021 19:44:05 number of iterations : 18 23-Apr-2021 19:44:05 number of reformations : 18 23-Apr-2021 19:44:05 ------- converged at time : 0.651552 23-Apr-2021 19:44:05 number of iterations : 19 23-Apr-2021 19:44:05 number of reformations : 19 23-Apr-2021 19:44:05 ------- converged at time : 0.768994 23-Apr-2021 19:44:05 number of iterations : 19 23-Apr-2021 19:44:05 number of reformations : 19 23-Apr-2021 19:44:05 ------- converged at time : 0.892695 23-Apr-2021 19:44:05 number of iterations : 16 23-Apr-2021 19:44:05 number of reformations : 16 23-Apr-2021 19:44:05 ------- converged at time : 1 23-Apr-2021 19:44:05 number of iterations : 5 23-Apr-2021 19:44:05 number of reformations : 5 23-Apr-2021 19:44:05 ------- converged at time : 1.04 23-Apr-2021 19:44:05 number of iterations : 5 23-Apr-2021 19:44:05 number of reformations : 5 23-Apr-2021 19:44:05 ------- converged at time : 1.104 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 1.1872 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 1.28576 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 1.39661 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 1.51729 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 1.64583 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 1.78066 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 1.92053 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 2 23-Apr-2021 19:44:05 number of iterations : 5 23-Apr-2021 19:44:05 number of reformations : 5 23-Apr-2021 19:44:05 ------- converged at time : 2.04 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 2.104 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 2.1872 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 2.28576 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 2.39661 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 2.51729 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 2.64583 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 2.78066 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 2.92053 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 3 23-Apr-2021 19:44:05 number of iterations : 5 23-Apr-2021 19:44:05 number of reformations : 5 23-Apr-2021 19:44:05 ------- converged at time : 3.04 23-Apr-2021 19:44:05 number of iterations : 5 23-Apr-2021 19:44:05 number of reformations : 5 23-Apr-2021 19:44:05 ------- converged at time : 3.104 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 3.1872 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 3.28576 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 3.39661 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 3.51729 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 3.64583 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 3.78066 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 3.92053 23-Apr-2021 19:44:05 number of iterations : 4 23-Apr-2021 19:44:05 number of reformations : 4 23-Apr-2021 19:44:05 ------- converged at time : 4 23-Apr-2021 19:44:05 Elapsed time : 0:00:37 23-Apr-2021 19:44:05 N O R M A L T E R M I N A T I O N # Done 23-Apr-2021 19:44:05 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Import FEBio results
if runFlag==1 %i.e. a succesful run
Importing nodal displacements from a log file
dataStruct=importFEBio_logfile(fullfile(savePath,febioLogFileName_disp),1,1); %Access data N_disp_mat=dataStruct.data; %Displacement timeVec=dataStruct.time; %Time %Create deformed coordinate set V_DEF=N_disp_mat+repmat(V,[1 1 size(N_disp_mat,3)]);
Importing nodal displacements from a log file
dataStructForce=importFEBio_logfile(fullfile(savePath,febioLogFileName_force),1,1); timeData=dataStructForce.time(:); forceData=dataStructForce.data(:); logicPush=timeData>=1 & timeData<=2; timePush=timeData(logicPush); timePush=timePush-1; forcePush=forceData(logicPush); dispPush=timePush.*prescribedDisplacement_X; logicPull=timeData>=3 & timeData<=4; timePull=timeData(logicPull); timePull=timePull-3; forcePull=forceData(logicPull); dispPull=timePull.*-prescribedDisplacement_X; Ux=[flipud(dispPull(2:end)); dispPush(2:end)]; Fx=[flipud(forcePull(2:end)); forcePush(2:end)]; dF=diff(Fx); du=diff(Ux); S_diff=dF./du; u_diff=Ux(1:end-1)+du/2; ui=linspace(-prescribedDisplacement_X,prescribedDisplacement_X,100); S=interp1(u_diff,S_diff,ui,'linear','extrap');
cFigure; hold on; xlabel('U_x [mm]'); ylabel('F_x [N]'); hp1=plot(dispPull,forcePull,'b.-','LineWidth',3,'MarkerSize',25); hp2=plot(dispPush,forcePush,'r.-','LineWidth',3,'MarkerSize',25); legend([hp1 hp2],{'Pull','Push'},'Location','NorthEastOutside') set(gca,'FontSize',fontSize); axis tight; axis square; grid on; box on; drawnow

cFigure; hold on; xlabel('U_x [mm]'); ylabel('S_x [N/mm]'); hp1=plot(u_diff,S_diff,'k.','MarkerSize',50); hp2=plot(ui,S,'k-','LineWidth',3); legend([hp1 hp2],{'FEA','Interpolated'},'Location','NorthEastOutside') set(gca,'FontSize',fontSize); axis tight; axis square; grid on; box on; drawnow

Plotting the simulated results using anim8 to visualize and animate deformations
DN_magnitude=sqrt(sum(N_disp_mat(:,:,end).^2,2)); %Current displacement magnitude % Create basic view and store graphics handle to initiate animation hf=cFigure; %Open figure gtitle([febioFebFileNamePart,': Press play to animate']); title('Displacement magnitude [mm]','Interpreter','Latex') hp=gpatch(Fb,V_DEF(:,:,end),DN_magnitude,'k',1); %Add graphics object to animate % hp.Marker='.'; % hp.MarkerSize=markerSize2; hp.FaceColor='interp'; axisGeom(gca,fontSize); colormap(gjet(250)); colorbar; caxis([0 max(DN_magnitude)]); axis(axisLim(V_DEF)); %Set axis limits statically camlight headlight; % Set up animation features animStruct.Time=timeVec; %The time vector for qt=1:1:size(N_disp_mat,3) %Loop over time increments DN_magnitude=sqrt(sum(N_disp_mat(:,:,qt).^2,2)); %Current displacement magnitude %Set entries in animation structure animStruct.Handles{qt}=[hp hp]; %Handles of objects to animate animStruct.Props{qt}={'Vertices','CData'}; %Properties of objects to animate animStruct.Set{qt}={V_DEF(:,:,qt),DN_magnitude}; %Property values for to set in order to animate end anim8(hf,animStruct); %Initiate animation feature drawnow;

Importing element stress from a log file
dataStruct=importFEBio_logfile(fullfile(savePath,febioLogFileName_stress),1,1);
%Access data
E_stress_mat=dataStruct.data;
E_stress_mat(isnan(E_stress_mat))=0;
Plotting the simulated results using anim8 to visualize and animate deformations
[CV]=faceToVertexMeasure(E,V,E_stress_mat(:,:,end)); % Create basic view and store graphics handle to initiate animation hf=cFigure; %Open figure gtitle([febioFebFileNamePart,': Press play to animate']); title('$\sigma_{1}$ [MPa]','Interpreter','Latex') hp=gpatch(Fb,V_DEF(:,:,end),CV,'k',1); %Add graphics object to animate % hp.Marker='.'; % hp.MarkerSize=markerSize2; hp.FaceColor='interp'; axisGeom(gca,fontSize); colormap(gjet(250)); colorbar; caxis([min(E_stress_mat(:)) max(E_stress_mat(:))]); axis(axisLim(V_DEF)); %Set axis limits statically camlight headlight; % Set up animation features animStruct.Time=timeVec; %The time vector for qt=1:1:size(N_disp_mat,3) %Loop over time increments [CV]=faceToVertexMeasure(E,V,E_stress_mat(:,:,qt)); %Set entries in animation structure animStruct.Handles{qt}=[hp hp]; %Handles of objects to animate animStruct.Props{qt}={'Vertices','CData'}; %Properties of objects to animate animStruct.Set{qt}={V_DEF(:,:,qt),CV}; %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) 2006-2020 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/.