######################################################################                                                       
# ABAQUS Benchmarks                                                  #
# Script to generate 2-dimensional model for non-symmetric multi-    #
# delamination analysis in ABAQUS/Standard.                          #
######################################################################
#
#*******************
# NOTE: ---->
#*******************
# Before submitting the job for analysis material for *COHESIVE SECTION
# should be edited to incorporate coupled (*ELASTIC, TYPE=COUPLED TRACTION)
# or uncoupled (*ELASTIC, TYPE=TRACTION) traction separation behavior. 
#
from abaqus import *
from abaqusConstants import *
session.Viewport(name='Viewport: 1', origin=(0.0, 0.0), width=261.9140625, 
    height=208.0078125)
session.viewports['Viewport: 1'].makeCurrent()
session.viewports['Viewport: 1'].maximize()
from caeModules import *
from driverUtils import executeOnCaeStartup
executeOnCaeStartup()
Mdb()
session.viewports['Viewport: 1'].setValues(displayedObject=None)
#
# Change model name to Alfano
#
mdb.models.changeKey(fromName='Model-1', toName='Alfano')
session.viewports['Viewport: 1'].setValues(displayedObject=None)
#
# a. Sketch and create Plate part
# b. Create profile of two cracks and the cohesive layers
# c. Remove material from the plate part by extruding two cuts (a. - b.)
#
s = mdb.models['Alfano'].Sketch(name='__profile__', sheetSize=0.4)
g, v, d = s.geometry, s.vertices, s.dimensions
s.setPrimaryObject(option=STANDALONE)
s.rectangle(point1=(0.0, 0.0), point2=(0.2, 0.0032065))
session.viewports['Viewport: 1'].view.fitView()
s.HorizontalDimension(vertex1=v[2], vertex2=v[6], textPoint=(
    0.0984431654214859, -0.00427813269197941))
s.VerticalDimension(vertex1=v[6], vertex2=v[4], textPoint=(0.202578246593475, 
    0.00281412387266755))
p = mdb.models['Alfano'].Part(name='Plate', dimensionality=TWO_D_PLANAR, 
    type=DEFORMABLE_BODY)
p = mdb.models['Alfano'].parts['Plate']
p.BaseShell(sketch=s)
s.unsetPrimaryObject()
p = mdb.models['Alfano'].parts['Plate']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
del mdb.models['Alfano'].sketches['__profile__']
s1 = mdb.models['Alfano'].Sketch(name='__profile__', sheetSize=0.4, 
    gridSpacing=0.01)
g, v, d = s1.geometry, s1.vertices, s1.dimensions
s1.setPrimaryObject(option=SUPERIMPOSE)
p = mdb.models['Alfano'].parts['Plate']
p.projectReferencesOntoSketch(sketch=s1, filter=COPLANAR_EDGES)
r, r1 = s1.referenceGeometry, s1.referenceVertices
s1.rectangle(point1=(0.0, 0.001325), point2=(0.2, 0.00133825))
s1.rectangle(point1=(0.0, 0.00160325), point2=(0.2, 0.0016165))
s1.VerticalDimension(vertex1=r1.findAt((0.0, 0.0), 1), vertex2=v[1], 
    textPoint=(-0.000283258472336456, 0.000579574669245631))
s1.VerticalDimension(vertex1=v[1], vertex2=v[2], textPoint=(
    -0.000283258472336456, 0.00133974268101156))
s1.VerticalDimension(vertex1=v[2], vertex2=v[9], textPoint=(
    -0.000282614288153127, 0.00158178864512593))
s1.VerticalDimension(vertex1=v[9], vertex2=v[10], textPoint=(
    -0.000280867592664436, 0.00161846936680377))
session.viewports['Viewport: 1'].view.fitView()
s1.HorizontalDimension(vertex1=v[2], vertex2=v[6], textPoint=(
    0.0364299640059471, -0.00248509692028165))
session.viewports['Viewport: 1'].view.fitView()
s1.HorizontalDimension(vertex1=v[9], vertex2=v[14], textPoint=(
    0.0737828761339188, -0.0102748423814774))
s1.constraintReferences(vertex1=r1.findAt((0.0, 0.0), 1))
p = mdb.models['Alfano'].parts['Plate']
p.Cut(sketch=s1)
s1.unsetPrimaryObject()
del mdb.models['Alfano'].sketches['__profile__']
#
# Create top and bottom cohesive layers with finite thickness
#
s = mdb.models['Alfano'].Sketch(name='__profile__', sheetSize=0.4)
g, v, d = s.geometry, s.vertices, s.dimensions
s.setPrimaryObject(option=STANDALONE)
s.rectangle(point1=(0.0, 0.0), point2=(0.2, 1.325e-05))
session.viewports['Viewport: 1'].view.fitView()
s.HorizontalDimension(vertex1=v[1], vertex2=v[6], textPoint=(0.199968859553337, 
    -1.13981359390891e-05))
s.VerticalDimension(vertex1=v[6], vertex2=v[4], textPoint=(0.20001257956028, 
    6.92647154210135e-06))
p = mdb.models['Alfano'].Part(name='Bottom_cohesive', 
    dimensionality=TWO_D_PLANAR, type=DEFORMABLE_BODY)
p = mdb.models['Alfano'].parts['Bottom_cohesive']
p.BaseShell(sketch=s)
s.unsetPrimaryObject()
p = mdb.models['Alfano'].parts['Bottom_cohesive']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
del mdb.models['Alfano'].sketches['__profile__']
p = mdb.models['Alfano'].Part(name='Top_cohesive', 
    objectToCopy=mdb.models['Alfano'].parts['Bottom_cohesive'])
session.viewports['Viewport: 1'].setValues(displayedObject=p)
s1 = mdb.models['Alfano'].Sketch(name='__profile__', sheetSize=0.4, 
    gridSpacing=0.01)
g, v, d = s1.geometry, s1.vertices, s1.dimensions
s1.setPrimaryObject(option=SUPERIMPOSE)
p = mdb.models['Alfano'].parts['Top_cohesive']
p.projectReferencesOntoSketch(sketch=s1, filter=COPLANAR_EDGES)
r, r1 = s1.referenceGeometry, s1.referenceVertices
s1.rectangle(point1=(0.0, 0.0), point2=(0.04, 1.325e-05))
session.viewports['Viewport: 1'].view.fitView()
s1.VerticalDimension(vertex1=v[4], vertex2=v[2], textPoint=(
    -0.0020265206694603, 0.000235539977438748))
s1.constraintReferences(vertex1=r1.findAt((0.0, 0.0), 1))
p = mdb.models['Alfano'].parts['Top_cohesive']
p.Cut(sketch=s1)
s1.unsetPrimaryObject()
del mdb.models['Alfano'].sketches['__profile__']
p = mdb.models['Alfano'].parts['Bottom_cohesive']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
s = mdb.models['Alfano'].Sketch(name='__profile__', sheetSize=0.4, 
    gridSpacing=0.01)
g, v, d = s.geometry, s.vertices, s.dimensions
s.setPrimaryObject(option=SUPERIMPOSE)
p = mdb.models['Alfano'].parts['Bottom_cohesive']
p.projectReferencesOntoSketch(sketch=s, filter=COPLANAR_EDGES)
r, r1 = s.referenceGeometry, s.referenceVertices
s.rectangle(point1=(0.06, 0.0), point2=(0.08, 1.325e-05))
s.VerticalDimension(vertex1=v[1], vertex2=v[2], textPoint=(
    -0.00223228335380554, 0.00554204545915127))
s.HorizontalDimension(vertex1=v[2], vertex2=v[6], textPoint=(
    0.0799874514341354, -7.78574940341059e-06))
session.viewports['Viewport: 1'].view.fitView()
p = mdb.models['Alfano'].parts['Bottom_cohesive']
p.Cut(sketch=s)
s.unsetPrimaryObject()
del mdb.models['Alfano'].sketches['__profile__']
#
# Create material Mat1 for the plate and Coh for the cohesive layer
# 
session.viewports['Viewport: 1'].partDisplay.setValues(sectionAssignments=ON, 
    engineeringFeatures=ON)
