DEMO_febio_0071_pneunet_actuator_set_01

Below is a demonstration for:

Contents

Keywords

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 force

%Load
appliedPressure1=0.05;
appliedPressure2=appliedPressure1/10;

%Define applied force
appliedForce=[1e-12 0 0];

%Material parameter set
c1=1; %Shear-modulus-like parameter
m1=2; %Material parameter setting degree of non-linearity
k_factor=50; %Bulk modulus factor
k=c1*k_factor; %Bulk modulus

c2=50*c1; %Shear-modulus-like parameter
m2=2; %Material parameter setting degree of non-linearity
k2=c2*k_factor; %Bulk modulus

% FEA control settings
numTimeSteps=10; %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); %Maximum time step size

runMode='internal';%'internal';
actuatorLengthDesired=100;

pointSpacing=2;
periodSize=10;
numPeriods=ceil(actuatorLengthDesired/periodSize);

heightFactor=1.4; %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;

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)=plotV(V(bcPrescribeList,:),'g.','MarkerSize',markerSize);

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';

%Control section
febio_spec.Control.analysis='STATIC';
febio_spec.Control.time_steps=numTimeSteps;
febio_spec.Control.step_size=1/numTimeSteps;
febio_spec.Control.solver.max_refs=max_refs;
febio_spec.Control.solver.max_ups=max_ups;
febio_spec.Control.time_stepper.dtmin=dtmin;
febio_spec.Control.time_stepper.dtmax=dtmax;
febio_spec.Control.time_stepper.max_retries=max_retries;
febio_spec.Control.time_stepper.opt_iter=opt_iter;

%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;

%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

% -> 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;

%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;

%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];

%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);

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:42:43
FEBio path: /home/kevin/FEBioStudio/bin/febio3
# Attempt removal of existing log files                23-Apr-2021 19:42:43
 * Removal succesful                                   23-Apr-2021 19:42:43
# Attempt removal of existing .xplt files              23-Apr-2021 19:42:43
 * Removal succesful                                   23-Apr-2021 19:42:43
