Greeting
I want print to the Von Mises stress into a txt file using Abaqus script that I have. Attached is the script. I'm able to print the s values at different frames but the point is that I'm not able to print Von Mises stress. I think I should modify this line of the code:
field = frame.fieldOutputs['S'].getSubset(position=INTEGRATION_POINT, readOnly=False)
S = field.bulkDataBlocks[0].data
I tried different commands, but it did not work. So, I'm wondering if there is a way to print Von Mises stress using this script.
Thanks!
Sajjad
Hello Sajjad Izadpanah ,
To calculate the Von Mises stress, you will need to modify the script to extract the stress tensor components and then compute the Von Mises stress based on these components. You can then write the computed Von Mises stress values to a text file. Here is the modified script:
Add a function to compute the Von Mises stress:
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
def compute_von_mises_stress(stress_components):
s11, s22, s33, s12, s13, s23 = stress_components
von_mises_stress = np.sqrt(0.5 * ((s11 - s22) ** 2 + (s22 - s33) ** 2 + (s33 - s11) ** 2 + 6 * (s12 ** 2 + s13 ** 2 + s23 ** 2)))
return von_mises_stress
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Modify the script to extract the stress tensor components and compute the Von Mises stress:
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
field = frame.fieldOutputs['S'].getSubset(position=INTEGRATION_POINT, readOnly=False)
s11 = field.getScalarField(componentLabel='S11')
s22 = field.getScalarField(componentLabel='S22')
s33 = field.getScalarField(componentLabel='S33')
s12 = field.getScalarField(componentLabel='S12')
s13 = field.getScalarField(componentLabel='S13')
s23 = field.getScalarField(componentLabel='S23')
stress_components = [s11.values[0].data, s22.values[0].data, s33.values[0].data, s12.values[0].data, s13.values[0].data, s23.values[0].data]
von_mises_stress = compute_von_mises_stress(stress_components)
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Write the Von Mises stress values to a text file:
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
with open('von_mises_stress.txt', 'w') as f:
for stress_value in von_mises_stress:
f.write(str(stress_value) + '\n')
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Now, when you run the modified script, it should compute the Von Mises stress and write the values to a text file called von_mises_stress.txt.
Regards
- Badges
- Science topic
- Similar topics
- Biological Science
- Histology