mdb.models['Alfano'].Material(name='Coh')
mdb.models['Alfano'].materials['Coh'].Elastic(table=((200.0, 0.3), ))
mdb.models['Alfano'].Material(name='Mat1')
mdb.models['Alfano'].materials['Mat1'].Density(table=((2038800.0, ), ))
mdb.models['Alfano'].materials['Mat1'].Elastic(type=ENGINEERING_CONSTANTS, 
    table=((115000000000.0, 8500000000.0, 8500000000.0, 0.29, 0.29, 0.3, 
    4500000000.0, 3300000000.0, 4500000000.0), ))
mdb.models['Alfano'].HomogeneousSolidSection(name='plate', material='Mat1', 
    thickness=0.02)
mdb.models['Alfano'].CohesiveSection(name='coh', material='Coh', 
    response=TRACTION_SEPARATION, initialThicknessType=GEOMETRY, 
    outOfPlaneThickness=0.02)
p = mdb.models['Alfano'].parts['Plate']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
p2 = mdb.models['Alfano'].parts['Plate']
p2.DatumCsysByThreePoints(name='Global', coordSysType=CARTESIAN, origin=(0.0, 
    0.0, 0.0), point1=(1.0, 0.0, 0.0), line2=(0.0, 1.0, 0.0))
f = p2.faces
faces = f[0:3]
region = regionToolset.Region(faces=faces)
p1 = mdb.models['Alfano'].parts['Plate']
datums = p1.datums[3]
p1.assignMaterialOrientation(region=region, localCsys=datums, axis=AXIS_3)
#: Specified material orientation has been assigned to the selected regions.
p1.SectionAssignment(region=region, sectionName='plate')
p = mdb.models['Alfano'].parts['Bottom_cohesive']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
p2 = mdb.models['Alfano'].parts['Bottom_cohesive']
f = p2.faces
faces = f[0:2]
region = regionToolset.Region(faces=faces)
p1 = mdb.models['Alfano'].parts['Bottom_cohesive']
p1.SectionAssignment(region=region, sectionName='coh')
p = mdb.models['Alfano'].parts['Top_cohesive']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
p2 = mdb.models['Alfano'].parts['Top_cohesive']
f = p2.faces
faces = f[0:1]
region = regionToolset.Region(faces=faces)
p1 = mdb.models['Alfano'].parts['Top_cohesive']
p1.SectionAssignment(region=region, sectionName='coh')
a = mdb.models['Alfano'].rootAssembly
session.viewports['Viewport: 1'].setValues(displayedObject=a)
a = mdb.models['Alfano'].rootAssembly
a.DatumCsysByDefault(CARTESIAN)
p = mdb.models['Alfano'].parts['Bottom_cohesive']
a.Instance(name='Bottom_cohesive-1', part=p, dependent=ON)
a = mdb.models['Alfano'].rootAssembly
p = mdb.models['Alfano'].parts['Plate']
a.Instance(name='Plate-1', part=p, dependent=ON)
p3 = a.instances['Plate-1']
p3.translate(vector=(0.22, 0.0, 0.0))
session.viewports['Viewport: 1'].view.fitView()
a = mdb.models['Alfano'].rootAssembly
p = mdb.models['Alfano'].parts['Top_cohesive']
a.Instance(name='Top_cohesive-1', part=p, dependent=ON)
p3 = a.instances['Top_cohesive-1']
p3.translate(vector=(0.396, 0.0, 0.0))
session.viewports['Viewport: 1'].view.fitView()
a = mdb.models['Alfano'].rootAssembly
v11 = a.instances['Top_cohesive-1'].vertices
v12 = a.instances['Plate-1'].vertices
a.CoincidentPoint(movablePoint=v11[2], fixedPoint=v12[10])
session.viewports['Viewport: 1'].view.fitView()
a = mdb.models['Alfano'].rootAssembly
v11 = a.instances['Bottom_cohesive-1'].vertices
v12 = a.instances['Plate-1'].vertices
a.CoincidentPoint(movablePoint=v11[2], fixedPoint=v12[4])
session.viewports['Viewport: 1'].view.fitView()
session.viewports['Viewport: 1'].assemblyDisplay.setValues(visibleInstances=(
    'Plate-1', ))