# Starting FEBio...                                    23-Apr-2021 19:42:43
  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

 ]0;(0%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 1 : 0.1 =====

Setting parameter "pressure" to : 0.1
Setting parameter "pressure" to : 0.1

Reforming stiffness matrix: reformation #1

===== reforming stiffness matrix:
	Nr of equations ........................... : 11868
	Nr of nonzeroes in stiffness matrix ....... : 329070
 1
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 0.063443
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    3.551850e+00    0.000000e+00 
	   energy              1.824882e-01    1.013225e-02    1.824882e-03 
	   displacement        1.221090e+04    4.914907e+01    4.914907e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.598867
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    9.129209e+01    0.000000e+00 
	   energy              1.824882e-01    1.216891e-02    1.824882e-03 
	   displacement        1.221090e+04    3.415322e+02    6.489781e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    4.345802e-01    0.000000e+00 
	   energy              1.824882e-01    2.039877e-03    1.824882e-03 
	   displacement        1.221090e+04    2.798812e+01    9.459962e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.286993
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    3.264270e+00    0.000000e+00 
	   energy              1.824882e-01    1.059595e-03    1.824882e-03 
	   displacement        1.221090e+04    3.552995e+01    1.346656e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 0.1
	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.261073e-01    1.524637e+01    0.000000e+00 
	   energy              1.824882e-01    1.076180e-03    1.824882e-03 
	   displacement        1.221090e+04    1.097710e+02    2.223956e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 0.1
	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.261073e-01    2.480959e-01    0.000000e+00 
	   energy              1.824882e-01    6.133403e-04    1.824882e-03 
	   displacement        1.221090e+04    1.244576e+01    2.568206e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 12
	stiffness matrix reformations = 7
	step from line search         = 0.583239
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    1.348358e+00    0.000000e+00 
	   energy              1.824882e-01    1.498448e-04    1.824882e-03 
	   displacement        1.221090e+04    2.273313e+01    3.072643e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 13
	stiffness matrix reformations = 8
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    2.498408e-01    0.000000e+00 
	   energy              1.824882e-01    1.431020e-04    1.824882e-03 
	   displacement        1.221090e+04    1.195274e+01    3.466854e-03 
Reforming stiffness matrix: reformation #9

 9
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 14
	stiffness matrix reformations = 9
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    6.299071e-02    0.000000e+00 
	   energy              1.824882e-01    2.878421e-05    1.824882e-03 
	   displacement        1.221090e+04    5.799786e+00    3.755433e-03 
Reforming stiffness matrix: reformation #10

 10
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 15
	stiffness matrix reformations = 10
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    1.045630e-03    0.000000e+00 
	   energy              1.824882e-01    2.558025e-06    1.824882e-03 
	   displacement        1.221090e+04    7.447602e-01    3.861660e-03 
Reforming stiffness matrix: reformation #11

 11
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 16
	stiffness matrix reformations = 11
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    2.167983e-06    0.000000e+00 
	   energy              1.824882e-01    1.558825e-08    1.824882e-03 
	   displacement        1.221090e+04    3.370511e-02    3.884448e-03 
Reforming stiffness matrix: reformation #12

 12
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 17
	stiffness matrix reformations = 12
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.261073e-01    1.328946e-12    0.000000e+00 
	   energy              1.824882e-01    1.600941e-12    1.824882e-03 
	   displacement        1.221090e+04    2.454970e-05    3.885064e-03 

convergence summary
    number of iterations   : 12
    number of reformations : 12

------- converged at time : 0.1


Data Record #1
===========================================================================
Step = 1
Time = 0.1
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 1
Time = 0.1
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(10%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 2 : 0.2 =====

Setting parameter "pressure" to : 0.2
Setting parameter "pressure" to : 0.2

Reforming stiffness matrix: reformation #1

 1
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 0.062214
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    3.468647e+00    0.000000e+00 
	   energy              1.843426e-01    1.006262e-02    1.843426e-03 
	   displacement        1.252032e+04    4.846154e+01    4.846154e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.587952
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    9.122631e+01    0.000000e+00 
	   energy              1.843426e-01    1.268107e-02    1.843426e-03 
	   displacement        1.252032e+04    3.396452e+02    6.439014e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    4.544166e-01    0.000000e+00 
	   energy              1.843426e-01    2.074478e-03    1.843426e-03 
	   displacement        1.252032e+04    2.892757e+01    9.452040e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.287391
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    3.324551e+00    0.000000e+00 
	   energy              1.843426e-01    1.070376e-03    1.843426e-03 
	   displacement        1.252032e+04    3.577440e+01    1.347204e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 0.2
	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.271682e-01    1.541493e+01    0.000000e+00 
	   energy              1.843426e-01    1.081115e-03    1.843426e-03 
	   displacement        1.252032e+04    1.105846e+02    2.228320e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 0.2
	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.271682e-01    2.540562e-01    0.000000e+00 
	   energy              1.843426e-01    6.236412e-04    1.843426e-03 
	   displacement        1.252032e+04    1.263853e+01    2.575663e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 12
	stiffness matrix reformations = 7
	step from line search         = 0.581333
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    1.372637e+00    0.000000e+00 
	   energy              1.843426e-01    1.530616e-04    1.843426e-03 
	   displacement        1.252032e+04    2.293788e+01    3.083173e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 13
	stiffness matrix reformations = 8
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    2.595034e-01    0.000000e+00 
	   energy              1.843426e-01    1.464123e-04    1.843426e-03 
	   displacement        1.252032e+04    1.220119e+01    3.482246e-03 
Reforming stiffness matrix: reformation #9

 9
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 14
	stiffness matrix reformations = 9
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    6.572454e-02    0.000000e+00 
	   energy              1.843426e-01    2.985541e-05    1.843426e-03 
	   displacement        1.252032e+04    5.931827e+00    3.774792e-03 
Reforming stiffness matrix: reformation #10

 10
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 15
	stiffness matrix reformations = 10
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    1.140883e-03    0.000000e+00 
	   energy              1.843426e-01    2.719120e-06    1.843426e-03 
	   displacement        1.252032e+04    7.786809e-01    3.883706e-03 
Reforming stiffness matrix: reformation #11

 11
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 16
	stiffness matrix reformations = 11
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    2.504569e-06    0.000000e+00 
	   energy              1.843426e-01    1.763432e-08    1.843426e-03 
	   displacement        1.252032e+04    3.626448e-02    3.907412e-03 
Reforming stiffness matrix: reformation #12

 12
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 17
	stiffness matrix reformations = 12
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.271682e-01    2.209247e-12    0.000000e+00 
	   energy              1.843426e-01    6.337451e-12    1.843426e-03 
	   displacement        1.252032e+04    2.763688e-05    3.908068e-03 

convergence summary
    number of iterations   : 12
    number of reformations : 12

------- converged at time : 0.2


Data Record #1
===========================================================================
Step = 2
Time = 0.2
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 2
Time = 0.2
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(20%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 3 : 0.3 =====

Setting parameter "pressure" to : 0.3
Setting parameter "pressure" to : 0.3

Reforming stiffness matrix: reformation #1

 1
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 0.061388
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    3.419841e+00    0.000000e+00 
	   energy              1.859974e-01    1.003498e-02    1.859974e-03 
	   displacement        1.275359e+04    4.806217e+01    4.806217e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.580856
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    9.133135e+01    0.000000e+00 
	   energy              1.859974e-01    1.304323e-02    1.859974e-03 
	   displacement        1.275359e+04    3.387716e+02    6.412571e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    4.672572e-01    0.000000e+00 
	   energy              1.859974e-01    2.097225e-03    1.859974e-03 
	   displacement        1.275359e+04    2.953236e+01    9.453907e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.287907
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    3.365789e+00    0.000000e+00 
	   energy              1.859974e-01    1.077215e-03    1.859974e-03 
	   displacement        1.275359e+04    3.597349e+01    1.348639e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 0.3
	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.282647e-01    1.546101e+01    0.000000e+00 
	   energy              1.859974e-01    1.074322e-03    1.859974e-03 
	   displacement        1.275359e+04    1.109408e+02    2.231766e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 0.3
	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.282647e-01    2.580542e-01    0.000000e+00 
	   energy              1.859974e-01    6.291189e-04    1.859974e-03 
	   displacement        1.275359e+04    1.276988e+01    2.581239e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 12
	stiffness matrix reformations = 7
	step from line search         = 0.581941
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    1.385277e+00    0.000000e+00 
	   energy              1.859974e-01    1.538012e-04    1.859974e-03 
	   displacement        1.275359e+04    2.307024e+01    3.090802e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 13
	stiffness matrix reformations = 8
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    2.604020e-01    0.000000e+00 
	   energy              1.859974e-01    1.474475e-04    1.859974e-03 
	   displacement        1.275359e+04    1.223883e+01    3.490986e-03 
Reforming stiffness matrix: reformation #9

 9
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 14
	stiffness matrix reformations = 9
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    6.627101e-02    0.000000e+00 
	   energy              1.859974e-01    2.994621e-05    1.859974e-03 
	   displacement        1.275359e+04    5.963747e+00    3.784692e-03 
Reforming stiffness matrix: reformation #10

 10
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 15
	stiffness matrix reformations = 10
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    1.132599e-03    0.000000e+00 
	   energy              1.859974e-01    2.725356e-06    1.859974e-03 
	   displacement        1.275359e+04    7.766270e-01    3.893601e-03 
Reforming stiffness matrix: reformation #11

 11
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 16
	stiffness matrix reformations = 11
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    2.473894e-06    0.000000e+00 
	   energy              1.859974e-01    1.816204e-08    1.859974e-03 
	   displacement        1.275359e+04    3.609972e-02    3.917284e-03 
Reforming stiffness matrix: reformation #12

 12
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 17
	stiffness matrix reformations = 12
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.282647e-01    3.098814e-12    0.000000e+00 
	   energy              1.859974e-01    1.417323e-11    1.859974e-03 
	   displacement        1.275359e+04    2.758151e-05    3.917940e-03 

convergence summary
    number of iterations   : 12
    number of reformations : 12

------- converged at time : 0.3


Data Record #1
===========================================================================
Step = 3
Time = 0.3
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 3
Time = 0.3
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(30%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 4 : 0.4 =====

Setting parameter "pressure" to : 0.4
Setting parameter "pressure" to : 0.4

Reforming stiffness matrix: reformation #1

 1
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 0.060956
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    3.404377e+00    0.000000e+00 
	   energy              1.874370e-01    1.004977e-02    1.874370e-03 
	   displacement        1.290437e+04    4.794789e+01    4.794789e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.577620
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    9.168375e+01    0.000000e+00 
	   energy              1.874370e-01    1.325379e-02    1.874370e-03 
	   displacement        1.290437e+04    3.390491e+02    6.412272e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    4.713490e-01    0.000000e+00 
	   energy              1.874370e-01    2.106966e-03    1.874370e-03 
	   displacement        1.290437e+04    2.975739e+01    9.466053e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.288377
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    3.383868e+00    0.000000e+00 
	   energy              1.874370e-01    1.080168e-03    1.874370e-03 
	   displacement        1.290437e+04    3.611555e+01    1.350953e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 0.4
	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.293971e-01    1.540293e+01    0.000000e+00 
	   energy              1.874370e-01    1.060021e-03    1.874370e-03 
	   displacement        1.290437e+04    1.109070e+02    2.234591e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 0.4
	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.293971e-01    2.593112e-01    0.000000e+00 
	   energy              1.874370e-01    6.293282e-04    1.874370e-03 
	   displacement        1.290437e+04    1.282062e+01    2.584993e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 12
	stiffness matrix reformations = 7
	step from line search         = 0.584709
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    1.384709e+00    0.000000e+00 
	   energy              1.874370e-01    1.521126e-04    1.874370e-03 
	   displacement        1.290437e+04    2.311868e+01    3.095465e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 13
	stiffness matrix reformations = 8
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    2.527834e-01    0.000000e+00 
	   energy              1.874370e-01    1.461183e-04    1.874370e-03 
	   displacement        1.290437e+04    1.207177e+01    3.493110e-03 
Reforming stiffness matrix: reformation #9

 9
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 14
	stiffness matrix reformations = 9
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    6.449077e-02    0.000000e+00 
	   energy              1.874370e-01    2.906551e-05    1.874370e-03 
	   displacement        1.290437e+04    5.890068e+00    3.785042e-03 
Reforming stiffness matrix: reformation #10

 10
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 15
	stiffness matrix reformations = 10
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    1.024479e-03    0.000000e+00 
	   energy              1.874370e-01    2.576030e-06    1.874370e-03 
	   displacement        1.290437e+04    7.394058e-01    3.891294e-03 
Reforming stiffness matrix: reformation #11

 11
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 16
	stiffness matrix reformations = 11
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    2.087092e-06    0.000000e+00 
	   energy              1.874370e-01    1.712389e-08    1.874370e-03 
	   displacement        1.290437e+04    3.323251e-02    3.914009e-03 
Reforming stiffness matrix: reformation #12

 12
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 17
	stiffness matrix reformations = 12
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.293971e-01    3.753573e-12    0.000000e+00 
	   energy              1.874370e-01    2.295793e-11    1.874370e-03 
	   displacement        1.290437e+04    2.463373e-05    3.914629e-03 

convergence summary
    number of iterations   : 12
    number of reformations : 12

------- converged at time : 0.4


Data Record #1
===========================================================================
Step = 4
Time = 0.4
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 4
Time = 0.4
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(40%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 5 : 0.5 =====

Setting parameter "pressure" to : 0.5
Setting parameter "pressure" to : 0.5

Reforming stiffness matrix: reformation #1

 1
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 0.060910
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    3.421733e+00    0.000000e+00 
	   energy              1.886526e-01    1.010727e-02    1.886526e-03 
	   displacement        1.296937e+04    4.811645e+01    4.811645e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.578234
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    9.231915e+01    0.000000e+00 
	   energy              1.886526e-01    1.331167e-02    1.886526e-03 
	   displacement        1.296937e+04    3.405357e+02    6.438847e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    4.660002e-01    0.000000e+00 
	   energy              1.886526e-01    2.103221e-03    1.886526e-03 
	   displacement        1.296937e+04    2.958545e+01    9.488676e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.288738
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    3.377126e+00    0.000000e+00 
	   energy              1.886526e-01    1.079268e-03    1.886526e-03 
	   displacement        1.296937e+04    3.619572e+01    1.354142e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 0.5
	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.305652e-01    1.524961e+01    0.000000e+00 
	   energy              1.886526e-01    1.040123e-03    1.886526e-03 
	   displacement        1.296937e+04    1.105125e+02    2.236915e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 0.5
	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.305652e-01    2.575122e-01    0.000000e+00 
	   energy              1.886526e-01    6.241384e-04    1.886526e-03 
	   displacement        1.296937e+04    1.278295e+01    2.586950e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 12
	stiffness matrix reformations = 7
	step from line search         = 0.589504
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    1.370310e+00    0.000000e+00 
	   energy              1.886526e-01    1.480783e-04    1.886526e-03 
	   displacement        1.296937e+04    2.307715e+01    3.097118e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 13
	stiffness matrix reformations = 8
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    2.374481e-01    0.000000e+00 
	   energy              1.886526e-01    1.424150e-04    1.886526e-03 
	   displacement        1.296937e+04    1.170987e+01    3.488671e-03 
Reforming stiffness matrix: reformation #9

 9
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 14
	stiffness matrix reformations = 9
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    6.047211e-02    0.000000e+00 
	   energy              1.886526e-01    2.729306e-05    1.886526e-03 
	   displacement        1.296937e+04    5.710158e+00    3.775835e-03 
Reforming stiffness matrix: reformation #10

 10
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 15
	stiffness matrix reformations = 10
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    8.431802e-04    0.000000e+00 
	   energy              1.886526e-01    2.293814e-06    1.886526e-03 
	   displacement        1.296937e+04    6.715493e-01    3.876936e-03 
Reforming stiffness matrix: reformation #11

 11
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 16
	stiffness matrix reformations = 11
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    1.498860e-06    0.000000e+00 
	   energy              1.886526e-01    1.493884e-08    1.886526e-03 
	   displacement        1.296937e+04    2.824906e-02    3.897838e-03 
Reforming stiffness matrix: reformation #12

 12
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 17
	stiffness matrix reformations = 12
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.305652e-01    4.106421e-12    0.000000e+00 
	   energy              1.886526e-01    3.006556e-11    1.886526e-03 
	   displacement        1.296937e+04    2.016708e-05    3.898399e-03 

convergence summary
    number of iterations   : 12
    number of reformations : 12

------- converged at time : 0.5


Data Record #1
===========================================================================
Step = 5
Time = 0.5
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 5
Time = 0.5
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(50%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 6 : 0.6 =====

Setting parameter "pressure" to : 0.6
Setting parameter "pressure" to : 0.6

Reforming stiffness matrix: reformation #1

 1
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 0.061243
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    3.471903e+00    0.000000e+00 
	   energy              1.896434e-01    1.020756e-02    1.896434e-03 
	   displacement        1.294856e+04    4.856624e+01    4.856624e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.582640
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    9.323257e+01    0.000000e+00 
	   energy              1.896434e-01    1.321572e-02    1.896434e-03 
	   displacement        1.294856e+04    3.432115e+02    6.491969e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    4.516956e-01    0.000000e+00 
	   energy              1.896434e-01    2.086201e-03    1.896434e-03 
	   displacement        1.294856e+04    2.902894e+01    9.521718e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.289025
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    3.346518e+00    0.000000e+00 
	   energy              1.896434e-01    1.074512e-03    1.896434e-03 
	   displacement        1.294856e+04    3.621619e+01    1.358213e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 0.6
	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.317690e-01    1.499990e+01    0.000000e+00 
	   energy              1.896434e-01    1.014107e-03    1.896434e-03 
	   displacement        1.294856e+04    1.097458e+02    2.238684e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 0.6
	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.317690e-01    2.528526e-01    0.000000e+00 
	   energy              1.896434e-01    6.137413e-04    1.896434e-03 
	   displacement        1.294856e+04    1.266118e+01    2.587102e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 12
	stiffness matrix reformations = 7
	step from line search         = 0.596418
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    1.342472e+00    0.000000e+00 
	   energy              1.896434e-01    1.418074e-04    1.896434e-03 
	   displacement        1.294856e+04    2.294548e+01    3.095750e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 13
	stiffness matrix reformations = 8
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    2.156347e-01    0.000000e+00 
	   energy              1.896434e-01    1.364179e-04    1.896434e-03 
	   displacement        1.294856e+04    1.116577e+01    3.477736e-03 
Reforming stiffness matrix: reformation #9

 9
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 14
	stiffness matrix reformations = 9
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    5.453566e-02    0.000000e+00 
	   energy              1.896434e-01    2.476916e-05    1.896434e-03 
	   displacement        1.294856e+04    5.428791e+00    3.757158e-03 
Reforming stiffness matrix: reformation #10

 10
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 15
	stiffness matrix reformations = 10
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    6.286314e-04    0.000000e+00 
	   energy              1.896434e-01    1.920180e-06    1.896434e-03 
	   displacement        1.294856e+04    5.805656e-01    3.850884e-03 
Reforming stiffness matrix: reformation #11

 11
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 16
	stiffness matrix reformations = 11
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    9.118131e-07    0.000000e+00 
	   energy              1.896434e-01    1.224688e-08    1.896434e-03 
	   displacement        1.294856e+04    2.212419e-02    3.869318e-03 
Reforming stiffness matrix: reformation #12

 12
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 17
	stiffness matrix reformations = 12
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.317690e-01    4.255151e-12    0.000000e+00 
	   energy              1.896434e-01    3.398803e-11    1.896434e-03 
	   displacement        1.294856e+04    1.575378e-05    3.869812e-03 

convergence summary
    number of iterations   : 12
    number of reformations : 12

------- converged at time : 0.6


Data Record #1
===========================================================================
Step = 6
Time = 0.6
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 6
Time = 0.6
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(60%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 7 : 0.7 =====

Setting parameter "pressure" to : 0.7
Setting parameter "pressure" to : 0.7

Reforming stiffness matrix: reformation #1

 1
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 0.061949
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    3.555339e+00    0.000000e+00 
	   energy              1.904160e-01    1.035048e-02    1.904160e-03 
	   displacement        1.284520e+04    4.929592e+01    4.929592e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.590733
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    9.437842e+01    0.000000e+00 
	   energy              1.904160e-01    1.296504e-02    1.904160e-03 
	   displacement        1.284520e+04    3.469789e+02    6.570236e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    4.299875e-01    0.000000e+00 
	   energy              1.904160e-01    2.056794e-03    1.904160e-03 
	   displacement        1.284520e+04    2.812816e+01    9.564838e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.289367
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    3.295421e+00    0.000000e+00 
	   energy              1.904160e-01    1.065873e-03    1.904160e-03 
	   displacement        1.284520e+04    3.618560e+01    1.363175e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 0.7
	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.330082e-01    1.464434e+01    0.000000e+00 
	   energy              1.904160e-01    9.793800e-04    1.904160e-03 
	   displacement        1.284520e+04    1.085600e+02    2.239672e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 0.7
	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.330082e-01    2.459755e-01    0.000000e+00 
	   energy              1.904160e-01    5.986218e-04    1.904160e-03 
	   displacement        1.284520e+04    1.247063e+01    2.585413e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 12
	stiffness matrix reformations = 7
	step from line search         = 0.605752
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    1.302464e+00    0.000000e+00 
	   energy              1.904160e-01    1.334416e-04    1.904160e-03 
	   displacement        1.284520e+04    2.272874e+01    3.091372e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 13
	stiffness matrix reformations = 8
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    1.889465e-01    0.000000e+00 
	   energy              1.904160e-01    1.282961e-04    1.904160e-03 
	   displacement        1.284520e+04    1.045555e+01    3.460394e-03 
Reforming stiffness matrix: reformation #9

 9
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 14
	stiffness matrix reformations = 9
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    4.719137e-02    0.000000e+00 
	   energy              1.904160e-01    2.168368e-05    1.904160e-03 
	   displacement        1.284520e+04    5.055835e+00    3.729195e-03 
Reforming stiffness matrix: reformation #10

 10
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 15
	stiffness matrix reformations = 10
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    4.213471e-04    0.000000e+00 
	   energy              1.904160e-01    1.506946e-06    1.904160e-03 
	   displacement        1.284520e+04    4.759825e-01    3.813694e-03 
Reforming stiffness matrix: reformation #11

 11
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 16
	stiffness matrix reformations = 11
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    4.675712e-07    0.000000e+00 
	   energy              1.904160e-01    9.604478e-09    1.904160e-03 
	   displacement        1.284520e+04    1.592939e-02    3.829259e-03 
Reforming stiffness matrix: reformation #12

 12
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 17
	stiffness matrix reformations = 12
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.330082e-01    4.380441e-12    0.000000e+00 
	   energy              1.904160e-01    3.485714e-11    1.904160e-03 
	   displacement        1.284520e+04    1.236520e-05    3.829694e-03 

convergence summary
    number of iterations   : 12
    number of reformations : 12

------- converged at time : 0.7


Data Record #1
===========================================================================
Step = 7
Time = 0.7
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 7
Time = 0.7
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(70%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 8 : 0.8 =====

Setting parameter "pressure" to : 0.8
Setting parameter "pressure" to : 0.8

Reforming stiffness matrix: reformation #1

 1
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 0.063021
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    3.672880e+00    0.000000e+00 
	   energy              1.909848e-01    1.053557e-02    1.909848e-03 
	   displacement        1.266561e+04    5.030396e+01    5.030396e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.602352
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    9.567103e+01    0.000000e+00 
	   energy              1.909848e-01    1.256010e-02    1.909848e-03 
	   displacement        1.266561e+04    3.516627e+02    6.671153e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    4.031914e-01    0.000000e+00 
	   energy              1.909848e-01    2.016538e-03    1.909848e-03 
	   displacement        1.266561e+04    2.694650e+01    9.617328e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.289968
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    3.229089e+00    0.000000e+00 
	   energy              1.909848e-01    1.053313e-03    1.909848e-03 
	   displacement        1.266561e+04    3.611778e+01    1.369024e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 0.8
	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.342825e-01    1.416904e+01    0.000000e+00 
	   energy              1.909848e-01    9.320891e-04    1.909848e-03 
	   displacement        1.266561e+04    1.068845e+02    2.239514e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 0.8
	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.342825e-01    2.378200e-01    0.000000e+00 
	   energy              1.909848e-01    5.794935e-04    1.909848e-03 
	   displacement        1.266561e+04    1.223478e+01    2.581810e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 12
	stiffness matrix reformations = 7
	step from line search         = 0.617975
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    1.252129e+00    0.000000e+00 
	   energy              1.909848e-01    1.231717e-04    1.909848e-03 
	   displacement        1.266561e+04    2.243568e+01    3.083999e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 13
	stiffness matrix reformations = 8
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    1.592822e-01    0.000000e+00 
	   energy              1.909848e-01    1.183010e-04    1.909848e-03 
	   displacement        1.266561e+04    9.600465e+00    3.436764e-03 
Reforming stiffness matrix: reformation #9

 9
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 14
	stiffness matrix reformations = 9
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    3.906250e-02    0.000000e+00 
	   energy              1.909848e-01    1.826278e-05    1.909848e-03 
	   displacement        1.266561e+04    4.605395e+00    3.692221e-03 
Reforming stiffness matrix: reformation #10

 10
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 15
	stiffness matrix reformations = 10
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    2.512111e-04    0.000000e+00 
	   energy              1.909848e-01    1.105686e-06    1.909848e-03 
	   displacement        1.266561e+04    3.681650e-01    3.766116e-03 
Reforming stiffness matrix: reformation #11

 11
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 16
	stiffness matrix reformations = 11
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    2.018370e-07    0.000000e+00 
	   energy              1.909848e-01    7.309189e-09    1.909848e-03 
	   displacement        1.266561e+04    1.053586e-02    3.778695e-03 
Reforming stiffness matrix: reformation #12

 12
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 17
	stiffness matrix reformations = 12
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.342825e-01    4.648704e-12    0.000000e+00 
	   energy              1.909848e-01    3.362981e-11    1.909848e-03 
	   displacement        1.266561e+04    1.019181e-05    3.779087e-03 

convergence summary
    number of iterations   : 12
    number of reformations : 12

------- converged at time : 0.8


Data Record #1
===========================================================================
Step = 8
Time = 0.8
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 8
Time = 0.8
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(80%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 9 : 0.9 =====

Setting parameter "pressure" to : 0.9
Setting parameter "pressure" to : 0.9

Reforming stiffness matrix: reformation #1

 1
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 0.064452
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    3.825663e+00    0.000000e+00 
	   energy              1.913705e-01    1.076201e-02    1.913705e-03 
	   displacement        1.241865e+04    5.158798e+01    5.158798e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.617268
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    9.698663e+01    0.000000e+00 
	   energy              1.913705e-01    1.200448e-02    1.913705e-03 
	   displacement        1.241865e+04    3.570111e+02    6.791123e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    3.739750e-01    0.000000e+00 
	   energy              1.913705e-01    1.967556e-03    1.913705e-03 
	   displacement        1.241865e+04    2.556329e+01    9.678004e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.291079
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    3.153892e+00    0.000000e+00 
	   energy              1.913705e-01    1.036816e-03    1.913705e-03 
	   displacement        1.241865e+04    3.602993e+01    1.375720e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 0.9
	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.355913e-01    1.356109e+01    0.000000e+00 
	   energy              1.913705e-01    8.682425e-04    1.913705e-03 
	   displacement        1.241865e+04    1.046419e+02    2.237739e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 0.9
	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.355913e-01    2.294310e-01    0.000000e+00 
	   energy              1.913705e-01    5.572081e-04    1.913705e-03 
	   displacement        1.241865e+04    1.198146e+01    2.576181e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 12
	stiffness matrix reformations = 7
	step from line search         = 0.633665
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    1.193506e+00    0.000000e+00 
	   energy              1.913705e-01    1.112560e-04    1.913705e-03 
	   displacement        1.241865e+04    2.207676e+01    3.073618e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 13
	stiffness matrix reformations = 8
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    1.287134e-01    0.000000e+00 
	   energy              1.913705e-01    1.067600e-04    1.913705e-03 
	   displacement        1.241865e+04    8.628362e+00    3.407008e-03 
Reforming stiffness matrix: reformation #9

 9
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 14
	stiffness matrix reformations = 9
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    3.079935e-02    0.000000e+00 
	   energy              1.913705e-01    1.475012e-05    1.913705e-03 
	   displacement        1.241865e+04    4.094848e+00    3.646606e-03 
Reforming stiffness matrix: reformation #10

 10
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 15
	stiffness matrix reformations = 10
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    1.314331e-04    0.000000e+00 
	   energy              1.913705e-01    7.574546e-07    1.913705e-03 
	   displacement        1.241865e+04    2.668953e-01    3.709087e-03 
Reforming stiffness matrix: reformation #11

 11
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 16
	stiffness matrix reformations = 11
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    7.408204e-08    0.000000e+00 
	   energy              1.913705e-01    5.421012e-09    1.913705e-03 
	   displacement        1.241865e+04    6.423546e-03    3.718834e-03 
Reforming stiffness matrix: reformation #12

 12
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 17
	stiffness matrix reformations = 12
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.355913e-01    5.110382e-12    0.000000e+00 
	   energy              1.913705e-01    3.101882e-11    1.913705e-03 
	   displacement        1.241865e+04    8.904010e-06    3.719198e-03 

convergence summary
    number of iterations   : 12
    number of reformations : 12

------- converged at time : 0.9


Data Record #1
===========================================================================
Step = 9
Time = 0.9
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 9
Time = 0.9
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(90%) tempModel.feb - FEBio 3.2.0  
===== beginning time step 10 : 1 =====

Setting parameter "pressure" to : 1
Setting parameter "pressure" to : 1

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.066231
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.369342e-01    4.014996e+00    0.000000e+00 
	   energy              1.915995e-01    1.102850e-02    1.915995e-03 
	   displacement        1.211515e+04    5.314403e+01    5.314403e-05 
Reforming stiffness matrix: reformation #2

 2
 Nonlinear solution status: time= 1
	stiffness updates             = 0
	right hand side evaluations   = 5
	stiffness matrix reformations = 2
	step from line search         = 0.635176
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.369342e-01    9.816863e+01    0.000000e+00 
	   energy              1.915995e-01    1.130699e-02    1.915995e-03 
	   displacement        1.211515e+04    3.627021e+02    6.925473e-04 
Reforming stiffness matrix: reformation #3

 3
 Nonlinear solution status: time= 1
	stiffness updates             = 0
	right hand side evaluations   = 6
	stiffness matrix reformations = 3
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.369342e-01    3.449486e-01    0.000000e+00 
	   energy              1.915995e-01    1.912434e-03    1.915995e-03 
	   displacement        1.211515e+04    2.406537e+01    9.745118e-04 
Reforming stiffness matrix: reformation #4

 4
 Nonlinear solution status: time= 1
	stiffness updates             = 0
	right hand side evaluations   = 8
	stiffness matrix reformations = 4
	step from line search         = 0.292972
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.369342e-01    3.076492e+00    0.000000e+00 
	   energy              1.915995e-01    1.016416e-03    1.915995e-03 
	   displacement        1.211515e+04    3.594069e+01    1.383165e-03 
Reforming stiffness matrix: reformation #5

 5
 Nonlinear solution status: time= 1
	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.369342e-01    1.281433e+01    0.000000e+00 
	   energy              1.915995e-01    7.849624e-04    1.915995e-03 
	   displacement        1.211515e+04    1.017680e+02    2.233831e-03 
Reforming stiffness matrix: reformation #6

 6
 Nonlinear solution status: time= 1
	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.369342e-01    2.217891e-01    0.000000e+00 
	   energy              1.915995e-01    5.326498e-04    1.915995e-03 
	   displacement        1.211515e+04    1.173881e+01    2.568384e-03 
Reforming stiffness matrix: reformation #7

 7
 Nonlinear solution status: time= 1
	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.369342e-01    4.433787e+00    0.000000e+00 
	   energy              1.915995e-01    6.857805e-04    1.915995e-03 
	   displacement        1.211515e+04    5.073116e+01    3.338586e-03 
Reforming stiffness matrix: reformation #8

 8
 Nonlinear solution status: time= 1
	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.369342e-01    1.557046e-03    0.000000e+00 
	   energy              1.915995e-01    2.931451e-05    1.915995e-03 
	   displacement        1.211515e+04    9.223717e-01    3.450122e-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         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            2.369342e-01    1.344739e-02    0.000000e+00 
	   energy              1.915995e-01    1.787195e-06    1.915995e-03 
	   displacement        1.211515e+04    2.716077e+00    3.645863e-03 
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            2.369342e-01    2.175314e-09    0.000000e+00 
	   energy              1.915995e-01    9.403052e-09    1.915995e-03 
	   displacement        1.211515e+04    8.946896e-04    3.649464e-03 

convergence summary
    number of iterations   : 10
    number of reformations : 10

------- converged at time : 1


Data Record #1
===========================================================================
Step = 10
Time = 1
Data = ux;uy;uz
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_disp_out.txt

Data Record #2
===========================================================================
Step = 10
Time = 1
Data = s1
File = /mnt/data/MATLAB/GIBBON/data/temp/tempModel_stress_out.txt
 ]0;(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 ............ : 118

	Average number of equilibrium iterations .......... : 11.8

	Total number of right hand evaluations ............ : 167

	Total number of stiffness reformations ............ : 118


 L I N E A R   S O L V E R   S T A T S

	Total calls to linear solver ........ : 118

	Avg iterations per solve ............ : 1

	Time in linear solver: 0:00:05

 ]0;(0%) tempModel.feb - FEBio 3.2.0  
 Elapsed time : 0:00:14


 N O R M A L   T E R M I N A T I O N

 * Log file found.                                     23-Apr-2021 19:42:57
# Parsing log file...                                  23-Apr-2021 19:42:57
    number of iterations   : 12                        23-Apr-2021 19:42:58
    number of reformations : 12                        23-Apr-2021 19:42:58
------- converged at time : 0.1                        23-Apr-2021 19:42:58
    number of iterations   : 12                        23-Apr-2021 19:42:58
    number of reformations : 12                        23-Apr-2021 19:42:58
------- converged at time : 0.2                        23-Apr-2021 19:42:58
    number of iterations   : 12                        23-Apr-2021 19:42:58
    number of reformations : 12                        23-Apr-2021 19:42:58
------- converged at time : 0.3                        23-Apr-2021 19:42:58
    number of iterations   : 12                        23-Apr-2021 19:42:58
    number of reformations : 12                        23-Apr-2021 19:42:58
------- converged at time : 0.4                        23-Apr-2021 19:42:58
    number of iterations   : 12                        23-Apr-2021 19:42:58
    number of reformations : 12                        23-Apr-2021 19:42:58
------- converged at time : 0.5                        23-Apr-2021 19:42:58
    number of iterations   : 12                        23-Apr-2021 19:42:58
    number of reformations : 12                        23-Apr-2021 19:42:58
------- converged at time : 0.6                        23-Apr-2021 19:42:58
    number of iterations   : 12                        23-Apr-2021 19:42:58
    number of reformations : 12                        23-Apr-2021 19:42:58
------- converged at time : 0.7                        23-Apr-2021 19:42:58
    number of iterations   : 12                        23-Apr-2021 19:42:58
    number of reformations : 12                        23-Apr-2021 19:42:58
------- converged at time : 0.8                        23-Apr-2021 19:42:58
    number of iterations   : 12                        23-Apr-2021 19:42:58
    number of reformations : 12                        23-Apr-2021 19:42:58
------- converged at time : 0.9                        23-Apr-2021 19:42:58
    number of iterations   : 10                        23-Apr-2021 19:42:58
    number of reformations : 10                        23-Apr-2021 19:42:58
------- converged at time : 1                          23-Apr-2021 19:42:58
 Elapsed time : 0:00:14                                23-Apr-2021 19:42:58
 N O R M A L   T E R M I N A T I O N
# Done                                                 23-Apr-2021 19:42:58
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

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)]);

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, gibbon.toolbox@gmail.com

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/.