DEMO_FEBio_block_biaxial_loading_stiffness_analysis

Below is a demonstration for: 1) Building an FEBio model for uniaxial compression 2) Running the model 3) Importing analysis results

Contents

clear; close all; clc;

Plot settings

fontSize=20;
faceAlpha1=0.8;
faceAlpha2=1;
edgeColor=0.25*ones(1,3);
edgeWidth=1.5;
markerSize=40;
lineWidth=3;

Control parameters

% path names
defaultFolder = fileparts(fileparts(mfilename('fullpath')));
savePath=fullfile(defaultFolder,'data','temp');

modelNameEnd='tempModel';
modelName=fullfile(savePath,modelNameEnd);

%Specifying dimensions and number of elements
sampleWidth=10;
sampleThickness=10;
sampleHeight=10;
pointSpacings=2*ones(1,3);
initialArea=sampleWidth*sampleThickness;

numElementsWidth=round(sampleWidth/pointSpacings(1));
numElementsThickness=round(sampleThickness/pointSpacings(2));
numElementsHeight=round(sampleHeight/pointSpacings(3));

boxDim=[sampleWidth sampleThickness sampleHeight]; %Dimensions
boxEl=[numElementsWidth numElementsThickness numElementsHeight]; %Number of elements

stretchLoads=[nan 1.4 0.7]; %Applied stretches, use NaN to leave face free
displacementMagnitudes=(boxDim.*stretchLoads)-boxDim;

%Material parameter set
c1=1e-3;
m1=3;
k_factor=1e2;

CREATING MESHED BOX

%Create box 1
[box1]=hexMeshBox(boxDim,boxEl);
E=box1.E;
V=box1.V;
Fb=box1.Fb;
faceBoundaryMarker=box1.faceBoundaryMarker;

X=V(:,1); Y=V(:,2); Z=V(:,3);
VE=[mean(X(E),2) mean(Y(E),2) mean(Z(E),2)];

elementMaterialIndices=ones(size(E,1),1);
% Plotting boundary surfaces
hf=cFigure;
title('Model surfaces','FontSize',fontSize);
xlabel('X','FontSize',fontSize); ylabel('Y','FontSize',fontSize); zlabel('Z','FontSize',fontSize);
hold on;
patch('Faces',Fb,'Vertices',V,'FaceColor','flat','CData',faceBoundaryMarker,'FaceAlpha',faceAlpha2,'lineWidth',edgeWidth,'edgeColor',edgeColor);

colormap(gjet(6)); icolorbar;
set(gca,'FontSize',fontSize);
view(3); axis tight;  axis equal;  grid on;
drawnow;

DEFINE BC's

%Define supported node sets
logicFace=faceBoundaryMarker==1;
Fr=Fb(logicFace,:);
bcSupportList_X=unique(Fr(:));

logicFace=faceBoundaryMarker==3;
Fr=Fb(logicFace,:);
bcSupportList_Y=unique(Fr(:));

logicFace=faceBoundaryMarker==5;
Fr=Fb(logicFace,:);
bcSupportList_Z=unique(Fr(:));

%Prescribed displacement nodes
logicPrescribe=faceBoundaryMarker==6;
Fr=Fb(logicPrescribe,:);
bcPrescribeList_Z=unique(Fr(:));
bcPrescribeMagnitudes_Z=displacementMagnitudes(ones(1,numel(bcPrescribeList_Z)),3);

logicPrescribe=faceBoundaryMarker==4;
Fr=Fb(logicPrescribe,:);
bcPrescribeList_Y=unique(Fr(:));
bcPrescribeMagnitudes_Y=displacementMagnitudes(ones(1,numel(bcPrescribeList_Y)),2);

logicPrescribe=faceBoundaryMarker==2;
Fr=Fb(logicPrescribe,:);
bcPrescribeList_X=unique(Fr(:));
bcPrescribeMagnitudes_X=displacementMagnitudes(ones(1,numel(bcPrescribeList_X)),1);

Visualize BC's

hf=cFigure;
title('Complete model','FontSize',fontSize);
xlabel('X','FontSize',fontSize); ylabel('Y','FontSize',fontSize); zlabel('Z','FontSize',fontSize);
hold on;

patch('Faces',Fb,'Vertices',V,'FaceColor','flat','CData',faceBoundaryMarker,'FaceAlpha',faceAlpha2,'lineWidth',edgeWidth,'edgeColor',edgeColor);
plotV(V(bcSupportList_X,:),'r.','MarkerSize',markerSize);
plotV(V(bcSupportList_Y,:),'g.','MarkerSize',markerSize);
plotV(V(bcSupportList_Z,:),'b.','MarkerSize',markerSize);
plotV(V(bcPrescribeList_Z,:),'k.','MarkerSize',markerSize);
plotV(V(bcPrescribeList_Y,:),'k+','MarkerSize',markerSize);
plotV(V(bcPrescribeList_X,:),'k*','MarkerSize',markerSize);
set(gca,'FontSize',fontSize);

