×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

Hi, I have the following script,

Hi, I have the following script,

Hi, I have the following script,

(OP)
Hi,

I have the following script,which generates an ODB File.
My model is a simple Structure,consisting of circular beams.

When I am running the above script,its generating an ODB of size 10 GB.
From the ODB, I only need to retrieve Vonmisses stress, Displacement.

Script:

# -*- coding: mbcs -*-
from part import *
from material import *
from section import *
from assembly import *
from step import *
from interaction import *
from load import *
from mesh import *
from optimization import *
from job import *
from sketch import *
from visualization import *
from connectorBehavior import *
import sys
import numpy
index=0
parameter=[[]]
for x in range(1,101):
radius_beam=0.21+float(x)*0.003
mdb.models['Model-1'].ConstrainedSketch(name='__profile__', sheetSize=200.0)
mdb.models['Model-1'].sketches['__profile__'].rectangle(point1=(-3.0, 3.0),
point2=(3.0, -3.0))
mdb.models['Model-1'].Part(dimensionality=THREE_D, name='Part-1', type=
DEFORMABLE_BODY)
mdb.models['Model-1'].parts['Part-1'].BaseSolidExtrude(depth=8.25, sketch=
mdb.models['Model-1'].sketches['__profile__'])
del mdb.models['Model-1'].sketches['__profile__']

mdb.models['Model-1'].parts['Part-1'].DatumPlaneByOffset(flip=SIDE1, offset=2.5
, plane=mdb.models['Model-1'].parts['Part-1'].faces[5])
mdb.models['Model-1'].parts['Part-1'].Mirror(keepOriginal=ON, mirrorPlane=
mdb.models['Model-1'].parts['Part-1'].datums[2])

mdb.models['Model-1'].parts['Part-1'].DatumPlaneByOffset(flip=SIDE2, offset=1.0
, plane=mdb.models['Model-1'].parts['Part-1'].faces[6])
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=1.27, name='__sweep__',
sheetSize=50.85, transform=
mdb.models['Model-1'].parts['Part-1'].MakeSketchTransform(
sketchPlane=mdb.models['Model-1'].parts['Part-1'].datums[4],
sketchPlaneSide=SIDE1,
sketchUpEdge=mdb.models['Model-1'].parts['Part-1'].edges[14],
sketchOrientation=TOP, origin=(2.0, 0.0, -2.5)))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__sweep__'])
mdb.models['Model-1'].sketches['__sweep__'].Line(point1=(-2.0, 2.5), point2=(
-2.0, -2.5))
mdb.models['Model-1'].sketches['__sweep__'].VerticalConstraint(addUndoState=
False, entity=mdb.models['Model-1'].sketches['__sweep__'].geometry[2])
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=0.49, name='__profile__',
sheetSize=19.69, transform=
mdb.models['Model-1'].parts['Part-1'].MakeSketchTransform(
sketchPlane=mdb.models['Model-1'].parts['Part-1'].faces[11],
sketchPlaneSide=SIDE1,
sketchUpEdge=mdb.models['Model-1'].parts['Part-1'].edges[17],
sketchOrientation=BOTTOM, origin=(0.0, 0.0, -5.0)))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=(
2.0, -2.0), point1=(2.3275, -2.205))
#####################first radius parameter##########################################


print (radius_beam)
mdb.models['Model-1'].sketches['__profile__'].RadialDimension(curve=
mdb.models['Model-1'].sketches['__profile__'].geometry[6], radius=radius_beam,
textPoint=(3.0, -1.29795169830322))
####################################################################################
mdb.models['Model-1'].parts['Part-1'].SolidSweep(path=
mdb.models['Model-1'].sketches['__sweep__'], pathOrientation=TOP,
pathPlane=mdb.models['Model-1'].parts['Part-1'].datums[4], pathUpEdge=
mdb.models['Model-1'].parts['Part-1'].edges[14], profile=
mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=BOTTOM,
sketchPlane=mdb.models['Model-1'].parts['Part-1'].faces[11], sketchUpEdge=
mdb.models['Model-1'].parts['Part-1'].edges[17])
del mdb.models['Model-1'].sketches['__sweep__']
del mdb.models['Model-1'].sketches['__profile__']

