gkapitanov
Bioengineer
- Aug 29, 2015
- 4
Hi,
A simple project - we have a rectangle, a load is evenly applied on a centered line on its top, and we have to record the point coordinates and the axial strain at those coordinates. Everything seems right, the abasus picture output is a symmetrical strain field, as expected(attached picture). However, when I extract the data points' coordinates and corresponding strains and plot them on matlab, the result is off (third post picture). I want data that resembles the abaqus output - what am I doing wrong? MATLAB output is not the problem - the problem is that the coordinates do not match the strains. The code I use for extracting the data is below. Thanks for any help.
from odbAccess import *
#open file to write on and the odb
f = open('/space/gkapitanov/Desktop/gkapitanov_ayati_lab/Pseudo_ABAQUS_Plot/Full_data2.txt', 'wab')
odb = openOdb(path='/space/gkapitanov/Desktop/gkapitanov_ayati_lab/Pseudo_ABAQUS_Plot/Blunt_result2.odb')
# set up the extract
pointData = odb.rootAssembly.instances['DISK-1'].nodes;
#last frame
modelFrame = odb.steps['Step-2'].frames[-1]
strain = modelFrame.fieldOutputs['E']
# define a variable for the values of the strain data
field Values = strain.values
for v in fieldValues:
labelVar = v.elementLabel
x_var = pointData[labelVar].coordinates[0]
y_var = pointData[labelVar].coordinates[1]
f.write('%6.9f %6.9f %6.9f\n' % (x_var, y_var, v.data[2]))
odb.close()
f.close()
A simple project - we have a rectangle, a load is evenly applied on a centered line on its top, and we have to record the point coordinates and the axial strain at those coordinates. Everything seems right, the abasus picture output is a symmetrical strain field, as expected(attached picture). However, when I extract the data points' coordinates and corresponding strains and plot them on matlab, the result is off (third post picture). I want data that resembles the abaqus output - what am I doing wrong? MATLAB output is not the problem - the problem is that the coordinates do not match the strains. The code I use for extracting the data is below. Thanks for any help.
from odbAccess import *
#open file to write on and the odb
f = open('/space/gkapitanov/Desktop/gkapitanov_ayati_lab/Pseudo_ABAQUS_Plot/Full_data2.txt', 'wab')
odb = openOdb(path='/space/gkapitanov/Desktop/gkapitanov_ayati_lab/Pseudo_ABAQUS_Plot/Blunt_result2.odb')
# set up the extract
pointData = odb.rootAssembly.instances['DISK-1'].nodes;
#last frame
modelFrame = odb.steps['Step-2'].frames[-1]
strain = modelFrame.fieldOutputs['E']
# define a variable for the values of the strain data
field Values = strain.values
for v in fieldValues:
labelVar = v.elementLabel
x_var = pointData[labelVar].coordinates[0]
y_var = pointData[labelVar].coordinates[1]
f.write('%6.9f %6.9f %6.9f\n' % (x_var, y_var, v.data[2]))
odb.close()
f.close()