colormap(gjet(6)); icolorbar;
set(gca,'FontSize',fontSize);
view(3); axis tight;  axis equal;  grid on;
drawnow;

CONSTRUCTING FEB MODEL

FEB_struct.febio_spec.version='2.0';
FEB_struct.Module.Type='solid';

% Defining file names
FEB_struct.run_filename=[modelName,'.feb']; %FEB file name
FEB_struct.run_logname=[modelName,'.txt']; %FEBio log file name

%Geometry section
FEB_struct.Geometry.Nodes=V;
FEB_struct.Geometry.Elements={E}; %The element sets
FEB_struct.Geometry.ElementType={'hex8'}; %The element types
FEB_struct.Geometry.ElementMat={elementMaterialIndices};
FEB_struct.Geometry.ElementsPartName={'Block'};

%Material section
k=c1*k_factor;

%Material section
FEB_struct.Materials{1}.Type='Ogden';
FEB_struct.Materials{1}.Name='Block_material';
FEB_struct.Materials{1}.Properties={'c1','m1','c2','m2','k'};
FEB_struct.Materials{1}.Values={c1,m1,c1,-m1,k};

%Step specific control sections
FEB_struct.Control.AnalysisType='static';
FEB_struct.Control.Properties={'time_steps','step_size',...
    'max_refs','max_ups',...
    'dtol','etol','rtol','lstol'};
numSteps=20;
FEB_struct.Control.Values={numSteps,1/numSteps,25,5,0.001,0.01,0,0.9};
FEB_struct.Control.TimeStepperProperties={'dtmin','dtmax','max_retries','opt_iter','aggressiveness'};
FEB_struct.Control.TimeStepperValues={(1/(100*numSteps)),1/numSteps,5,10,1};

%Defining node sets
FEB_struct.Geometry.NodeSet{1}.Set=bcSupportList_X;
FEB_struct.Geometry.NodeSet{1}.Name='bcSupportList_X';
FEB_struct.Geometry.NodeSet{2}.Set=bcSupportList_Y;
FEB_struct.Geometry.NodeSet{2}.Name='bcSupportList_Y';
FEB_struct.Geometry.NodeSet{3}.Set=bcSupportList_Z;
FEB_struct.Geometry.NodeSet{3}.Name='bcSupportList_Z';


%Adding BC information
FEB_struct.Boundary.Fix{1}.bc='x';
FEB_struct.Boundary.Fix{1}.SetName=FEB_struct.Geometry.NodeSet{1}.Name;
FEB_struct.Boundary.Fix{2}.bc='y';
FEB_struct.Boundary.Fix{2}.SetName=FEB_struct.Geometry.NodeSet{2}.Name;
FEB_struct.Boundary.Fix{3}.bc='z';
FEB_struct.Boundary.Fix{3}.SetName=FEB_struct.Geometry.NodeSet{3}.Name;

%Prescribed BC's
indBC=1;
if ~any(isnan(bcPrescribeMagnitudes_X(:)))
    FEB_struct.Boundary.Prescribe{indBC}.Set=bcPrescribeList_X;
    FEB_struct.Boundary.Prescribe{indBC}.bc='x';
    FEB_struct.Boundary.Prescribe{indBC}.lc=1;
    FEB_struct.Boundary.Prescribe{indBC}.nodeScale=bcPrescribeMagnitudes_X;
    FEB_struct.Boundary.Prescribe{indBC}.Type='relative';
    indBC=indBC+1;
end

if ~any(isnan(bcPrescribeMagnitudes_Y(:)))
    FEB_struct.Boundary.Prescribe{indBC}.Set=bcPrescribeList_Y;
    FEB_struct.Boundary.Prescribe{indBC}.bc='y';
    FEB_struct.Boundary.Prescribe{indBC}.lc=1;
    FEB_struct.Boundary.Prescribe{indBC}.nodeScale=bcPrescribeMagnitudes_Y;
    FEB_struct.Boundary.Prescribe{indBC}.Type='relative';
    indBC=indBC+1;
end

if ~any(isnan(bcPrescribeMagnitudes_Z(:)))
    FEB_struct.Boundary.Prescribe{indBC}.Set=bcPrescribeList_Z;
    FEB_struct.Boundary.Prescribe{indBC}.bc='z';
    FEB_struct.Boundary.Prescribe{indBC}.lc=1;
    FEB_struct.Boundary.Prescribe{indBC}.nodeScale=bcPrescribeMagnitudes_Z;
    FEB_struct.Boundary.Prescribe{indBC}.Type='relative';
end

%Load curves
FEB_struct.LoadData.LoadCurves.id=1;
FEB_struct.LoadData.LoadCurves.type={'linear'};
FEB_struct.LoadData.LoadCurves.loadPoints={[0 0;1 1;]};

%Adding output requests
FEB_struct.Output.VarTypes={'displacement','stress','relative volume'};

%Specify log file output
run_disp_output_name=[modelNameEnd,'_node_out.txt'];
run_force_output_name=[modelNameEnd,'_force_out.txt'];
run_stiffness_output_name=[modelNameEnd,'_stiffness_out.txt'];
FEB_struct.run_output_names={run_disp_output_name,run_force_output_name,run_stiffness_output_name};
FEB_struct.output_types={'node_data','node_data','element_data'};
FEB_struct.data_types={'ux;uy;uz','Rx;Ry;Rz','cxxxx;cxxyy;cyyyy;cxxzz;cyyzz;czzzz;cxxxy;cyyxy;czzxy;cxyxy;cxxyz;cyyyz;czzyz;cxyyz;cyzyz;cxxxz;cyyxz;czzxz;cxyxz;cyzxz;cxzxz'};

SAVING .FEB FILE

FEB_struct.disp_opt=0; %Display waitbars
febStruct2febFile(FEB_struct);
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
--- Writing FEBio XML object --- 18-Aug-2017 14:36:18
Adding Module level
Adding Control level
Adding Globals level
Adding Material level
Adding Geometry level
----> Adding node field
----> Adding element field
----> Adding hex8 element entries....
----> Adding NodeSet field
Adding Boundary level
----> Defining fix type boundary conditions
----> Defining prescribe type boundary conditions
Adding LoadData level
----> Defining load curves
Adding Output level
----> Adding plotfile field
----> Adding logfile field
Writing .feb file
--- Done --- 18-Aug-2017 14:36:18

RUNNING FEBIO JOB

FEBioRunStruct.run_filename=FEB_struct.run_filename;
FEBioRunStruct.run_logname=FEB_struct.run_logname;
FEBioRunStruct.disp_on=1;
FEBioRunStruct.disp_log_on=1;
FEBioRunStruct.runMode='internal';%'internal';
FEBioRunStruct.t_check=0.25; %Time for checking log file (dont set too small)
FEBioRunStruct.maxtpi=1e99; %Max analysis time
FEBioRunStruct.maxLogCheckTime=3; %Max log file checking time

[runFlag]=runMonitorFEBio(FEBioRunStruct);%START FEBio NOW!!!!!!!!
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
--- STARTING FEBIO JOB --- 18-Aug-2017 14:36:18
===========================================================================
         ________    _________   _________     __     _________            
        |        |\ |        |\ |        |\   |  |\  /         \\          
        |    ____|| |    ____|| |    __  ||   |__|| |    ___    ||         
        |   |\___\| |   |\___\| |   |\_| ||    \_\| |   //  \   ||         
        |   ||      |   ||      |   || | ||    __   |  ||    |  ||         
        |   ||__    |   ||__    |   ||_| ||   |  |\ |  ||    |  ||         
        |       |\  |       |\  |         \\  |  || |  ||    |  ||         
        |    ___||  |    ___||  |    ___   || |  || |  ||    |  ||         
        |   |\__\|  |   |\__\|  |   |\__|  || |  || |  ||    |  ||         
        |   ||      |   ||      |   ||  |  || |  || |  ||    |  ||         
        |   ||      |   ||___   |   ||__|  || |  || |   \\__/   ||         
        |   ||      |        |\ |          || |  || |           ||         
        |___||      |________|| |__________|| |__||  \_________//          
                                                                           
      F I N I T E   E L E M E N T S   F O R   B I O M E C H A N I C S      
                                                                           
                 --- v e r s i o n - 2 . 5 . 2 . 8980 ---                 
                                                                           
                                                                           
  Musculoskeletal Research Laboratory                                      
  University of Utah                                                       
  http://febio.org                                                         
                                                                           
  FEBio is a registered trademark.                                         
  copyright (c) 2006-2015 - All rights reserved                            
                                                                           
 This is the NON-COMMERCIAL version of FEBio. This version may only be          
 used for non-commercial purposes as described in the license agreement.        
 To obtain a valid commercial license file, please contact the developers.      
                                                                           
===========================================================================


]0;(0%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 1 : 0.05 =====
===== reforming stiffness matrix:
	Nr of equations ........................... : 540
	Nr of nonzeroes in stiffness matrix ....... : 14670

 1
 Nonlinear solution status: time= 0.05
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.620219e-02    3.375603e-07    0.000000e+00 
	   energy              2.557677e-01    1.582034e-04    2.557677e-03 
	   displacement        2.890261e+00    2.890261e+00    2.890261e-06 
 2
 Nonlinear solution status: time= 0.05
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.620219e-02    1.113326e-10    0.000000e+00 
	   energy              2.557677e-01    8.743459e-08    2.557677e-03 
	   displacement        2.890261e+00    8.016046e-04    2.866578e-06 
 3
 Nonlinear solution status: time= 0.05
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.620219e-02    9.398996e-14    0.000000e+00 
	   energy              2.557677e-01    1.689544e-12    2.557677e-03 
	   displacement        2.890261e+00    6.022960e-08    2.866585e-06 


------- converged at time : 0.05

]0;(5%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 2 : 0.1 =====
 1
 Nonlinear solution status: time= 0.1
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.478305e-02    3.259777e-07    0.000000e+00 
	   energy              2.516524e-01    1.355890e-04    2.516524e-03 
	   displacement        2.844594e+00    2.844594e+00    2.844594e-06 
 2
 Nonlinear solution status: time= 0.1
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.478305e-02    9.897860e-11    0.000000e+00 
	   energy              2.516524e-01    7.328997e-08    2.516524e-03 
	   displacement        2.844594e+00    7.602255e-04    2.824590e-06 
 3
 Nonlinear solution status: time= 0.1
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.478305e-02    8.682706e-14    0.000000e+00 
	   energy              2.516524e-01    1.075551e-12    2.516524e-03 
	   displacement        2.844594e+00    4.481773e-08    2.824549e-06 


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

]0;(10%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 3 : 0.15 =====
 1
 Nonlinear solution status: time= 0.15
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.370704e-02    3.161010e-07    0.000000e+00 
	   energy              2.481273e-01    1.143785e-04    2.481273e-03 
	   displacement        2.806098e+00    2.806098e+00    2.806098e-06 
 2
 Nonlinear solution status: time= 0.15
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.370704e-02    8.913089e-11    0.000000e+00 
	   energy              2.481273e-01    6.064284e-08    2.481273e-03 
	   displacement        2.806098e+00    7.257085e-04    2.789474e-06 
 3
 Nonlinear solution status: time= 0.15
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.370704e-02    8.154085e-14    0.000000e+00 
	   energy              2.481273e-01    4.885812e-13    2.481273e-03 
	   displacement        2.806098e+00    3.271278e-08    2.789395e-06 


------- converged at time : 0.15

]0;(15%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 4 : 0.2 =====
 1
 Nonlinear solution status: time= 0.2
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.295625e-02    3.078250e-07    0.000000e+00 
	   energy              2.451691e-01    9.438794e-05    2.451691e-03 
	   displacement        2.774182e+00    2.774182e+00    2.774182e-06 
 2
 Nonlinear solution status: time= 0.2
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.295625e-02    8.145823e-11    0.000000e+00 
	   energy              2.451691e-01    4.921258e-08    2.451691e-03 
	   displacement        2.774182e+00    6.973424e-04    2.760691e-06 
 3
 Nonlinear solution status: time= 0.2
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.295625e-02    7.798462e-14    0.000000e+00 
	   energy              2.451691e-01    5.044864e-14    2.451691e-03 
	   displacement        2.774182e+00    2.330556e-08    2.760581e-06 


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

]0;(20%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 5 : 0.25 =====
 1
 Nonlinear solution status: time= 0.25
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.251733e-02    3.010627e-07    0.000000e+00 
	   energy              2.427587e-01    7.539920e-05    2.427587e-03 
	   displacement        2.748387e+00    2.748387e+00    2.748387e-06 
 2
 Nonlinear solution status: time= 0.25
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.251733e-02    7.568965e-11    0.000000e+00 
	   energy              2.427587e-01    3.875812e-08    2.427587e-03 
	   displacement        2.748387e+00    6.745582e-04    2.737828e-06 
 3
 Nonlinear solution status: time= 0.25
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.251733e-02    7.596286e-14    0.000000e+00 
	   energy              2.427587e-01    5.573621e-13    2.427587e-03 
	   displacement        2.748387e+00    1.613520e-08    2.737694e-06 


------- converged at time : 0.25

]0;(25%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 6 : 0.3 =====
 1
 Nonlinear solution status: time= 0.3
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.238128e-02    2.957427e-07    0.000000e+00 
	   energy              2.408806e-01    5.720898e-05    2.408806e-03 
	   displacement        2.728344e+00    2.728344e+00    2.728344e-06 
 2
 Nonlinear solution status: time= 0.3
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.238128e-02    7.162155e-11    0.000000e+00 
	   energy              2.408806e-01    2.906662e-08    2.408806e-03 
	   displacement        2.728344e+00    6.569061e-04    2.720560e-06 
 3
 Nonlinear solution status: time= 0.3
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.238128e-02    7.535562e-14    0.000000e+00 
	   energy              2.408806e-01    1.046312e-12    2.408806e-03 
	   displacement        2.728344e+00    1.085383e-08    2.720408e-06 


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

]0;(30%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 7 : 0.35 =====
 1
 Nonlinear solution status: time= 0.35
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.254311e-02    2.918080e-07    0.000000e+00 
	   energy              2.395230e-01    3.962567e-05    2.395230e-03 
	   displacement        2.713767e+00    2.713767e+00    2.713767e-06 
 2
 Nonlinear solution status: time= 0.35
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.254311e-02    6.910876e-11    0.000000e+00 
	   energy              2.395230e-01    1.994714e-08    2.395230e-03 
	   displacement        2.713767e+00    6.440411e-04    2.708643e-06 
 3
 Nonlinear solution status: time= 0.35
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.254311e-02    7.610931e-14    0.000000e+00 
	   energy              2.395230e-01    1.530360e-12    2.395230e-03 
	   displacement        2.713767e+00    7.206477e-09    2.708477e-06 


------- converged at time : 0.35

]0;(35%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 8 : 0.4 =====
 1
 Nonlinear solution status: time= 0.4
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.300177e-02    2.892157e-07    0.000000e+00 
	   energy              2.386775e-01    2.246669e-05    2.386775e-03 
	   displacement        2.704450e+00    2.704450e+00    2.704450e-06 
 2
 Nonlinear solution status: time= 0.4
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.300177e-02    6.805804e-11    0.000000e+00 
	   energy              2.386775e-01    1.122515e-08    2.386775e-03 
	   displacement        2.704450e+00    6.357127e-04    2.701906e-06 
 3
 Nonlinear solution status: time= 0.4
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.300177e-02    7.823201e-14    0.000000e+00 
	   energy              2.386775e-01    2.022024e-12    2.386775e-03 
	   displacement        2.704450e+00    5.016476e-09    2.701732e-06 


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

]0;(40%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 9 : 0.45 =====
 1
 Nonlinear solution status: time= 0.45
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.376003e-02    2.879357e-07    0.000000e+00 
	   energy              2.383389e-01    5.555900e-06    2.383389e-03 
	   displacement        2.700259e+00    2.700259e+00    2.700259e-06 
 2
 Nonlinear solution status: time= 0.45
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.376003e-02    6.842401e-11    0.000000e+00 
	   energy              2.383389e-01    2.737616e-09    2.383389e-03 
	   displacement        2.700259e+00    6.317565e-04    2.700250e-06 
 3
 Nonlinear solution status: time= 0.45
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.376003e-02    8.179284e-14    0.000000e+00 
	   energy              2.383389e-01    2.533779e-12    2.383389e-03 
	   displacement        2.700259e+00    4.175422e-09    2.700073e-06 


------- converged at time : 0.45

]0;(45%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 10 : 0.5 =====
 1
 Nonlinear solution status: time= 0.5
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.482451e-02    2.879508e-07    0.000000e+00 
	   energy              2.385055e-01    1.127878e-05    2.385055e-03 
	   displacement        2.701129e+00    2.701129e+00    2.701129e-06 
 2
 Nonlinear solution status: time= 0.5
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.482451e-02    7.020698e-11    0.000000e+00 
	   energy              2.385055e-01    5.671394e-09    2.385055e-03 
	   displacement        2.701129e+00    6.320895e-04    2.703647e-06 
 3
 Nonlinear solution status: time= 0.5
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.482451e-02    8.692546e-14    0.000000e+00 
	   energy              2.385055e-01    3.078537e-12    2.385055e-03 
	   displacement        2.701129e+00    4.636914e-09    2.703471e-06 


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

]0;(50%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 11 : 0.55 =====
 1
 Nonlinear solution status: time= 0.55
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.620581e-02    2.892561e-07    0.000000e+00 
	   energy              2.391785e-01    2.820769e-05    2.391785e-03 
	   displacement        2.707066e+00    2.707066e+00    2.707066e-06 
 2
 Nonlinear solution status: time= 0.55
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.620581e-02    7.345281e-11    0.000000e+00 
	   energy              2.391785e-01    1.415490e-08    2.391785e-03 
	   displacement        2.707066e+00    6.367071e-04    2.712136e-06 
 3
 Nonlinear solution status: time= 0.55
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.620581e-02    9.383572e-14    0.000000e+00 
	   energy              2.391785e-01    3.670139e-12    2.391785e-03 
	   displacement        2.707066e+00    6.413701e-09    2.711965e-06 


------- converged at time : 0.55

]0;(55%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 12 : 0.6 =====
 1
 Nonlinear solution status: time= 0.6
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.791873e-02    2.918588e-07    0.000000e+00 
	   energy              2.403627e-01    4.540168e-05    2.403627e-03 
	   displacement        2.718143e+00    2.718143e+00    2.718143e-06 
 2
 Nonlinear solution status: time= 0.6
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.791873e-02    7.825463e-11    0.000000e+00 
	   energy              2.403627e-01    2.286701e-08    2.403627e-03 
	   displacement        2.718143e+00    6.456826e-04    2.725826e-06 
 3
 Nonlinear solution status: time= 0.6
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.791873e-02    1.028142e-13    0.000000e+00 
	   energy              2.403627e-01    4.323873e-12    2.403627e-03 
	   displacement        2.718143e+00    9.577989e-09    2.725665e-06 


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

]0;(60%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 13 : 0.65 =====
 1
 Nonlinear solution status: time= 0.65
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.998251e-02    2.957781e-07    0.000000e+00 
	   energy              2.420661e-01    6.303435e-05    2.420661e-03 
	   displacement        2.734505e+00    2.734505e+00    2.734505e-06 
 2
 Nonlinear solution status: time= 0.65
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.998251e-02    8.475653e-11    0.000000e+00 
	   energy              2.420661e-01    3.196695e-08    2.420661e-03 
	   displacement        2.734505e+00    6.591690e-04    2.744895e-06 
 3
 Nonlinear solution status: time= 0.65
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            5.998251e-02    1.142546e-13    0.000000e+00 
	   energy              2.420661e-01    5.057050e-12    2.420661e-03 
	   displacement        2.734505e+00    1.426491e-08    2.744749e-06 


------- converged at time : 0.65

]0;(65%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 14 : 0.7 =====
 1
 Nonlinear solution status: time= 0.7
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            6.242133e-02    3.010453e-07    0.000000e+00 
	   energy              2.443000e-01    8.128434e-05    2.443000e-03 
	   displacement        2.756367e+00    2.756367e+00    2.756367e-06 
 2
 Nonlinear solution status: time= 0.7
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            6.242133e-02    9.315948e-11    0.000000e+00 
	   energy              2.443000e-01    4.162321e-08    2.443000e-03 
	   displacement        2.756367e+00    6.774028e-04    2.769597e-06 
 3
 Nonlinear solution status: time= 0.7
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            6.242133e-02    1.286793e-13    0.000000e+00 
	   energy              2.443000e-01    5.889666e-12    2.443000e-03 
	   displacement        2.756367e+00    2.067937e-08    2.769473e-06 


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

]0;(70%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 15 : 0.75 =====
 1
 Nonlinear solution status: time= 0.75
	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.526478e-02    3.077043e-07    0.000000e+00 
	   energy              2.470794e-01    1.003377e-04    2.470794e-03 
	   displacement        2.784024e+00    2.784024e+00    2.784024e-06 
 2
 Nonlinear solution status: time= 0.75
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            6.526478e-02    1.037297e-10    0.000000e+00 
	   energy              2.470794e-01    5.201799e-08    2.470794e-03 
	   displacement        2.784024e+00    7.007108e-04    2.800268e-06 
 3
 Nonlinear solution status: time= 0.75
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            6.526478e-02    1.467745e-13    0.000000e+00 
	   energy              2.470794e-01    6.845177e-12    2.470794e-03 
	   displacement        2.784024e+00    2.910664e-08    2.800170e-06 


------- converged at time : 0.75

]0;(75%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 16 : 0.8 =====
 1
 Nonlinear solution status: time= 0.8
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            6.854855e-02    3.158114e-07    0.000000e+00 
	   energy              2.504232e-01    1.203904e-04    2.504232e-03 
	   displacement        2.817853e+00    2.817853e+00    2.817853e-06 
 2
 Nonlinear solution status: time= 0.8
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            6.854855e-02    1.168101e-10    0.000000e+00 
	   energy              2.504232e-01    6.335201e-08    2.504232e-03 
	   displacement        2.817853e+00    7.295187e-04    2.837328e-06 
 3
 Nonlinear solution status: time= 0.8
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            6.854855e-02    1.694382e-13    0.000000e+00 
	   energy              2.504232e-01    7.951429e-12    2.504232e-03 
	   displacement        2.817853e+00    3.992742e-08    2.837264e-06 


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

]0;(80%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 17 : 0.85 =====
 1
 Nonlinear solution status: time= 0.85
	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.231531e-02    3.254358e-07    0.000000e+00 
	   energy              2.543544e-01    1.416508e-04    2.543544e-03 
	   displacement        2.858321e+00    2.858321e+00    2.858321e-06 
 2
 Nonlinear solution status: time= 0.85
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            7.231531e-02    1.328349e-10    0.000000e+00 
	   energy              2.543544e-01    7.584997e-08    2.543544e-03 
	   displacement        2.858321e+00    7.643635e-04    2.881292e-06 
 3
 Nonlinear solution status: time= 0.85
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            7.231531e-02    1.978445e-13    0.000000e+00 
	   energy              2.543544e-01    9.241794e-12    2.543544e-03 
	   displacement        2.858321e+00    5.363830e-08    2.881271e-06 


------- converged at time : 0.85

]0;(85%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 18 : 0.9 =====
 1
 Nonlinear solution status: time= 0.9
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            7.661573e-02    3.366606e-07    0.000000e+00 
	   energy              2.589003e-01    1.643427e-04    2.589003e-03 
	   displacement        2.905998e+00    2.905998e+00    2.905998e-06 
 2
 Nonlinear solution status: time= 0.9
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            7.661573e-02    1.523494e-10    0.000000e+00 
	   energy              2.589003e-01    8.976672e-08    2.589003e-03 
	   displacement        2.905998e+00    8.059088e-04    2.932783e-06 
 3
 Nonlinear solution status: time= 0.9
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            7.661573e-02    2.335312e-13    0.000000e+00 
	   energy              2.589003e-01    1.075655e-11    2.589003e-03 
	   displacement        2.905998e+00    7.087891e-08    2.932814e-06 


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

]0;(90%) tempModel.feb - FEBio 2.5.2.8980 
===== beginning time step 19 : 0.95 =====
 1
 Nonlinear solution status: time= 0.95
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            8.150979e-02    3.495826e-07    0.000000e+00 
	   energy              2.640932e-01    1.887085e-04    2.640932e-03 
	   displacement        2.961565e+00    2.961565e+00    2.961565e-06 
 2
 Nonlinear solution status: time= 0.95
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            8.150979e-02    1.760343e-10    0.000000e+00 
	   energy              2.640932e-01    1.053944e-07    2.640932e-03 
	   displacement        2.961565e+00    8.549645e-04    2.992543e-06 
 3
 Nonlinear solution status: time= 0.95
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            8.150979e-02    2.785186e-13    0.000000e+00 
	   energy              2.640932e-01    1.254460e-11    2.640932e-03 
	   displacement        2.961565e+00    9.246764e-08    2.992639e-06 


