ich bin gerade dabei, ein python Skript für meine Optimierungsschleife zu schreiben. Nur ist das natürlich nicht ganz ohne Probleme
Das Problem ist in der 104. Zeile:
Zwar werden die Flächen, wenn ich es mir anschaue, richtig gewählt aber erst wenn ich sie allerdings editiere und erneut wähle, rechnet der Solver....
Ich tippe auf eine falsche Instanz. Aber egal was ich probiert habe=>ohne Erfolg.
Ich freue mich über jeden Tipp..
Viele Grüße aus Stuttgart
Da ich das Skript nicht hochladen kann, füge ich hier ein:
"""
Abaqus Skript
Rohrbogen
"""
from abaqus import *
from abaqusConstants import *
backwardCompatibility.setValues(includeDeprecated=True, reportDeprecated=False)
# Create a model
komplett = mdb.Model(name='Model-1')
mdb.saveAs('c:/Dokumente und Einstellungen/ykosarc/Desktop/Abaqus Scripting/model.cae');
#Import part
import part
acis = mdb.openAcis('c:/Dokumente und Einstellungen/ykosarc/Desktop/Abaqus Scripting/import.sat',
scaleFromFile=OFF)
komplett.PartFromGeometryFile(name='rohrbogen', geometryFile=acis,
dimensionality=THREE_D, type=DEFORMABLE_BODY)
rohrbogen = komplett.parts['rohrbogen']
#Ignore the imported edges
import mesh
rohrbogen.ignoreEntity((rohrbogen.edges[1],rohrbogen.edges[5]))
import part
#Define Datum plane
datum_xz=rohrbogen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE,offset=0)
#Partition the cell
main_cell=rohrbogen.cells[0]
datum_plane=rohrbogen.datums[3]
rohrbogen.PartitionCellByDatumPlane(main_cell,datum_plane)
# Create the material.
import material
Stahl = komplett.Material(name='Steel')
# Create the elastic properties: youngsModulus is 209.E3
# and poissonsRatio is 0.3
elasticProperties = (210.E3, 0.3)
densityProperties =(7800.E-9, 0)
Stahl.Elastic(table=(elasticProperties, ) )
Stahl.Density(table=(densityProperties, ) )
import section
# Create a section.
mySection = komplett.HomogeneousSolidSection(name='SteelSection',material='Steel')
# Assign the section to the region. The region refers
# to the single cell in this model.
region = (rohrbogen.cells,)
rohrbogen.SectionAssignment(region=region,sectionName='SteelSection')
#Mesh the part
import mesh
# Assign an element type to the part instance.
region = (rohrbogen.cells,)
elemType = mesh.ElemType(elemCode=C3D8R, elemLibrary=STANDARD)
rohrbogen.setElementType(regions=region, elemTypes=(elemType,))
# Seed the part instance.
rohrbogen.seedPart(size=13)
# Mesh the part instance.
rohrbogen.generateMesh()
# Create assembly
import assembly
myAssembly = komplett.rootAssembly
myInstance = myAssembly.Instance(name='Instance',
part=rohrbogen,dependent=ON)
#Create a step
import step
komplett.StaticStep(name='Statisch', previous='Initial',
timePeriod=1.0, initialInc=0.1,
description='Statische Belastung')
#Create loads and BC's
import load
# Create a boundary condition
bc_surface=(myInstance.faces[4],myInstance.faces[7])
komplett.EncastreBC(name='Fixed',createStepName='Statisch',region=bc_surface)
# Create a surface traction
load_surface=(myInstance.faces[2],myInstance.faces[9])
#komplett.SurfaceTraction(name='SurfaceTrac',createStepName='Statisch',region=load_surface,magnitude=1000,directionVector=((0,0,0),(0,0,-1)))
komplett.Pressure(name='Pressure',createStepName='Statisch',region=load_surface,magnitude=10)
import job
# Create an analysis job for the model and submit it.
jobName = 'Job_scripted'
myJob = mdb.Job(name=jobName, model=komplett,
description='Process the job',numCpus=2,numDomains=2)
# Wait for the job to complete.
myJob.submit()
myJob.waitForCompletion()
[Diese Nachricht wurde von solocan am 05. Jun. 2012 editiert.]
[Diese Nachricht wurde von solocan am 05. Jun. 2012 editiert.]
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP