---

Hallo zusammen,

mit Abaqus/CAE bin ich seit einigen Monaten recht gut vertraut.
Allerdings stehe ich bei der Pythonprogrammierung sowie Subroutinenerstellung noch am Anfang.

Zu meinem aktuellen Problem:
Ich möchte einen Zylinder modellieren und diesem Surfaces zuweisen. Da es sich um einen axialsymm. Körper handelt, habe ich ein 2D Modell (Rectangle) erstellt und ihm bei der Netzerstellung DCAX4-Elemente zugewiesen.

Skript.py:

print ' '
print 'Erstellen der Geometrie'
print ' '

#***************************************************************************************
# Modellparameter
#
# Einheiten:
#    Laengeneinheit    [mm]
#***************************************************************************************
innenradius =  0.0
aussenradius = 10.0
hoehe = 80.0

elementanzahl = 20
verfeinerung  = 5.0

#***************************************************************************************
# Erstellung der Geometrie
#***************************************************************************************
mdb.models[modelName].ConstrainedSketch(name='__profile__', sheetSize=aussenradius)
mdb.models[modelName].sketches['__profile__'].sketchOptions.setValues(viewStyle=AXISYM)
mdb.models[modelName].sketches['__profile__'].ConstructionLine(point1=(0.0, -50.0), point2=(0.0, 150.0))
mdb.models[modelName].sketches['__profile__'].rectangle(point1=(innenradius, 0.0), point2=(aussenradius, hoehe))
mdb.models[modelName].Part(name=partName, dimensionality=AXISYMMETRIC, type=DEFORMABLE_BODY)
mdb.models[modelName].parts[partName].BaseShell(sketch= mdb.models[modelName].sketches['__profile__'])
del mdb.models[modelName].sketches['__profile__']

#***************************************************************************************
# Erstellung des Netzes
#
# Elementtypen axialsymmetrisch, MassDiffusionStep
# linear: DCAX4
# quadratisch: DCAX8

Edges kann ich rel. einfach via findAt zuweisen. Aber kann ich dem Körper nun surfaces (surface sets) zuordnen? könnte ich bspw die äußere Kante (Symmetrie --> Mantelfläche) als surface definieren?!

Ich hoffe ich konnte euch die Rahmenbedingungen hinreichend erklären und jemand hat einen Tipp...

Viele Grüße

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Du kannst Surfaces erstellen und wählst dabei in axialsymmetrischen Modellen die Edges des Körpers aus.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hmm... , keiner der eine Antwort, Frage oder einen Korrekturvorschlag hat?!

Über die rpy habe ich folgendes bekommen:

mdb.models['Model-1'].rootAssembly.Surface(name='Surf-Mantel', side1Edges=
    mdb.models['Model-1'].rootAssembly.instances['Zylinder'].edges.getSequenceFromMask(
    ('[#2 ]', ), ))

Nur möchte ich dies über Python ansprechen können.
Kann ich also einfach edges via FindAt definieren, ggf. ganze sets aus edges und diese dann
als surface definieren?

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Ahh, Sorry, habe deine Info nicht gelesen.
Teste es gleich mal aus und gebe eine Rückmeldung.
Danke!

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hi,

habe folgenden Code verwendet:

surfMantel  = mdb.models[modelName].rootAssembly.Surface(side1Edges=mdb.models[modelName].rootAssembly.instances[partName].edges.findAt(10.0,40.0,0.0), name='Surf-Mantel')

Ausgangspunkt war:
.
.
Erstellung - Sets
#***************************************************************************************
mdb.models[modelName].rootAssembly.Set(name='Set-All', faces=mdb.models[modelName].rootAssembly.instances[partName].faces.getSequenceFromMask(mask=('[#1 ]', ), ))

surfMantel  = mdb.models[modelName].rootAssembly.Surface(side1Edges=mdb.models[modelName].rootAssembly.instances[partName].edges.getSequenceFromMask(mask=('[#6 ]', ), ), name='Surf-Mantel)
.
.

Ich muss findAt doch immer in 3D angeben, oder?!

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Bei findAt ist axialsymmetrisch die z-Komponente null.

Auf Assembly-Ebene holst du dir eine Edge z.B. so:

Code:

---

x = mdb.models['Model-1'].rootAssembly.instances['Part-1-1'].edges.findAt((15.0, 20.0, 0))

---

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP