NORMALIZED CONCENTRATION (MASS DIFFUSION ANALYSIS) output from Abaqus to Matlab (Record key 221)

In this example a simple Abaqus model is analysed and results are retrieved by postprocessing the results *.fil file generated by Abaqus using Matlab. For more information please see the Documentation of Abaqus2Matlab toolbox.

Contents

Run Abaqus model

Change current directory to Abaqus working directory

a = strfind(S,'\');
cd(S(1:a(end)-1))

Copy the input file to be run by Abaqus into the Abaqus working directory

copyfile([S,'\AbaqusInputFiles\221.inp'],[S(1:a(end)-1),'\221.inp'],'f')

Run the input file 221.inp with Abaqus

!abaqus job=221

Pause Matlab execution to give Abaqus enough time to create the lck file

pause(10)

If the lck file exists then halt Matlab execution

while exist('221.lck','file')==2
  pause(0.1)
end

Postprocess Abaqus results file with Matlab

Assign all lines of the fil file in an one-row string (after Abaqus analysis terminates)

Rec = Fil2str('221.fil');

Obtain the desired output data

out = Rec221(Rec)
out =

   1.0e+04 *

   0.000100000000000                   0
   0.000200000000000                   0
   0.000300000000000   4.874005305039780
   0.000400000000000   1.624668435013260
   0.000100000000000   2.785145888594160
   0.000200000000000   2.785145888594160
   0.000300000000000   0.000000000000001
   0.000400000000000   0.000000000000000
   0.000100000000000                   0
   0.000200000000000   5.570291777188320
   0.000300000000000   5.570291777188321
   0.000400000000000                   0
   0.000100000000000                   0
   0.000200000000000                   0
   0.000300000000000   2.785145888594159
   0.000400000000000   2.785145888594159
   0.000100000000000   5.570291777188320
   0.000200000000000                   0
   0.000300000000000                   0
   0.000400000000000   5.570291777188319

Verify output

Check number of attributes

nAttr=size(out,2)
nAttr =

     2

Check the number of entries

nEntr=size(out,1)
nEntr =

    20

Check class of output

cOut=class(out)
cOut =

double

____________________________________________________
Copyright (c) 2016 by George Papazafeiropoulos
Captain, Infrastructure Engineer, Hellenic Air Force
Civil Engineer, M.Sc., Ph.D. candidate, NTUA
Email: gpapazafeiropoulos@yahoo.gr
Website: http://users.ntua.gr/gpapazaf/