from abaqus import *
from abaqusConstants import *
session.Viewport(name='Viewport: 1', origin=(0.0, 0.0), width=227.098129272461, 
    height=234.260009765625)
session.viewports['Viewport: 1'].makeCurrent()
session.viewports['Viewport: 1'].maximize()
from caeModules import *
from driverUtils import executeOnCaeStartup
executeOnCaeStartup()
Mdb()
#: A new model database has been created.
#: The model "Model-1" has been created.
session.viewports['Viewport: 1'].setValues(displayedObject=None)
#############################################################################################
#                     einige Variablen
#############################################################################################
l1 =     2.5000000
l2 =     2.0000000
l3 =     3.0000000

elsize=0.5
#############################################################################################
#                     Skizze wird erstellt
#############################################################################################
p = mdb.models['Model-1'].Part(name='Part-1', dimensionality=THREE_D, 
    type=DEFORMABLE_BODY)
p.ReferencePoint(point=(0.0, 0.0, 0.0))
p = mdb.models['Model-1'].parts['Part-1']

p.DatumPointByCoordinate(coords=(    -2.5000000,    -2.0000000,     0.0000000))
p.DatumPointByCoordinate(coords=(     0.5000000,    -2.0000000,     0.0000000))
p.DatumPointByCoordinate(coords=(     1.5000000,     1.0000000,     0.0000000))
p.DatumPointByCoordinate(coords=(    -1.5000000,     1.0000000,     0.0000000))
p.DatumPointByCoordinate(coords=(    -1.5000000,    -1.0000000,     3.0000000))
p.DatumPointByCoordinate(coords=(     1.5000000,    -1.0000000,     3.0000000))
p.DatumPointByCoordinate(coords=(     2.5000000,     2.0000000,     3.0000000))
p.DatumPointByCoordinate(coords=(    -0.5000000,     2.0000000,     3.0000000))

#
d1 = p.datums
p.DatumPlaneByThreePoints(point1=d1[2], point2=d1[3], point3=d1[6])
p.DatumPlaneByThreePoints(point1=d1[5], point2=d1[4], point3=d1[9])
p.DatumPlaneByThreePoints(point1=d1[2], point2=d1[5], point3=d1[6])
p.DatumPlaneByThreePoints(point1=d1[3], point2=d1[4], point3=d1[7])
p.DatumPlaneByThreePoints(point1=d1[2], point2=d1[3], point3=d1[4])
p.DatumPlaneByThreePoints(point1=d1[6], point2=d1[7], point3=d1[8])
session.viewports['Viewport: 1'].view.fitView()
p = mdb.models['Model-1'].parts['Part-1']
p.DatumAxisByPrincipalAxis(principalAxis=YAXIS)
p.DatumAxisByPrincipalAxis(principalAxis=XAXIS)
p.DatumAxisByPrincipalAxis(principalAxis=ZAXIS)
# klappt bis hierher
#############################################################################################
#                     Wire- Features
#############################################################################################
p = mdb.models['Model-1'].parts['Part-1']
d = p.datums
p.WirePolyLine(points=((d[2], d[3]), (d[3], d[4]), (d[4], d[5]), (d[5], d[2])), 
    mergeWire=OFF, meshable=ON)
p.WirePolyLine(points=((d[6], d[7]), (d[7], d[8]), (d[8], d[9]), (d[9], 
    d[6])), mergeWire=OFF, meshable=ON)
p.WirePolyLine(points=((d[3], d[2]), (d[2], d[6]), (d[6], d[7]), (d[7], d[3])), 
    mergeWire=OFF, meshable=ON)
p.WirePolyLine(points=((d[4], d[5]), (d[5], d[9]), (d[9], d[8]), (d[8], 
    d[4])), mergeWire=OFF, meshable=ON)
e = p.edges
p.AddFaces(edgeList = e[1:2]+e[6:7]+e[9:10]+e[11:12], tryAnalytical=True)
p.AddFaces(edgeList = e[1:5], tryAnalytical=True)
p.AddFaces(edgeList = e[0:2]+e[7:8]+e[10:11], tryAnalytical=True)
p.AddFaces(edgeList = e[1:3]+e[6:7]+e[11:12], tryAnalytical=True)
p.AddFaces(edgeList = e[0:1]+e[6:7]+e[9:10]+e[11:12], tryAnalytical=True)
p.AddFaces(edgeList = e[0:1]+e[4:5]+e[10:12], tryAnalytical=True)
f = p.faces
p.AddCells(faceList = f[0:6])
#
#############################################################################################
#                     Mesh, Assembly, Section, Material
#############################################################################################
p = mdb.models['Model-1'].parts['Part-1']
p.seedPart(size=elsize, deviationFactor=0.1)
p = mdb.models['Model-1'].parts['Part-1']
p.generateMesh()
#
mdb.models['Model-1'].Material(name='Material-1')
mdb.models['Model-1'].materials['Material-1'].Elastic(table=((210000.0, 0.3), 
    ))
mdb.models['Model-1'].HomogeneousSolidSection(name='Section-1', 
    material='Material-1', thickness=1.0)
p = mdb.models['Model-1'].parts['Part-1']
c = p.cells
cells = c.getSequenceFromMask(mask=('[#1 ]', ), )
region = regionToolset.Region(cells=cells)
p = mdb.models['Model-1'].parts['Part-1']
p.SectionAssignment(region=region, sectionName='Section-1', offset=0.0)
a = mdb.models['Model-1'].rootAssembly
session.viewports['Viewport: 1'].setValues(displayedObject=a)
a = mdb.models['Model-1'].rootAssembly
a.DatumCsysByDefault(CARTESIAN)
p = mdb.models['Model-1'].parts['Part-1']
a.Instance(name='Part-1-1', part=p, dependent=ON)
a = mdb.models['Model-1'].rootAssembly
a.translate(instanceList=('Part-1-1', ), vector=(0.0, 0.0, -l3/2))
#
#############################################################################################
#                     Speicherung
#############################################################################################
a = mdb.models['Model-1'].rootAssembly
session.viewports['Viewport: 1'].setValues(displayedObject=a)
session.viewports['Viewport: 1'].view.setProjection(projection=PARALLEL)
session.viewports['Viewport: 1'].assemblyDisplay.geometryOptions.setValues(
    datumAxes=OFF, datumPlanes=OFF)
mdb.saveAs(
    pathName='/home_mfk/nlfkm/branke/ABAQUS/Homogenisierung/Modell/SpattestSkript')
p = mdb.models['Model-1'].parts['Part-1']