###############rounding of edges of beam with 0.2###################################
mdb.models['Model-1'].parts['Part-1'].Round(edgeList=(
mdb.models['Model-1'].parts['Part-1'].edges[0],
mdb.models['Model-1'].parts['Part-1'].edges[1]), radius=0.2)
mdb.models['Model-1'].parts['Part-1'].DatumPlaneByOffset(flip=SIDE2, offset=3.0
, plane=mdb.models['Model-1'].parts['Part-1'].faces[5])
mdb.models['Model-1'].parts['Part-1'].DatumPlaneByOffset(flip=SIDE2, offset=3.0
, plane=mdb.models['Model-1'].parts['Part-1'].faces[13])
mdb.models['Model-1'].parts['Part-1'].Mirror(keepOriginal=ON, mirrorPlane=
mdb.models['Model-1'].parts['Part-1'].datums[8])
mdb.models['Model-1'].parts['Part-1'].Mirror(keepOriginal=ON, mirrorPlane=
mdb.models['Model-1'].parts['Part-1'].datums[7])

#####################construction of slanted beams#################################
mdb.models['Model-1'].parts['Part-1'].DatumAxisByCylFace(face=
mdb.models['Model-1'].parts['Part-1'].faces[9])
mdb.models['Model-1'].parts['Part-1'].DatumPointByProjOnEdge(edge=
mdb.models['Model-1'].parts['Part-1'].datums[11], point=
mdb.models['Model-1'].parts['Part-1'].vertices[16])
mdb.models['Model-1'].parts['Part-1'].DatumAxisByCylFace(face=
mdb.models['Model-1'].parts['Part-1'].faces[12])
mdb.models['Model-1'].parts['Part-1'].DatumPointByProjOnEdge(edge=
mdb.models['Model-1'].parts['Part-1'].edges[23], point=
mdb.models['Model-1'].parts['Part-1'].vertices[17])
del mdb.models['Model-1'].parts['Part-1'].features['Datum pt-2']
mdb.models['Model-1'].parts['Part-1'].DatumPointByProjOnEdge(edge=
mdb.models['Model-1'].parts['Part-1'].datums[13], point=
mdb.models['Model-1'].parts['Part-1'].vertices[17])
mdb.models['Model-1'].parts['Part-1'].DatumAxisByCylFace(face=
mdb.models['Model-1'].parts['Part-1'].faces[1])
mdb.models['Model-1'].parts['Part-1'].DatumPointByProjOnEdge(edge=
mdb.models['Model-1'].parts['Part-1'].datums[16], point=
mdb.models['Model-1'].parts['Part-1'].vertices[10])
del mdb.models['Model-1'].parts['Part-1'].features['Datum pt-2']
del mdb.models['Model-1'].parts['Part-1'].features['Datum axis-2']
mdb.models['Model-1'].parts['Part-1'].DatumAxisByTwoPoint(point1=
mdb.models['Model-1'].parts['Part-1'].datums[17], point2=
mdb.models['Model-1'].parts['Part-1'].datums[12])
mdb.models['Model-1'].parts['Part-1'].DatumPlaneByPointNormal(normal=
mdb.models['Model-1'].parts['Part-1'].datums[18], point=
mdb.models['Model-1'].parts['Part-1'].datums[12])
mdb.models['Model-1'].parts['Part-1'].DatumPlaneByOffset(flip=SIDE2, offset=1.0
, plane=mdb.models['Model-1'].parts['Part-1'].faces[17])
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=1.59, name='__sweep__',
sheetSize=63.7, transform=
mdb.models['Model-1'].parts['Part-1'].MakeSketchTransform(
sketchPlane=mdb.models['Model-1'].parts['Part-1'].datums[20],
sketchPlaneSide=SIDE1,
sketchUpEdge=mdb.models['Model-1'].parts['Part-1'].edges[9],
sketchOrientation=TOP, origin=(0.0, 2.0, -2.5)))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__sweep__'])
mdb.models['Model-1'].sketches['__sweep__'].Line(point1=(-2.0, -2.5), point2=(
2.0, 2.5))
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=1.59, name='__profile__',
sheetSize=63.7, transform=
mdb.models['Model-1'].parts['Part-1'].MakeSketchTransform(
sketchPlane=mdb.models['Model-1'].parts['Part-1'].datums[19],
sketchPlaneSide=SIDE1,
sketchUpEdge=mdb.models['Model-1'].parts['Part-1'].edges[16],
sketchOrientation=BOTTOM, origin=(2.0, 0.0, 0.0)))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=(
2.0, 0.0), point1=(2.0, -1.60398920346051))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
addUndoState=False, entity1=
mdb.models['Model-1'].sketches['__profile__'].vertices[3], entity2=
mdb.models['Model-1'].sketches['__profile__'].geometry[2])