------- converged at time : 0.95

]0;(95%) tempModel.feb - FEBio 2.5.2.8980 MUST POINT CONTROLLER: adjusting time step. dt = 0.05


===== beginning time step 20 : 1 =====
 1
 Nonlinear solution status: time= 1
	stiffness updates             = 0
	right hand side evaluations   = 2
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            8.706834e-02    3.643139e-07    0.000000e+00 
	   energy              2.699706e-01    2.150125e-04    2.699706e-03 
	   displacement        3.025832e+00    3.025832e+00    3.025832e-06 
 2
 Nonlinear solution status: time= 1
	stiffness updates             = 1
	right hand side evaluations   = 3
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            8.706834e-02    2.047381e-10    0.000000e+00 
	   energy              2.699706e-01    1.230710e-07    2.699706e-03 
	   displacement        3.025832e+00    9.125107e-04    3.061452e-06 
 3
 Nonlinear solution status: time= 1
	stiffness updates             = 2
	right hand side evaluations   = 4
	stiffness matrix reformations = 1
	step from line search         = 1.000000
	convergence norms :     INITIAL         CURRENT         REQUIRED
	   residual            8.706834e-02    3.354715e-13    0.000000e+00 
	   energy              2.699706e-01    1.466562e-11    2.699706e-03 
	   displacement        3.025832e+00    1.194482e-07    3.061625e-06 


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

]0;(100%) tempModel.feb - FEBio 2.5.2.8980 

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 .................... : 20

	Total number of equilibrium iterations ............ : 60

	Average number of equilibrium iterations .......... : 3

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

	Total number of stiffness reformations ............ : 20

	Time in linear solver: 0:00:00


 Elapsed time : 0:00:02


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

Waiting for log file...
Proceeding to check log file...18-Aug-2017 14:36:20
------- converged at time : 0.05
------- converged at time : 0.1
------- converged at time : 0.15
------- converged at time : 0.2
------- converged at time : 0.25
------- converged at time : 0.3
------- converged at time : 0.35
------- converged at time : 0.4
------- converged at time : 0.45
------- converged at time : 0.5
------- converged at time : 0.55
------- converged at time : 0.6
------- converged at time : 0.65
------- converged at time : 0.7
------- converged at time : 0.75
------- converged at time : 0.8
------- converged at time : 0.85
------- converged at time : 0.9
------- converged at time : 0.95
------- converged at time : 1
--- Done --- 18-Aug-2017 14:36:20
if runFlag==1 %i.e. a succesful run

IMPORTING NODAL DISPLACEMENT RESULTS

Importing nodal displacements from a log file

    [~, N_disp_mat,~]=importFEBio_logfile(fullfile(savePath,FEB_struct.run_output_names{1})); %Nodal displacements

IMPORTING NODAL FORCES

Importing nodal forces from a log file

    [time_mat, N_force_mat,~]=importFEBio_logfile(fullfile(savePath,FEB_struct.run_output_names{2})); %Nodal forces
    time_mat=[0; time_mat(:)]; %Time

