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/