###################### 2nd Parameter#################################
mdb.models['Model-1'].sketches['__profile__'].RadialDimension(curve=
mdb.models['Model-1'].sketches['__profile__'].geometry[4], radius= radius_beam,
textPoint=(7.3201789855957, 2.236790697694))
#####################################################################

mdb.models['Model-1'].parts['Part-1'].SolidSweep(path=
mdb.models['Model-1'].sketches['__sweep__'], pathOrientation=TOP,
pathPlane=mdb.models['Model-1'].parts['Part-1'].datums[20], pathUpEdge=
mdb.models['Model-1'].parts['Part-1'].edges[9], profile=
mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=BOTTOM,
sketchPlane=mdb.models['Model-1'].parts['Part-1'].datums[19], sketchUpEdge=
mdb.models['Model-1'].parts['Part-1'].edges[16])
del mdb.models['Model-1'].sketches['__sweep__']
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].parts['Part-1'].DatumAxisByCylFace(face=
mdb.models['Model-1'].parts['Part-1'].faces[13])
mdb.models['Model-1'].parts['Part-1'].DatumPointByProjOnEdge(edge=
mdb.models['Model-1'].parts['Part-1'].datums[22], point=
mdb.models['Model-1'].parts['Part-1'].vertices[20])
mdb.models['Model-1'].parts['Part-1'].DatumPointByProjOnEdge(edge=
mdb.models['Model-1'].parts['Part-1'].datums[16], point=
mdb.models['Model-1'].parts['Part-1'].vertices[22])
mdb.models['Model-1'].parts['Part-1'].DatumAxisByTwoPoint(point1=
mdb.models['Model-1'].parts['Part-1'].datums[12], point2=
mdb.models['Model-1'].parts['Part-1'].datums[23])
del mdb.models['Model-1'].parts['Part-1'].features['Datum pt-4']
mdb.models['Model-1'].parts['Part-1'].DatumPlaneByPointNormal(normal=
mdb.models['Model-1'].parts['Part-1'].datums[25], point=
mdb.models['Model-1'].parts['Part-1'].datums[12])
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=1.78, name='__sweep__',
sheetSize=71.45, transform=
mdb.models['Model-1'].parts['Part-1'].MakeSketchTransform(
sketchPlane=mdb.models['Model-1'].parts['Part-1'].datums[4],
sketchPlaneSide=SIDE1,
sketchUpEdge=mdb.models['Model-1'].parts['Part-1'].edges[22],
sketchOrientation=TOP, origin=(2.0, 0.0, -2.5)))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__sweep__'])
mdb.models['Model-1'].sketches['__sweep__'].Line(point1=(-2.0, 2.5), point2=(
2.0, -2.5))
mdb.models['Model-1'].ConstrainedSketch(gridSpacing=1.78, name='__profile__',
sheetSize=71.45, transform=
mdb.models['Model-1'].parts['Part-1'].MakeSketchTransform(
sketchPlane=mdb.models['Model-1'].parts['Part-1'].datums[26],
sketchPlaneSide=SIDE1,
sketchUpEdge=mdb.models['Model-1'].parts['Part-1'].edges[24],
sketchOrientation=TOP, origin=(0.0, 2.0, 0.0)))
mdb.models['Model-1'].parts['Part-1'].projectReferencesOntoSketch(filter=
COPLANAR_EDGES, sketch=mdb.models['Model-1'].sketches['__profile__'])
mdb.models['Model-1'].sketches['__profile__'].CircleByCenterPerimeter(center=(
-2.0, 0.0), point1=(-1.31586981890723, 0.534215919906273))
mdb.models['Model-1'].sketches['__profile__'].CoincidentConstraint(
addUndoState=False, entity1=
mdb.models['Model-1'].sketches['__profile__'].vertices[4], entity2=
mdb.models['Model-1'].sketches['__profile__'].geometry[3])
#######################3rd radius Parameter####################################
mdb.models['Model-1'].sketches['__profile__'].RadialDimension(curve=
mdb.models['Model-1'].sketches['__profile__'].geometry[6], radius=radius_beam,
textPoint=(-7.5357780456543, 1.22744603251194))
################################################################################