IMPORTING ELEMENT STIFFNESS MATRICES

Importing element stiffness tensors from a log file

    [~,stiffness_mat,~]=importFEBio_logfile(fullfile(savePath,FEB_struct.run_output_names{3})); %Nodal forces

    stiffness_mat=stiffness_mat(:,2:end,end); %Final stiffness state

    stiffness_mat_voigt=stiffness_mat(:,[1  2  4  11 16 7;...
                                         2  3  5  12 17 8;...
                                         4  5  6  13 18 9;...
                                         11 12 13 15 20 14;...
                                         16 17 18 20 21 19;...
                                         7  8  9  14 19 10]);
    stiffness_mat_voigt=reshape(stiffness_mat_voigt',6,6,size(stiffness_mat_voigt,1));
    stiffness_mat_voigt=reshape(mat2cell(stiffness_mat_voigt,6,6,...
        ones(size(stiffness_mat_voigt,3),1)),[size(stiffness_mat,1),1]);

    stiffness_mat_kelvin=stiffness_mat_voigt;
    for q=1:1:numel(stiffness_mat_voigt)
        cVoigt=stiffness_mat_voigt{q};
        c=voigtUnMap(cVoigt);
        cKelvin=kelvinMap(c);
        stiffness_mat_kelvin{q}=cKelvin;

    end
    viewFourthOrderTensor(c); %Visualize tensor C

Plotting the deformed model

    N_disp_mat=N_disp_mat(:,2:end,:);
    sizImport=size(N_disp_mat);
    sizImport(3)=sizImport(3)+1;
    N_disp_mat_n=zeros(sizImport);
    N_disp_mat_n(:,:,2:end)=N_disp_mat;
    N_disp_mat=N_disp_mat_n;
    DN=N_disp_mat(:,:,end);
    DN_magnitude=sqrt(sum(DN(:,3).^2,2));
    V_def=V+DN;
    [CF]=vertexToFaceMeasure(Fb,DN_magnitude);
    hf=cFigure;
    xlabel('X','FontSize',fontSize); ylabel('Y','FontSize',fontSize); zlabel('Z','FontSize',fontSize); hold on;

    hp=gpatch(Fb,V_def,CF,'k',1);
    gpatch(Fb,V,0.5*ones(1,3),'k',0.25);

    view(3); axis tight;  axis equal;  grid on; box on;
    colormap(gjet(250)); colorbar;
    caxis([0 max(DN_magnitude)]);
    axis([min(V_def(:,1)) max(V_def(:,1)) min(V_def(:,2)) max(V_def(:,2)) min(V(:,3)) max(V(:,3))]);
    view(130,25);
    camlight headlight;
    set(gca,'FontSize',fontSize);
    drawnow;

    animStruct.Time=time_mat;

    for qt=1:1:size(N_disp_mat,3)

        DN=N_disp_mat(:,:,qt);
        DN_magnitude=sqrt(sum(DN(:,3).^2,2));
        V_def=V+DN;
        [CF]=vertexToFaceMeasure(Fb,DN_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,CF}; %Property values for to set in order to animate
    end

    anim8(hf,animStruct);
    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) 2017 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/.