Could you please help me on how to compute the averaged stresses (=sum of all stresses at all integration points divided by the model volume) once I already finish the Abaqus analysis?
Thanks.
P
If using this toolbox for research or industrial purposes, please cite:
Advances in Engineering Software. Vol 105. March 2017. Pages 9-16. (2017)
Hello Nguyen,
You can proceed in two ways:
1) Assuming that you are running an Abaqus/Standard analysis, specify the following options in the Abaqus input file:
*FILE FORMAT, ASCII
*EL FILE
S, EVOL
This will output the stresses (S, record key 11) and element volume (EVOL, record key 78) for all elements in your model into a results file (for example FILENAME.fil). After this you have to execute in Matlab:
out=readFil ( 'FILENAME.fil', [11 78], 'C:\Abaqus_Temp')
where 'FILENAME.fil' is the name of the results file which will be generated after the Abaqus analysis of the above input file, and 'C:\Abaqus_Temp' (or something similar, depending on your PC) is the directory in which the 'FILENAME.fil' file exists. 11 and 78 are the record keys of the stresses and element volumes respectively.
2) Specify the following options in the Abaqus input file:
*OUTPUT, FIELD
*ELEMENT OUTPUT
S, EVOL
or
*OUTPUT, HISTORY
*ELEMENT OUTPUT
S, EVOL
This will output the stresses (S) and element volume (EVOL) for all elements in your model into the odb file (for example FILENAME.odb). either as field output (first case) or history output (second case). After this you have to execute in Matlab:
odbOut1 = readElementFieldOdb ( 'FILENAME.odb', 'False', 'False', 'False', 'S', 'False', 'C:\Abaqus_Temp')
odbOut2 = readElementFieldOdb ( 'FILENAME.odb', 'False', 'False', 'False', 'EVOL', 'False', 'C:\Abaqus_Temp')
or
odbOut1 = readHistoryOdb ( 'FILENAME.odb', 'False', 'S', 'False', 'C:\Abaqus_Temp')
odbOut2 = readHistoryOdb ( 'FILENAME.odb', 'False', 'EVOL', 'False', 'C:\Abaqus_Temp')
for field and history output respectively. 'FILENAME.odb' is the name of the odb file which will be generated after the Abaqus analysis of the above input file, and 'C:\Abaqus_Temp' (or something similar, depending on your PC) is the directory in which the 'FILENAME.odb' file exists.
More details regarding the above can be found in the documentation of Abaqus2Matlab. After extracting the above output, you have to calculate the averaged stresses as follows:
1) For the first case from above:
AvgStress = sum(out{1})/sum(out{2})
2) For the second case from above:
AvgStress = sum(odbOut1)/sum(odbOut2)
For any further inquiries you are encouraged to contact me here or through PM.
Best regards,
George