a = mdb.models['Alfano'].rootAssembly
s1 = a.instances['Plate-1'].edges
side1Edges1 = s1[0:1]
#
# Create surfaces for defining *TIE and *CONTACT PAIRS
#
a.Surface(side1Edges=side1Edges1, name='top')
a = mdb.models['Alfano'].rootAssembly
s1 = a.instances['Plate-1'].edges
side1Edges1 = s1[10:11]
a.Surface(side1Edges=side1Edges1, name='middle-top')
a = mdb.models['Alfano'].rootAssembly
s1 = a.instances['Plate-1'].edges
side1Edges1 = s1[8:9]
a.Surface(side1Edges=side1Edges1, name='middle-bot')
a = mdb.models['Alfano'].rootAssembly
s1 = a.instances['Plate-1'].edges
side1Edges1 = s1[4:5]
a.Surface(side1Edges=side1Edges1, name='bottom')
session.viewports['Viewport: 1'].assemblyDisplay.setValues(visibleInstances=(
    'Bottom_cohesive-1', ))
session.viewports['Viewport: 1'].view.fitView()
a = mdb.models['Alfano'].rootAssembly
s1 = a.instances['Bottom_cohesive-1'].edges
side1Edges1 = s1[3:4]+s1[5:6]
a.Surface(side1Edges=side1Edges1, name='botCoh-Top')
a = mdb.models['Alfano'].rootAssembly
s1 = a.instances['Bottom_cohesive-1'].edges
side1Edges1 = s1[1:2]+s1[7:8]
a.Surface(side1Edges=side1Edges1, name='botCoh-Bot')
session.viewports['Viewport: 1'].assemblyDisplay.setValues(visibleInstances=(
    'Top_cohesive-1', ))
a = mdb.models['Alfano'].rootAssembly
s1 = a.instances['Top_cohesive-1'].edges
side1Edges1 = s1[3:4]
a.Surface(side1Edges=side1Edges1, name='topCoh-Top')
a = mdb.models['Alfano'].rootAssembly
s1 = a.instances['Top_cohesive-1'].edges
side1Edges1 = s1[1:2]
a.Surface(side1Edges=side1Edges1, name='topCoh-Bot')
session.viewports['Viewport: 1'].assemblyDisplay.setValues(visibleInstances=(
    'Bottom_cohesive-1', 'Plate-1', 'Top_cohesive-1'))
session.viewports['Viewport: 1'].view.fitView()
#
# Mesh the respective parts
#
session.viewports['Viewport: 1'].assemblyDisplay.setValues(mesh=ON)
session.viewports['Viewport: 1'].assemblyDisplay.meshOptions.setValues(
    meshTechnique=ON)
p = mdb.models['Alfano'].parts['Top_cohesive']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
session.viewports['Viewport: 1'].partDisplay.setValues(sectionAssignments=OFF, 
    engineeringFeatures=OFF, mesh=ON)
session.viewports['Viewport: 1'].partDisplay.meshOptions.setValues(
    meshTechnique=ON)
p = mdb.models['Alfano'].parts['Plate']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
p2 = mdb.models['Alfano'].parts['Plate']
f = p2.faces
pickedRegions = f[0:3]
p2.setMeshControls(regions=pickedRegions, technique=STRUCTURED)
p2 = mdb.models['Alfano'].parts['Plate']
e = p2.edges
pickedEdges = e[0:1]+e[2:3]+e[4:5]+e[6:7]+e[8:9]+e[10:11]
p2.seedEdgeByNumber(edges=pickedEdges, number=200, constraint=FIXED)
session.viewports['Viewport: 1'].view.fitView()
p2 = mdb.models['Alfano'].parts['Plate']
e = p2.edges
pickedEdges = e[1:2]+e[3:4]+e[5:6]+e[7:8]+e[9:10]+e[11:12]
p2.seedEdgeByNumber(edges=pickedEdges, number=1, constraint=FIXED)
session.viewports['Viewport: 1'].view.fitView()
elemType1 = mesh.ElemType(elemCode=CPE4I, elemLibrary=STANDARD)
elemType2 = mesh.ElemType(elemCode=CPE3, elemLibrary=STANDARD)
p3 = mdb.models['Alfano'].parts['Plate']
f = p3.faces
faces = f[0:3]
pickedRegions =(faces, )
p3.setElementType(regions=pickedRegions, elemTypes=(elemType1, elemType2))
p3 = mdb.models['Alfano'].parts['Plate']
p3.generateMesh()
#
# Only swept meshing technique can be used for cohesive layers
#
p = mdb.models['Alfano'].parts['Top_cohesive']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
p2 = mdb.models['Alfano'].parts['Top_cohesive']
f = p2.faces
pickedRegions = f[0:1]
p2.setMeshControls(regions=pickedRegions, technique=SWEEP)
p3 = mdb.models['Alfano'].parts['Top_cohesive']
f, e = p3.faces, p3.edges
p3.setSweepPath(region=f[0], edge=e[0], sense=REVERSE)
p2 = mdb.models['Alfano'].parts['Top_cohesive']
e = p2.edges
pickedEdges = e[0:1]
p2.seedEdgeByNumber(edges=pickedEdges, number=1, constraint=FIXED)
p2 = mdb.models['Alfano'].parts['Top_cohesive']
e = p2.edges
pickedEdges = e[1:2]
p2.seedEdgeBySize(edges=pickedEdges, size=0.001, constraint=FIXED)
elemType1 = mesh.ElemType(elemCode=COH2D4, elemLibrary=STANDARD)
elemType2 = mesh.ElemType(elemCode=UNKNOWN_TRI, elemLibrary=STANDARD)
p3 = mdb.models['Alfano'].parts['Top_cohesive']
f = p3.faces
faces = f[0:1]
pickedRegions =(faces, )
p3.setElementType(regions=pickedRegions, elemTypes=(elemType1, elemType2))
p3 = mdb.models['Alfano'].parts['Top_cohesive']
p3.generateMesh()
session.viewports['Viewport: 1'].view.fitView()
p = mdb.models['Alfano'].parts['Bottom_cohesive']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
session.viewports['Viewport: 1'].view.fitView()
p2 = mdb.models['Alfano'].parts['Bottom_cohesive']
f = p2.faces
pickedRegions = f[0:2]
p2.setMeshControls(regions=pickedRegions, technique=SWEEP)
p3 = mdb.models['Alfano'].parts['Bottom_cohesive']
f1, e1 = p3.faces, p3.edges
p3.setSweepPath(region=f1[0], edge=e1[0], sense=REVERSE)
p3 = mdb.models['Alfano'].parts['Bottom_cohesive']
f, e = p3.faces, p3.edges
p3.setSweepPath(region=f[1], edge=e[6], sense=REVERSE)
session.viewports['Viewport: 1'].view.fitView()
p2 = mdb.models['Alfano'].parts['Bottom_cohesive']
e = p2.edges
pickedEdges = e[0:1]+e[6:7]
p2.seedEdgeByNumber(edges=pickedEdges, number=1, constraint=FIXED)
session.viewports['Viewport: 1'].view.fitView()
p2 = mdb.models['Alfano'].parts['Bottom_cohesive']
e = p2.edges
pickedEdges = e[1:2]+e[7:8]
p2.seedEdgeBySize(edges=pickedEdges, size=0.001, constraint=FIXED)
session.viewports['Viewport: 1'].view.fitView()
elemType1 = mesh.ElemType(elemCode=COH2D4, elemLibrary=STANDARD)
elemType2 = mesh.ElemType(elemCode=UNKNOWN_TRI, elemLibrary=STANDARD)
p3 = mdb.models['Alfano'].parts['Bottom_cohesive']
f = p3.faces
faces = f[0:2]
pickedRegions =(faces, )
p3.setElementType(regions=pickedRegions, elemTypes=(elemType1, elemType2))
p3 = mdb.models['Alfano'].parts['Bottom_cohesive']
p3.generateMesh()
a = mdb.models['Alfano'].rootAssembly
session.viewports['Viewport: 1'].setValues(displayedObject=a)
a2 = mdb.models['Alfano'].rootAssembly
a2.regenerate()
a = mdb.models['Alfano'].rootAssembly
v1 = a.instances['Plate-1'].vertices
verts1 = v1[3:4]
a.Set(vertices=verts1, name='topBC')
a = mdb.models['Alfano'].rootAssembly
v1 = a.instances['Plate-1'].vertices
verts1 = v1[6:7]
a.Set(vertices=verts1, name='botBC')
session.viewports['Viewport: 1'].view.fitView()
session.viewports['Viewport: 1'].assemblyDisplay.setValues(mesh=OFF)
session.viewports['Viewport: 1'].assemblyDisplay.meshOptions.setValues(
    meshTechnique=OFF)
#
# Create *STATIC step
#
mdb.models['Alfano'].StaticStep(name='Peel', previous='Initial', 
    description='Pulling apart the plates', maxNumInc=5000, 
    stabilizationMethod=DAMPING_FACTOR, stabilizationMagnitude=700000.0, 
    initialInc=0.01, minInc=1e-08, nlgeom=ON)
session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Peel')
mdb.models['Alfano'].fieldOutputRequests['F-Output-1'].setValues(
    frequency=9999)
a = mdb.models['Alfano'].rootAssembly
v1 = a.instances['Plate-1'].vertices
verts1 = v1[3:4]+v1[6:7]
a.Set(vertices=verts1, name='dispBC')
regionDef=mdb.models['Alfano'].rootAssembly.sets['dispBC']
mdb.models['Alfano'].HistoryOutputRequest(name='H-Output-2', 
    createStepName='Peel', variables=('U2', 'RF2'), region=regionDef, 
    sectionPoints=DEFAULT, rebar=EXCLUDE)
mdb.models['Alfano'].fieldOutputRequests['F-Output-1'].setValues(variables=(
    'S', 'U'))
a = mdb.models['Alfano'].rootAssembly
v1 = a.instances['Plate-1'].vertices
verts1 = v1[3:4]
region=regionToolset.Region(vertices=verts1)
mdb.models['Alfano'].steps['Peel'].Monitor(dof=2, node=region, frequency=1)
mdb.models['Alfano'].steps['Peel'].control.setValues(allowPropagation=OFF, 
    resetDefaultValues=OFF, discontinuous=ON, lineSearch=(4.0, 4.0, 0.25, 0.25, 
    0.15))
session.viewports['Viewport: 1'].assemblyDisplay.setValues(loads=ON, bcs=ON, 
    fields=ON, connectors=ON)
session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Initial')
a = mdb.models['Alfano'].rootAssembly
region = a.sets['dispBC']
mdb.models['Alfano'].DisplacementBC(name='fix', createStepName='Initial', 
    region=region, u1=SET, u2=UNSET, ur3=UNSET, amplitude=UNSET, 
    distribution=UNIFORM, localCsys=None)
session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Peel')
a = mdb.models['Alfano'].rootAssembly
region = a.sets['topBC']
mdb.models['Alfano'].DisplacementBC(name='applyDisp', createStepName='Peel', 
    region=region, u1=UNSET, u2=0.02, ur3=UNSET, amplitude=UNSET, fixed=OFF, 
    distribution=UNIFORM, localCsys=None)
mdb.models['Alfano'].boundaryConditions.changeKey(fromName='applyDisp', 
    toName='topDisp')
a = mdb.models['Alfano'].rootAssembly
region = a.sets['botBC']
mdb.models['Alfano'].DisplacementBC(name='botDisp', createStepName='Peel', 
    region=region, u1=UNSET, u2=-0.02, ur3=UNSET, amplitude=UNSET, fixed=OFF, 
    distribution=UNIFORM, localCsys=None)