mdb.models['Model-1'].parts['Part-1'].SolidSweep(flipSweepDirection=ON, path=
mdb.models['Model-1'].sketches['__sweep__'], pathOrientation=TOP,
pathPlane=mdb.models['Model-1'].parts['Part-1'].datums[4], pathUpEdge=
mdb.models['Model-1'].parts['Part-1'].edges[22], profile=
mdb.models['Model-1'].sketches['__profile__'], sketchOrientation=TOP,
sketchPlane=mdb.models['Model-1'].parts['Part-1'].datums[26], sketchUpEdge=
mdb.models['Model-1'].parts['Part-1'].edges[24])
del mdb.models['Model-1'].sketches['__sweep__']
del mdb.models['Model-1'].sketches['__profile__']
mdb.models['Model-1'].parts['Part-1'].Mirror(keepOriginal=ON, mirrorPlane=
mdb.models['Model-1'].parts['Part-1'].datums[2])
mdb.models['Model-1'].parts['Part-1'].Mirror(keepOriginal=ON, mirrorPlane=
mdb.models['Model-1'].parts['Part-1'].datums[8])
mdb.models['Model-1'].parts['Part-1'].Mirror(keepOriginal=ON, mirrorPlane=
mdb.models['Model-1'].parts['Part-1'].datums[7])

#################Partitioning for meshing###################################

#Partitioning of top and bottom blocks
mdb.models['Model-1'].parts['Part-1'].PartitionCellByPlaneThreePoints(cells=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask(('[#1 ]',
), ), point1=mdb.models['Model-1'].parts['Part-1'].vertices[70], point2=
mdb.models['Model-1'].parts['Part-1'].vertices[72], point3=
mdb.models['Model-1'].parts['Part-1'].vertices[71])
mdb.models['Model-1'].parts['Part-1'].PartitionCellByPlaneThreePoints(cells=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask(('[#2 ]',
), ), point1=mdb.models['Model-1'].parts['Part-1'].vertices[66], point2=
mdb.models['Model-1'].parts['Part-1'].vertices[69], point3=
mdb.models['Model-1'].parts['Part-1'].vertices[67])

