#
# Getting Started with ABAQUS
#
# Script for dynamic crane example
#
import string
def GetBlockPosition(modelName, blockPrefix):
if blockPrefix == '':
return len(mdb.models[modelName].keywordBlock.sieBlocks)-1
pos = 0
for block in mdb.models[modelName].keywordBlock.sieBlocks:
if string.lower(block[0:len(blockPrefix)])==string.lower(blockPrefix):
return pos
pos=pos+1
return -1
from abaqus import *
from abaqusConstants import *
session.viewports['Viewport: 1'].makeCurrent()
session.viewports['Viewport: 1'].maximize()
from caeModules import *
from driverUtils import executeOnCaeStartup
executeOnCaeStartup()
Mdb()
mdb.models.changeKey(fromName='Model-1', toName='standard')
##
## Sketch profile of truss
##
p = mdb.models['standard'].Part(name='Truss', dimensionality=THREE_D,
type=DEFORMABLE_BODY)
p.ReferencePoint(point=(0.0, 0.0, 0.0))
p = mdb.models['standard'].parts['Truss']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
d, r = p.datums, p.referencePoints
p.DatumPointByOffset(point=r[1], vector=(0.0, 1.0, 0.0))
p.DatumPointByOffset(point=r[1], vector=(8.0, 1.5, 0.9))
p.DatumPlaneByThreePoints(point1=r[1], point2=d[3], point3=d[2])
p.DatumAxisByPrincipalAxis(principalAxis=YAXIS)
session.viewports['Viewport: 1'].view.fitView()
t = p.MakeSketchTransform(sketchPlane=d[4], sketchUpEdge=d[5],
sketchPlaneSide=SIDE1, sketchOrientation=LEFT, origin=(0.0, 0.0, 0.0))
s = mdb.models['standard'].ConstrainedSketch(name='__profile__',
sheetSize=683.38, gridSpacing=17.08, transform=t)
g, v, d1, c = s.geometry, s.vertices, s.dimensions, s.constraints
s.setPrimaryObject(option=SUPERIMPOSE)
p = mdb.models['standard'].parts['Truss']
p.projectReferencesOntoSketch(sketch=s, filter=COPLANAR_EDGES)
session.viewports['Viewport: 1'].view.setValues(nearPlane=648.075,
farPlane=719.271, width=13.6266, height=9.28111, cameraPosition=(-71.5628,
0.978423, 679.645), cameraTarget=(4.83613, 0.978423, 0.544065))
s.Line(point1=(0.0, 0.0), point2=(8.05046582503149, 1.5))
s.CoincidentConstraint(entity1=v.findAt((0.0, 0.0)), entity2=g.findAt((0.0,
0.5)))
s.Line(point1=(8.05046582503149, 1.5), point2=(0.0, 1.0))
s.Line(point1=(0.0, 1.0), point2=(1.59409310185916, 0.297018794285658))
s.CoincidentConstraint(entity1=v.findAt((1.594093, 0.297019)),
entity2=g.findAt((4.025233, 0.75)))
s.Line(point1=(1.59409310185916, 0.297018794285658),
point2=(2.30592731982451, 1.14321701190598))
s.CoincidentConstraint(entity1=v.findAt((2.305927, 1.143217)),
entity2=g.findAt((4.025233, 1.25)))
s.Line(point1=(2.30592731982451, 1.14321701190598),
point2=(3.26033977276254, 0.607481575033539))
s.CoincidentConstraint(entity1=v.findAt((3.26034, 0.607482)),
entity2=g.findAt((4.025233, 0.75)))
s.Line(point1=(3.26033977276254, 0.607481575033539),
point2=(4.16525059699136, 1.25869624736747))
s.CoincidentConstraint(entity1=v.findAt((4.165251, 1.258696)),
entity2=g.findAt((4.025233, 1.25)))
s.Line(point1=(4.16525059699136, 1.25869624736747),
point2=(4.95150583156061, 0.922587451305882))
s.CoincidentConstraint(entity1=v.findAt((4.951506, 0.922587)),
entity2=g.findAt((4.025233, 0.75)))
s.breakCurve(curve1=g.findAt((4.025233, 0.75)),
point1=(0.733189206279801, 0.182644009590149),
curve2=g.findAt((1.95001, 0.720118)),
point2=(1.62787567272877, 0.433001518249512))
s.breakCurve(curve1=g.findAt((4.822279, 0.898509)),
point1=(2.34362044857835, 0.415118813514709),
curve2=g.findAt((2.783134, 0.875349)),
point2=(3.22040910681461, 0.701241612434387))
s.breakCurve(curve1=g.findAt((5.655403, 1.053741)),
point1=(4.04351950354929, 0.772772312164307),
curve2=g.findAt((4.558378, 1.090642)),
point2=(5.02767302451382, 0.95159912109375))
s.breakCurve(curve1=g.findAt((4.025233, 1.25)),
point1=(1.55629884485765, 1.09466052055359),
curve2=g.findAt((1.95001, 0.720118)),
point2=(2.32572264036564, 1.02312982082367))
s.breakCurve(curve1=g.findAt((5.178197, 1.321609)),
point1=(3.79300780849022, 1.23772192001343),
curve2=g.findAt((3.712795, 0.933089)),
point2=(4.168774592922, 1.18407392501831))
s.VerticalDimension(vertex1=v.findAt((0.0, 1.0)),
vertex2=v.findAt((0.0, 0.0)),
textPoint=(-0.304643334343144, 0.0216999053955078), value=1.0)
s.HorizontalDimension(vertex1=v.findAt((0.0, 0.0)),
vertex2=v.findAt((8.050466, 1.5)),
textPoint=(6.26233672422377, -0.443249702453613), value=8.05046582503149)
s.EqualLengthConstraint(entity1=g.findAt((0.797047, 0.148509)),
entity2=g.findAt((2.427216, 0.45225)))
s.EqualLengthConstraint(entity1=g.findAt((2.39114, 0.445528)),
entity2=g.findAt((4.033769, 0.751591)))
s.EqualLengthConstraint(entity1=g.findAt((3.985233, 0.742547)),
entity2=g.findAt((6.416373, 1.195528)))
s.EqualLengthConstraint(entity1=g.findAt((1.152964, 1.071609)),
entity2=g.findAt((3.235589, 1.200957)))
s.EqualLengthConstraint(entity1=g.findAt((3.458891, 1.214826)),
entity2=g.findAt((6.33116, 1.393217)))
mdb.models['standard'].ConstrainedSketch(name='Truss', objectToCopy=s)
p = mdb.models['standard'].parts['Truss']
d = p.datums
p.Wire(sketchPlane=d[4], sketchUpEdge=d[5], sketchPlaneSide=SIDE1,
sketchOrientation=LEFT, sketch=s)
s.unsetPrimaryObject()
del mdb.models['standard'].sketches['__profile__']
session.viewports['Viewport: 1'].view.fitView()
p = mdb.models['standard'].parts['Truss']
d, r = p.datums, p.referencePoints
p.DatumPointByOffset(point=d[2], vector=(0.0, 0.0, 2.0))
p.DatumPointByOffset(point=r[1], vector=(0.0, 0.0, 2.0))
p.DatumPointByOffset(point=d[3], vector=(0.0, 0.0, 0.2))
p.DatumPlaneByThreePoints(point1=d[8], point2=d[9], point3=d[7])
t = p.MakeSketchTransform(sketchPlane=d[10], sketchUpEdge=d[5],
sketchPlaneSide=SIDE1, sketchOrientation=LEFT, origin=(4.122203, 0.75,
1.536252))
s = mdb.models['standard'].ConstrainedSketch(name='__profile__',
sheetSize=683.38, gridSpacing=17.08, transform=t)
g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraints
s.setPrimaryObject(option=SUPERIMPOSE)
p.projectReferencesOntoSketch(sketch=s, filter=COPLANAR_EDGES)
s.retrieveSketch(sketch=mdb.models['standard'].sketches['Truss'])
s.move(vector=(-4.14820681508444, -0.75), objectList=(g.findAt((0.0, 0.5)),
g.findAt((1.0063, 0.687498)), g.findAt((2.3482, 0.770841)), g.findAt((
3.3545, 0.95834)), g.findAt((4.6959, 1.041652)), g.findAt((5.7022,
1.22915)), g.findAt((1.0063, 0.187498)), g.findAt((3.0189, 0.562495)),
g.findAt((5.0315, 0.937492)), g.findAt((7.044133, 1.312495)), g.findAt((
1.3419, 1.083343)), g.findAt((6.708533, 1.416655)), g.findAt((4.0252,
1.249998))))
d = p.datums
p.Wire(sketchPlane=d[10], sketchUpEdge=d[5], sketchPlaneSide=SIDE1,
sketchOrientation=LEFT, sketch=s)
s.unsetPrimaryObject()
del mdb.models['standard'].sketches['__profile__']
##
## dummy part
##
p = mdb.models['standard'].Part(name='Truss-all',
objectToCopy=mdb.models['standard'].parts['Truss'])
session.viewports['Viewport: 1'].setValues(displayedObject=p)
v = p.vertices
p.WirePolyLine(points=((v.findAt(coordinates=(2.0, 0.375, 1.775)), v.findAt(
coordinates=(2.0, 0.375, 0.225))), (v.findAt(coordinates=(2.0, 0.375,
0.225)), v.findAt(coordinates=(2.666667, 1.166667, 1.7))), (v.findAt(
coordinates=(2.666667, 1.166667, 1.7)), v.findAt(coordinates=(2.666667,
1.166667, 0.3))), (v.findAt(coordinates=(4.0, 0.75, 1.55)), v.findAt(
coordinates=(4.0, 0.75, 0.45))), (v.findAt(coordinates=(4.0, 0.75, 0.45)),
v.findAt(coordinates=(5.333333, 1.333333, 1.4))), (v.findAt(coordinates=(
5.333333, 1.333333, 1.4)), v.findAt(coordinates=(5.333333, 1.333333,
0.6))), (v.findAt(coordinates=(5.333333, 1.333333, 0.6)), v.findAt(
coordinates=(6.0, 1.125, 1.325))), (v.findAt(coordinates=(6.0, 1.125,
1.325)), v.findAt(coordinates=(6.0, 1.125, 0.675)))), mergeWire=OFF,
meshable=ON)
##
## create cross brace part in the assembly module
##
a = mdb.models['standard'].rootAssembly
session.viewports['Viewport: 1'].setValues(displayedObject=a)
## Set coordinate system (done by default)
a.DatumCsysByDefault(CARTESIAN)
p = mdb.models['standard'].parts['Truss']
a.Instance(name='Truss-1', part=p, dependent=ON)
p = mdb.models['standard'].parts['Truss-all']
a.Instance(name='Truss-all-1', part=p, dependent=ON)
a.PartFromBooleanCut(name='Cross brace',
instanceToBeCut=a.instances['Truss-all-1'],
cuttingInstances=(a.instances['Truss-1'], ))
p = mdb.models['standard'].parts['Cross brace']
a.Instance(name='Cross brace-1', part=p, dependent=ON)
a.suppressFeatures(('Truss-1', 'Truss-all-1', ))
a.features['Truss-1'].resume()
p = mdb.models['standard'].parts['Cross brace']
session.viewports['Viewport: 1'].setValues(displayedObject=p)
##
## Create beam profiles and beam sections
##
mdb.models['standard'].BoxProfile(name='MainBoxProfile', b=0.05,
a=0.1, uniformThickness=ON, t1=0.005)
mdb.models['standard'].BoxProfile(name='BraceBoxProfile', b=0.03,
a=0.03, uniformThickness=ON, t1=0.003)
mdb.models['standard'].BeamSection(name='MainMemberSection',
profile='MainBoxProfile',
integration=BEFORE_ANALYSIS, table=((2.e11, 8.e10),))
mdb.models['standard'].sections['MainMemberSection'].setValues(density=7800.0)
mdb.models['standard'].BeamSection(name='BracingSection',
profile='BraceBoxProfile',
integration=BEFORE_ANALYSIS, table=((2.e11, 8.e10),))
mdb.models['standard'].sections['BracingSection'].setValues(density=7800.0)
##
## Flip tangent directions
##
p = mdb.models['standard'].parts['Truss']
e = p.edges
edges = e
regions = regionToolset.Region(edges=edges)
p.flipTangent(regions=regions)
##
## Assign beam sections to the cross brace
##
p = mdb.models['standard'].parts['Cross brace']
e = p.edges
edges = e.findAt(((2.666667, 1.166667, 1.35), ),
((2.166667, 0.572917, 0.59375), ), ((2.0, 0.375, 1.3875), ),
((6.0, 1.125, 1.1625), ), ((5.5, 1.28125, 0.78125), ),
((5.333333, 1.333333, 1.2), ), ((4.333333, 0.895833, 0.6875), ),
((4.0, 0.75, 1.275), ))
region = regionToolset.Region(edges=edges)
p.SectionAssignment(region=region, sectionName='BracingSection')
##
## Assign beam sections to the internal bracing
##
p = mdb.models['standard'].parts['Truss']
e = p.edges
edges = e.findAt(((3.0, 1.0625, 0.3375), ), ((5.5, 1.28125, 0.61875), ),
((4.333333, 0.895833, 0.4875), ), ((0.5, 0.84375, 0.05625), ),
((2.166667, 0.572917, 0.24375), ), ((3.0, 1.0625, 1.6625), ),
((5.5, 1.28125, 1.38125), ), ((4.333333, 0.895833, 1.5125), ),
((0.5, 0.84375, 1.94375), ), ((2.166667, 0.572917, 1.75625), ))
region = regionToolset.Region(edges=edges)
p.SectionAssignment(region=region, sectionName='BracingSection')
##
## Assign beam sections to the main members
##
p = mdb.models['standard'].parts['Truss']
e = p.edges
edges = e.findAt(((4.5, 0.84375, 0.50625), ), ((4.666667, 1.291667, 0.525), ),
((6.5, 1.21875, 0.73125), ), ((7.333333, 1.458333, 0.825), ),
((2.5, 0.46875, 0.28125), ), ((2.0, 1.125, 0.225), ),
((0.5, 0.09375, 0.05625), ), ((4.5, 0.84375, 1.49375), ),
((4.666667, 1.291667, 1.475), ), ((6.5, 1.21875, 1.26875), ),
((7.333333, 1.458333, 1.175), ), ((2.5, 0.46875, 1.71875), ),
((2.0, 1.125, 1.775), ), ((0.5, 0.09375, 1.94375), ))
region = regionToolset.Region(edges=edges)
p.SectionAssignment(region=region, sectionName='MainMemberSection')
##
## Assign beam section orientations to truss 'B'
##
p = mdb.models['standard'].parts['Truss']
e = p.edges
edges = e.findAt(((4.5, 0.84375, 0.50625), ), ((3.0, 1.0625, 0.3375), ),
((4.666667, 1.291667, 0.525), ), ((5.5, 1.28125, 0.61875), ),
((6.5, 1.21875, 0.73125), ), ((7.333333, 1.458333, 0.825), ),
((4.333333, 0.895833, 0.4875), ), ((2.5, 0.46875, 0.28125), ),
((0.5, 0.84375, 0.05625), ), ((2.0, 1.125, 0.225), ),
((2.166667, 0.572917, 0.24375), ), ((0.5, 0.09375, 0.05625), ))
region=regionToolset.Region(edges=edges)
p.assignBeamSectionOrientation(region=region, method=N1_COSINES, n1=(-0.1118,
0.0, 0.9936))
##
## Assign beam section orientations to truss 'A'
##
p = mdb.models['standard'].parts['Truss']
e = p.edges
edges = e.findAt(((4.5, 0.84375, 1.49375), ), ((3.0, 1.0625, 1.6625), ),
((4.666667, 1.291667, 1.475), ), ((5.5, 1.28125, 1.38125), ),
((6.5, 1.21875, 1.26875), ), ((7.333333, 1.458333, 1.175), ),
((4.333333, 0.895833, 1.5125), ), ((2.5, 0.46875, 1.71875), ),
((0.5, 0.84375, 1.94375), ), ((2.0, 1.125, 1.775), ),
((2.166667, 0.572917, 1.75625), ), ((0.5, 0.09375, 1.94375), ))
region=regionToolset.Region(edges=edges)
p.assignBeamSectionOrientation(region=region, method=N1_COSINES, n1=(-0.1118,
0.0, -0.9936))
##
## Assign beam section orientations to 'Cross brace'
##
p = mdb.models['standard'].parts['Cross brace']
e = p.edges
edges = e.findAt(((2.666667, 1.166667, 1.35), ), ((2.166667, 0.572917,
0.59375), ), ((2.0, 0.375, 1.3875), ), ((6.0, 1.125, 1.1625), ), ((5.5,
1.28125, 0.78125), ), ((5.333333, 1.333333, 1.2), ), ((4.333333, 0.895833,
0.6875), ), ((4.0, 0.75, 1.275), ))
region=regionToolset.Region(edges=edges)
p.assignBeamSectionOrientation(region=region, method=N1_COSINES,
n1=(0.0, 1.0, 0.0))
session.viewports['Viewport: 1'].setValues(displayedObject=a)
a.regenerate()
##
##
session.viewports['Viewport: 1'].assemblyDisplay.setValues(
visibleInstances=('Truss-1', 'Cross brace-1'))
session.viewports['Viewport: 1'].assemblyDisplay.geometryOptions.setValues(
geometryEdgesInShaded=OFF, datumPoints=OFF, datumAxes=OFF, datumPlanes=OFF,
datumCoordSystems=OFF, referencePointLabels=OFF, referencePointSymbols=OFF)
##
## Translate the truss
##
p = a.instances['Truss-1']
p.translate(vector=(0.0, -0.5, -1.0))
##
## Translate the cross brace
##
p = a.instances['Cross brace-1']
p.translate(vector=(0.0, -0.5, -1.0))
##
## Create geometry sets
##
session.viewports['Viewport: 1'].assemblyDisplay.setValues(
visibleInstances=('Truss-1', ))
session.viewports['Viewport: 1'].assemblyDisplay.geometryOptions.setValues(
datumPoints=OFF)
v = a.instances['Truss-1'].vertices
verts = v.findAt(((0.0, 0.5, -1.0), ), ((0.0, -0.5, -1.0), ), ((0.0, 0.5,
1.0), ), ((0.0, -0.5, 1.0), ))
a.Set(vertices=verts, name='Attach')
verts = v.findAt(((8.0, 1.0, 0.1), ))
a.Set(vertices=verts, name='Tip-a')
verts = v.findAt(((8.0, 1.0, -0.1), ))
a.Set(vertices=verts, name='Tip-b')
e = a.instances['Truss-1'].edges
edges = e.findAt(((6.5, 0.71875, 0.26875), ))
a.Set(edges=edges, name='Leg-a')
edges = e.findAt(((6.5, 0.71875, -0.26875), ))
a.Set(edges=edges, name='Leg-b')
edges = e.findAt(((4.333333, 0.395833, 0.5125), ))
a.Set(edges=edges, name='Inner-a')
edges = e.findAt(((4.333333, 0.395833, -0.5125), ))
a.Set(edges=edges, name='Inner-b')
session.viewports['Viewport: 1'].assemblyDisplay.setValues(
visibleInstances=('Cross brace-1', ))
session.viewports['Viewport: 1'].view.setValues(
cameraPosition=(7.4736, 8.3559, 20.002),
cameraUpVector=(-0.34023, 0.72044, -0.60432))
session.viewports['Viewport: 1'].assemblyDisplay.setValues(
visibleInstances=('Truss-1', 'Cross brace-1'))
session.viewports['Viewport: 1'].view.setValues(session.views['Iso'])
##
## Create a frequency extraction step
##
mdb.models['standard'].FrequencyStep(name='Extract frequencies',
previous='Initial',
description='First 30 modes', numEigen=30,
maxEigen=None, shift=None, eigensolver=LANCZOS, minEigen=None,
blockSize=DEFAULT, maxBlocks=DEFAULT)
session.viewports['Viewport: 1'].assemblyDisplay.setValues(
step='Extract frequencies')
##
## Create a modal dynamics step
##
mdb.models['standard'].ModalDynamicsStep(name='Transient modal dynamics',
previous='Extract frequencies',
description='Simulation of Load Dropped on Crane',
continueAnalysis=OFF, timePeriod=0.5, incSize=0.005, directDamping=((1,
30, 0.05), ), compositeDamping=None, rayleighDamping=None, amplitude=STEP)
session.viewports['Viewport: 1'].assemblyDisplay.setValues(
step='Transient modal dynamics')
##
## Modify output requests
##
mdb.models['standard'].fieldOutputRequests['F-Output-2'].setValues(
variables=('U', ), frequency=5)
mdb.models['standard'].historyOutputRequests['H-Output-1'].setValues(
variables=('ALLIE', 'ALLKE', 'ALLVD'), frequency=1)
a0=mdb.models['standard'].rootAssembly
regionDef=a0.sets['Tip-a']
mdb.models['standard'].HistoryOutputRequest(name='H-Output-2',
createStepName='Transient modal dynamics',
variables=('U1', 'U2', 'U3'),
region=regionDef, frequency=1)
a0=mdb.models['standard'].rootAssembly
regionDef=a0.sets['Attach']
mdb.models['standard'].HistoryOutputRequest(name='H-Output-3',
createStepName='Transient modal dynamics',
variables=('RF1', 'RF2', 'RF3'),
region=regionDef, frequency=1)
##
## Create constraints between the tips of the truss
##
session.viewports['Viewport: 1'].assemblyDisplay.setValues(interactions=ON)
mdb.models['standard'].Equation(name='TipConstraint-1',
terms=((1.0, 'Tip-a', 1), (-1.0, 'Tip-b', 1)))
mdb.models['standard'].Constraint('TipConstraint-2',
mdb.models['standard'].constraints['TipConstraint-1'])
mdb.models['standard'].constraints['TipConstraint-2'].setValues(
terms=((1.0, 'Tip-a', 2), (-1.0, 'Tip-b', 2)))
##
## Create JOIN connectors between truss and bracing
##
v1 = a.instances['Cross brace-1'].vertices
v2 = a.instances['Truss-1'].vertices
a.WirePolyLine(points=((v1.findAt(coordinates=(2.0, -0.125, 0.775)), v2.findAt(
coordinates=(2.0, -0.125, 0.775))), (v1.findAt(coordinates=(2.666667,
0.666667, 0.7)), v2.findAt(coordinates=(2.666667, 0.666667, 0.7))), (
v1.findAt(coordinates=(2.0, -0.125, -0.775)), v2.findAt(coordinates=(2.0,
-0.125, -0.775))), (v1.findAt(coordinates=(2.666667, 0.666667, -0.7)),
v2.findAt(coordinates=(2.666667, 0.666667, -0.7))), (v1.findAt(
coordinates=(4.0, 0.25, -0.55)), v2.findAt(coordinates=(4.0, 0.25,
-0.55))), (v1.findAt(coordinates=(4.0, 0.25, 0.55)), v2.findAt(
coordinates=(4.0, 0.25, 0.55))), (v1.findAt(coordinates=(5.333333,
0.833333, 0.4)), v2.findAt(coordinates=(5.333333, 0.833333, 0.4))), (
v1.findAt(coordinates=(5.333333, 0.833333, -0.4)), v2.findAt(coordinates=(
5.333333, 0.833333, -0.4))), (v1.findAt(coordinates=(6.0, 0.625, -0.325)),
v2.findAt(coordinates=(6.0, 0.625, -0.325))), (v1.findAt(coordinates=(6.0,
0.625, 0.325)), v2.findAt(coordinates=(6.0, 0.625, 0.325)))),
mergeWire=OFF, meshable=OFF)
e1 = a.edges
edges1 = e1.findAt(((6.000025, 0.625, 0.325), ), ((6.000025, 0.625, -0.325), ),
((5.333358, 0.833333, -0.4), ), ((5.333358, 0.833333, 0.4), ), ((4.000025,
0.25, 0.55), ), ((4.000025, 0.25, -0.55), ), ((2.666692, 0.666667, -0.7),
), ((2.000025, -0.125, -0.775), ), ((2.666692, 0.666667, 0.7), ), ((
2.000025, -0.125, 0.775), ))
a.Set(edges=edges1, name='Wire-1-Edge-1')
mdb.models['standard'].ConnectorSection(name='ConnSect-1',
translationalType=JOIN)
region=a.sets['Wire-1-Edge-1']
a.SectionAssignment(sectionName='ConnSect-1', region=region)
session.viewports['Viewport: 1'].assemblyDisplay.setValues(loads=ON, bcs=ON,
predefinedFields=ON, interactions=OFF)
##
## Apply encastre bc to set 'Attach'
##
session.viewports['Viewport: 1'].assemblyDisplay.setValues(
step='Extract frequencies')
a = mdb.models['standard'].rootAssembly
region = a.sets['Attach']
mdb.models['standard'].EncastreBC(name='Fixed end',
createStepName='Extract frequencies', region=region)
##
## Create amplitude curve 'Bounce'
##
session.viewports['Viewport: 1'].assemblyDisplay.setValues(
step='Transient modal dynamics')
mdb.models['standard'].TabularAmplitude(name='Bounce', timeSpan=STEP,
smooth=0.25, data=((0.0, 0.0), (0.01, 1.0), (0.2, 1.0), (0.21, 0.0)))
##
## Apply concentrated force to set 'Tip-b'
##
a = mdb.models['standard'].rootAssembly
region = a.sets['Tip-b']
mdb.models['standard'].ConcentratedForce(name='Tip load',
createStepName='Transient modal dynamics', region=region, cf2=-10000.0,
amplitude='Bounce')
##
## Assign global seed
##
p = mdb.models['standard'].parts['Truss']
p.seedPart(size=2.0)
p = mdb.models['standard'].parts['Cross brace']
p.seedPart(size=2.0)
##
## Assign element type
##
elemType1 = mesh.ElemType(elemCode=B33)
p = mdb.models['standard'].parts['Cross brace']
e = p.edges
edges = e
pickedRegions =(edges, )
p.setElementType(regions=pickedRegions, elemTypes=(elemType1, ))
p = mdb.models['standard'].parts['Truss']
e = p.edges
edges = e
pickedRegions =(edges, )
p.setElementType(regions=pickedRegions, elemTypes=(elemType1, ))
##
## Generate mesh
##
p = mdb.models['standard'].parts['Truss']
p.generateMesh()
p = mdb.models['standard'].parts['Cross brace']
p.generateMesh()
##
## Add keywords
##
mdb.models['standard'].keywordBlock.synchVersions()
mdb.models['standard'].keywordBlock.insert(GetBlockPosition('standard', '*End Assembly')-1, """*NORMAL, TYPE=ELEMENT
Inner-a, Inner-a, -0.3986, 0.9114, 0.1025
Inner-b, Inner-b, 0.3986, -0.9114, 0.1025
Leg-a, Leg-a, -0.1820, 0.9829, 0.0205
Leg-b, Leg-b, 0.1820, -0.9829, 0.0205""")
##
## Create job
##
mdb.Job(name='DynCrane', model='standard',
description='3-D model of light-service cargo crane')
a = mdb.models['standard'].rootAssembly
a.regenerate()
##
## Save model database
##
mdb.saveAs('DynCrane.cae')
mdb.Model(name='explicit', objectToCopy=mdb.models['standard'])
a = mdb.models['explicit'].rootAssembly
session.viewports['Viewport: 1'].setValues(displayedObject=a)
##
## delete modal dynamics step
##
del mdb.models['explicit'].steps['Transient modal dynamics']
##
## replace frequency extraction step with explicit dynamics step
##
mdb.models['explicit'].ExplicitDynamicsStep(name='Extract frequencies',
previous='Initial', maintainAttributes=TRUE, timePeriod=0.5, nlgeom=OFF)
##
## rename explicit dynamics step
##
mdb.models['explicit'].steps.changeKey(fromName='Extract frequencies',
toName='Transient dynamics')
##
## Redefine history output
##
regionDef=mdb.models['explicit'].rootAssembly.sets['Tip-a']
mdb.models['explicit'].HistoryOutputRequest(name='H-Output-2',
createStepName='Transient dynamics', variables=('U1', 'U2', 'U3'),
region=regionDef, sectionPoints=DEFAULT, rebar=EXCLUDE)
regionDef=mdb.models['explicit'].rootAssembly.sets['Attach']
mdb.models['explicit'].HistoryOutputRequest(name='H-Output-3',
createStepName='Transient dynamics', variables=('RF1', 'RF2', 'RF3'),
region=regionDef, sectionPoints=DEFAULT, rebar=EXCLUDE)
##
## Add mass proportional damping
##
mdb.models['explicit'].sections['BracingSection'].setValues(poissonRatio=0.0,
density=7800.0, thermalExpansion=OFF, temperatureDependency=OFF,
dependencies=0, table=((2.e11, 8.e10), ),
alphaDamping=15.0, betaDamping=0.0, compositeDamping=0.0, centroid=(0.0,
0.0), shearCenter=(0.0, 0.0))
mdb.models['explicit'].sections['MainMemberSection'].setValues(
poissonRatio=0.0, density=7800.0, thermalExpansion=OFF,
temperatureDependency=OFF, dependencies=0, table=((2.e11,
8.e10), ), alphaDamping=15.0, betaDamping=0.0,
compositeDamping=0.0, centroid=(0.0, 0.0), shearCenter=(0.0, 0.0))
##
## Redefine tip load
##
a = mdb.models['explicit'].rootAssembly
region = a.sets['Tip-b']
mdb.models['explicit'].ConcentratedForce(name='Load-1',
createStepName='Transient dynamics', region=region, cf2=-10000.0,
amplitude='Bounce', localCsys=None)
##
## Change element library
##
elemType1 = mesh.ElemType(elemCode=B31, elemLibrary=EXPLICIT)
p = mdb.models['explicit'].parts['Cross brace']
e = p.edges
edges = e
pickedRegions =(edges, )
p.setElementType(regions=pickedRegions, elemTypes=(elemType1, ))
p = mdb.models['explicit'].parts['Truss']
elemType1 = mesh.ElemType(elemCode=B31, elemLibrary=EXPLICIT)
e = p.edges
edges = e
pickedRegions =(edges, )
p.setElementType(regions=pickedRegions, elemTypes=(elemType1, ))
a = mdb.models['explicit'].rootAssembly
a.regenerate()
##
## Create explicit job
##
mdb.Job(name='expDynCrane', model='explicit', type=ANALYSIS)
a = mdb.models['standard'].rootAssembly
a.regenerate()
a = mdb.models['explicit'].rootAssembly
a.regenerate()
session.viewports['Viewport: 1'].setValues(displayedObject=a)
mdb.save()
