1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
#!/user/bin/python
# -* - coding:UTF-8 -*-
#PRE
from abaqus import *
from abaqusConstants import *
import regionToolset
session.viewports['Viewport: 1'].setValues(displayedObject=None)
mdb.models.changeKey(fromName='Model-1',toName='Cantilever Beam')
beamModel = mdb.models['Cantilever Beam']
import sketch
import part
beamProfileSketch = beamModel.ConstrainedSketch(name='Beam CS Profile',
sheetSize=5)
beamProfileSketch.rectangle(point1=(0.1,0.1),point2=(0.3,-0.1))
beamPart = beamModel.Part(name='Beam',dimensionality=THREE_D,
type=DEFORMABLE_BODY)
beamPart.BaseSolidExtrude(sketch=beamProfileSketch,depth=5)
import material
beamMaterial = beamModel
.Material(name='AISI 1005 Steel')
beamMaterial.Density(table=((7872, ), ))
beamMaterial.Elastic(table=((200E9,0.29), ))
import section
beamSection = beamModel.HomogeneousSolidSection(name='Beam Section',
material='AISI 1005 Steel')
beam_region = regionToolset.Region(cells=beamPart.cells)
beamPart.SectionAssignment(region=beam_region,sectionName='Beam Section')
import assembly
beamAssembly = beamModel.rootAssembly
beamInstance = beamAssembly.Instance(name='Beam Instance',part=beamPart,
dependent=ON)
import step
beamModel.StaticStep(name='Apply Load',previous='Initial',
description='Load is applied during this step')
beamModel.fieldOutputRequests.changeKey(fromName='F-Output-1',
toName='Selected Field Outputs')
beamModel
.fieldOutputRequests['Selected Field Outputs'].setValues(variables=('S',
'E','PEMAG','U','RF','CF'))
beamModel.HistoryOutputRequest(name='Default History Outputs',
createStepName='Apply Load',variables=PRESELECT)
del beamModel.historyOutputRequests['H-Output-1']
top_face_pt_x = 0.2
top_face_pt_y = 0.1
top_face_pt_z = 2.5
top_face_pt = (top_face_pt_x,top_face_pt_y,top_face_pt_z)
top_face = beamInstance.faces.findAt((top_face_pt,))
top_face_region = regionToolset.Region(side1Faces=top_face)
beamModel.Pressure(name='Uniform Applied Pressure',createStepName='Apply Load',
region=top_face_region,distributionType=UNIFORM,
magnitude=10,amplitude=UNSET)
fixed_end_face_pt_x = 0.2
fixed_end_face_pt_y = 0
fixed_end_face_pt_z
= 0
fixed_end_face_pt = (fixed_end_face_pt_x,fixed_end_face_pt_y,fixed_end_face_pt_z)
fixed_end_face = beamInstance.faces.findAt((fixed_end_face_pt,))
fixed_end_face_region = regionToolset.Region(faces=fixed_end_face)
Cantilever Beam Example.py 2017/3/21 15:32
Page 1 of 2