#Partitioning of Beams with radius of beam as parameter
mdb.models['Model-1'].parts['Part-1'].DatumPlaneByOffset(flip=SIDE1, offset=radius_beam*2.9
, plane=mdb.models['Model-1'].parts['Part-1'].faces[36])
mdb.models['Model-1'].parts['Part-1'].DatumPlaneByOffset(flip=SIDE1, offset=radius_beam*2.9
, plane=mdb.models['Model-1'].parts['Part-1'].faces[37])
mdb.models['Model-1'].parts['Part-1'].PartitionCellByDatumPlane(cells=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask(('[#4 ]',
), ), datumPlane=mdb.models['Model-1'].parts['Part-1'].datums[33])
mdb.models['Model-1'].parts['Part-1'].PartitionCellByDatumPlane(cells=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask(('[#40 ]',
), ), datumPlane=mdb.models['Model-1'].parts['Part-1'].datums[34])
mdb.models['Model-1'].parts['Part-1'].PartitionCellByDatumPlane(cells=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#20121 ]', ), ), datumPlane=
mdb.models['Model-1'].parts['Part-1'].datums[2])
mdb.models['Model-1'].parts['Part-1'].PartitionCellByDatumPlane(cells=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#200019 ]', ), ), datumPlane=
mdb.models['Model-1'].parts['Part-1'].datums[7])
print(mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#12060 ]', ), ))
mdb.models['Model-1'].parts['Part-1'].PartitionCellByDatumPlane(cells=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#12060 ]', ), ), datumPlane=
mdb.models['Model-1'].parts['Part-1'].datums[8])

#####################Create sets for defining element type ##############################

mdb.models['Model-1'].parts['Part-1'].Set(cells=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#38573fff ]', ), ), name='Hexahedral_set')
mdb.models['Model-1'].parts['Part-1'].Set(cells=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#7a8c000 ]', ), ), name='Tetrahedral_set')
####################Assign Element Type#################################
##Linear hexahedral element without reduced integration
mdb.models['Model-1'].parts['Part-1'].setElementType(elemTypes=(ElemType(
elemCode=C3D8, elemLibrary=STANDARD, secondOrderAccuracy=OFF,
distortionControl=DEFAULT), ElemType(elemCode=C3D6, elemLibrary=STANDARD),
ElemType(elemCode=C3D4, elemLibrary=STANDARD)), regions=(
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#38573fff ]', ), ), ))
##Linear Tetrahedral element
mdb.models['Model-1'].parts['Part-1'].setElementType(elemTypes=(ElemType(
elemCode=C3D8R, elemLibrary=STANDARD, secondOrderAccuracy=OFF,
kinematicSplit=AVERAGE_STRAIN, hourglassControl=DEFAULT,
distortionControl=DEFAULT), ElemType(elemCode=C3D6, elemLibrary=STANDARD),
ElemType(elemCode=C3D4, elemLibrary=STANDARD, secondOrderAccuracy=OFF,
distortionControl=DEFAULT)), regions=(
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#7a8c000 ]', ), ), ))
####################Assign Mesh Control#################################
mdb.models['Model-1'].parts['Part-1'].setMeshControls(elemShape=TET, regions=
mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#7a8c000 ]', ), ), technique=FREE)
mdb.models['Model-1'].parts['Part-1'].setElementType(elemTypes=(ElemType(
elemCode=C3D20R, elemLibrary=STANDARD), ElemType(elemCode=C3D15,
elemLibrary=STANDARD), ElemType(elemCode=C3D10, elemLibrary=STANDARD)),
regions=(mdb.models['Model-1'].parts['Part-1'].cells.getSequenceFromMask((
'[#7a8c000 ]', ), ), ))
####################Assign Mesh seed Parameters ################################
mesh_size=0.1
mdb.models['Model-1'].parts['Part-1'].seedPart(deviationFactor=0.3,
minSizeFactor=0.1, size=mesh_size)
mdb.models['Model-1'].parts['Part-1'].generateMesh()
############################Material Description################################

Young_Mod=700000
Poisson_Ratio=0.3
mdb.models['Model-1'].Material(name='PLA')
mdb.models['Model-1'].materials['PLA'].Elastic(table=( (Young_Mod,Poisson_Ratio), ))

#########################Section Creation######################################
mdb.models['Model-1'].HomogeneousSolidSection(material='PLA', name=
'Continumm_Sec', thickness=None)
########################Section Assignment#####################################
mdb.models['Model-1'].parts['Part-1'].SectionAssignment(offset=0.0,
offsetField='', offsetType=MIDDLE_SURFACE, region=
mdb.models['Model-1'].parts['Part-1'].sets['Hexahedral_set'], sectionName=
'Continumm_Sec', thicknessAssignment=FROM_SECTION)
mdb.models['Model-1'].parts['Part-1'].SectionAssignment(offset=0.0,
offsetField='', offsetType=MIDDLE_SURFACE, region=
mdb.models['Model-1'].parts['Part-1'].sets['Tetrahedral_set'], sectionName=
'Continumm_Sec', thicknessAssignment=FROM_SECTION)
#########################Assembly Creation#####################################
mdb.models['Model-1'].rootAssembly.DatumCsysByDefault(CARTESIAN)
mdb.models['Model-1'].rootAssembly.Instance(dependent=ON, name='Part-1-1',
part=mdb.models['Model-1'].parts['Part-1'])
#########################Step Creation#####################################
mdb.models['Model-1'].StaticStep(initialInc=0.01, maxInc=33.0, minInc=0.00001,
name='Step-1', previous='Initial', timePeriod=33.0)
########################Fixed Boundary Condition#####################################

mdb.models['Model-1'].rootAssembly.Set(faces=
mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
('[#0:3 #8000 ]', ), ), name='Fixed_End_Set')
mdb.models['Model-1'].EncastreBC(createStepName='Step-1', localCsys=None, name=
'FIxed_End_BC', region=
mdb.models['Model-1'].rootAssembly.sets['Fixed_End_Set'])

########################Displacement Boundary Condition#####################################
mdb.models['Model-1'].rootAssembly.Set(faces=
mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
('[#0:3 #400 ]', ), ), name='Disp_End_Set')
mdb.models['Model-1'].VelocityBC(amplitude=UNSET, createStepName='Step-1',
distributionType=UNIFORM, fieldName='', localCsys=None, name='Disp_Cond',
region=mdb.models['Model-1'].rootAssembly.sets['Disp_End_Set'], v1=UNSET,
v2=UNSET, v3=-0.01, vr1=UNSET, vr2=UNSET, vr3=UNSET)
#####################################Create Sets For Parameter Extraction############################################
mdb.models['Model-1'].rootAssembly.Set(faces=
mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
('[#0:3 #1 ]', ), ), name='Straight_Beam')
mdb.models['Model-1'].rootAssembly.Set(cells=
mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].cells.getSequenceFromMask(
('[#8000 ]', ), ), name='Joint_1')
mdb.models['Model-1'].rootAssembly.Set(cells=
mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].cells.getSequenceFromMask(
('[#2000000 ]', ), ), name='Joint_2')
mdb.models['Model-1'].rootAssembly.Set(faces=
mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
('[#0:2 #1000000 ]', ), ), name='Cross_Beam1')
mdb.models['Model-1'].rootAssembly.Set(faces=
mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].faces.getSequenceFromMask(
('[#0:2 #10000000 ]', ), ), name='Cross_Beam2')
############################################Create Job################################################################
mdb.Job(atTime=None, contactPrint=OFF, description='', echoPrint=OFF,
explicitPrecision=SINGLE, getMemoryFromAnalysis=True, historyPrint=OFF,
memory=90, memoryUnits=PERCENTAGE, model='Model-1', modelPrint=OFF,
multiprocessingMode=DEFAULT, name='Job-1', nodalOutputPrecision=SINGLE,
numCpus=1, numGPUs=0, queue=None, resultsFormat=ODB, scratch='', type=
ANALYSIS, userSubroutine='', waitHours=0, waitMinutes=0)
mdb.jobs['Job-1'].submit(consistencyChecking=OFF)
mdb.jobs['Job-1'].waitForCompletion()
parameter.append([mesh_size,radius_beam])
import shutil
index=index+1
shutil.copy2('Job-1.odb','results1/result%s.odb'%index)
numpy.save('parameter1',parameter)


Can someone please let me know the changes to be done in my script,such that I obtain an ODB of size 500 MB at maximum.
From the ODB, I only need to retrieve Vonmisses stress, Displacement.
Where can I add few lines of script,so that my ODB File contains only the data of Vonmisses stress and Displacement.

Thanks in Advance!

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources