PORE OR ACOUSTIC PRESSURE output from Abaqus to Matlab (Record key 108)

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\108.inp'],[S(1:a(end)-1),'\108.inp'],'f')

Run the input file 108.inp with Abaqus

!abaqus job=108

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('108.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('108.fil');

Obtain the desired output data

out = Rec108(Rec)
out =

   1.000000000000000                   0
   3.000000000000000                   0
   7.000000000000000                   0
   9.000000000000000                   0
   1.000000000000000                   0
   3.000000000000000                   0
   7.000000000000000                   0
   9.000000000000000                   0
   1.000000000000000                   0
   3.000000000000000                   0
   7.000000000000000                   0
   9.000000000000000                   0
   1.000000000000000                   0
   3.000000000000000                   0
   7.000000000000000                   0
   9.000000000000000                   0
   1.000000000000000                   0
   3.000000000000000                   0
   7.000000000000000                   0
   9.000000000000000                   0
   1.000000000000000                   0
   3.000000000000000                   0
   7.000000000000000                   0
   9.000000000000000                   0
   1.000000000000000                   0
   3.000000000000000                   0
   7.000000000000000                   0
   9.000000000000000                   0
   1.000000000000000  13.363636363636360
   3.000000000000000  13.363636363636360
   7.000000000000000                   0
   9.000000000000000                   0
   1.000000000000000                   0
   3.000000000000000  12.600000000000000
   7.000000000000000                   0
   9.000000000000000  12.600000000000000
   1.000000000000000                   0
   3.000000000000000                   0
   7.000000000000000  13.363636363636360
   9.000000000000000  13.363636363636360
   1.000000000000000  12.590184459513329
   3.000000000000000                   0
   7.000000000000000  12.602058461213939
   9.000000000000000                   0
   1.000000000000000 -14.999999999999989
   3.000000000000000 -14.999999999999989
   7.000000000000000                   0
   9.000000000000000                   0
   1.000000000000000                   0
   3.000000000000000  -9.000000000000004
   7.000000000000000                   0
   9.000000000000000  -8.999999999999998
   1.000000000000000                   0
   3.000000000000000                   0
   7.000000000000000 -15.000000000000000
   9.000000000000000 -14.999999999999989
   1.000000000000000  -9.000000000000000
   3.000000000000000                   0
   7.000000000000000  -9.000000000000004
   9.000000000000000                   0

Verify output

Check number of attributes

nAttr=size(out,2)
nAttr =

     2

Check the number of entries

nEntr=size(out,1)
nEntr =

    60

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/