session.viewports['Viewport: 1'].assemblyDisplay.setValues(loads=OFF, bcs=OFF, 
    fields=OFF, interactions=ON, constraints=ON, engineeringFeatures=ON)
a = mdb.models['Alfano'].rootAssembly
region1=a.surfaces['top']
a = mdb.models['Alfano'].rootAssembly
region2=a.surfaces['topCoh-Top']
mdb.models['Alfano'].Tie(name='cohTop-Top', master=region1, slave=region2, 
    positionToleranceMethod=COMPUTED, adjust=ON, tieRotations=ON)
a = mdb.models['Alfano'].rootAssembly
region1=a.surfaces['middle-top']
a = mdb.models['Alfano'].rootAssembly
region2=a.surfaces['topCoh-Bot']
mdb.models['Alfano'].Tie(name='cohTop-Middle', master=region1, slave=region2, 
    positionToleranceMethod=COMPUTED, adjust=ON, tieRotations=ON)
a = mdb.models['Alfano'].rootAssembly
region1=a.surfaces['middle-bot']
a = mdb.models['Alfano'].rootAssembly
region2=a.surfaces['botCoh-Top']
mdb.models['Alfano'].Tie(name='botCoh-Middle', master=region1, slave=region2, 
    positionToleranceMethod=COMPUTED, adjust=ON, tieRotations=ON)
a = mdb.models['Alfano'].rootAssembly
region1=a.surfaces['bottom']
a = mdb.models['Alfano'].rootAssembly
region2=a.surfaces['botCoh-Bot']
mdb.models['Alfano'].Tie(name='botCoh-Bottom', master=region1, slave=region2, 
    positionToleranceMethod=COMPUTED, adjust=ON, tieRotations=ON)
mdb.models['Alfano'].ContactProperty('int')
mdb.models['Alfano'].interactionProperties['int'].TangentialBehavior(
    formulation=FRICTIONLESS)
#: The interaction property "int" has been created.
session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Initial')
a = mdb.models['Alfano'].rootAssembly
region1=a.surfaces['middle-bot']
a = mdb.models['Alfano'].rootAssembly
region2=a.surfaces['bottom']
mdb.models['Alfano'].SurfaceToSurfaceContactStd(name='lower', 
    createStepName='Initial', master=region1, slave=region2, sliding=FINITE, 
    interactionProperty='int', adjustMethod=NONE)
#: The interaction "lower" has been created.
a = mdb.models['Alfano'].rootAssembly
region1=a.surfaces['middle-top']
a = mdb.models['Alfano'].rootAssembly
region2=a.surfaces['top']
mdb.models['Alfano'].SurfaceToSurfaceContactStd(name='upper', 
    createStepName='Initial', master=region1, slave=region2, sliding=FINITE, 
    interactionProperty='int', adjustMethod=NONE)
#: The interaction "upper" has been created.
session.viewports['Viewport: 1'].assemblyDisplay.setValues(interactions=OFF, 
    constraints=OFF, connectors=OFF, engineeringFeatures=OFF)
mdb.Job(name='alfano_2d_tie_std', model='Alfano', type=ANALYSIS, 
    explicitPrecision=SINGLE, nodalOutputPrecision=SINGLE, 
    description='Alfano Crisfield non-symmetric multi-delamination analysis', 
    parallelizationMethodExplicit=DOMAIN, multiprocessingMode=DEFAULT, 
    numDomains=1, userSubroutine='', numCpus=1, preMemory=1024.0, 
    standardMemory=1024.0, standardMemoryPolicy=MODERATE, scratch='', 
    echoPrint=OFF, modelPrint=OFF, contactPrint=OFF, historyPrint=OFF)
session.viewports['Viewport: 1'].view.fitView()
mdb.jobs['alfano_2d_tie_std'].writeInput()
#: The job input file has been written to "alfano_2d_tie_std.inp".
mdb.saveAs(pathName='alfano')