Abaqus python实例003 sat文件导入并运算 立方体模型的导入(与实例

分类: Abaqus实例Python应用 |
下边代码分两部分,“{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{”之前的是python代码部分,“{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{”和“}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}”之间的部分是模型的sat文件,代码拷贝到txt文档中修改后缀为“.sat“即可。别忘了修改sat文件的存放目录,代码中有注释。注意查看
#!/user/bin/python
# -* - coding:UTF-8 -*-
from abaqus import *
from abaqusConstants import *
#session.Viewport(name='Viewport: 1', origin=(0.0, 0.0),
width=190.46875,height=140.649993896484)
myViewport=session.Viewport(name='Viewport: 1',origin=(0.0,
0.0), width=176.25, height=130.5)
session.viewports['Viewport: 1'].makeCurrent()
#窗口最大显示
session.viewports['Viewport:
1'].maximize()
from caeModules import *
from driverUtils import executeOnCaeStartup
executeOnCaeStartup()
session.viewports['Viewport:
1'].partDisplay.geometryOptions.setValues(
Mdb()
#: A new model database has been created.
#: The model "Model-1" has been created.
session.viewports['Viewport:
1'].setValues(displayedObject=None)
acis = mdb.openAcis('D:/AbaqusTemp/Cub/cub.sat',
scaleFromFile=OFF)
#修改sat文件的目录
#定义保存cae文件的目录
FilePath ='D:/AbaqusTemp/Cub/'
#定义保存cae文件名
FileName ='cub'
mdb.models['Model-1'].PartFromGeometryFile(name='Cube',
geometryFile=acis,
p = mdb.models['Model-1'].parts['Cube']
session.viewports['Viewport:
1'].setValues(displayedObject=p)
session.viewports['Viewport:
1'].partDisplay.setValues(mesh=ON)
session.viewports['Viewport:
1'].partDisplay.meshOptions.setValues(
session.viewports['Viewport:
1'].partDisplay.geometryOptions.setValues(
mm = mdb.models['Model-1']
#导入material模块
import material
#创建材料材料名为Steel
mySteel=mm.Material(name='Steel')
#定义弹性材料属性,杨氏模量为209E3,泊松比为0.3
elasticProperties=(209E3,0.3)
mySteel.Elastic(table=(elasticProperties,))
#导入section模块
import section
#创建实体截面
mySection=mm.HomogeneousSolidSection(name='cubeSection',material='Steel',thickness=1.0)
#为部件分配截面属性
region=(p.cells,)
p.SectionAssignment(region=region,sectionName='cubeSection')
#导入assembly模块
import assembly
#创建部件实例
myAssembly=mm.rootAssembly
myInstance=myAssembly.Instance(name='cubeInstance',part=p,dependent=OFF)
#导入step模块
import step
#在初始分析步Initial之后创建一个分析步。静力分析的时间为1.0,初始增量为0.1
mm.StaticStep(name='cubeLoad',previous='Initial',timePeriod=1.0,initialInc=0.1,description='Load
the top of the cube')
#导入load模块
import load
#通过坐标找出端部所在的面
endFaceCenter=(5.0,5.0,0)
endFace=myInstance.faces.findAt((endFaceCenter,))
#在端部创建固定端约束
endRegion=(endFace,)
mm.EncastreBC(name='Fixed',createStepName='cubeLoad',region=endRegion)
#通过坐标找到上表面
topFaceCenter=(5.0,5.0,10.0)
topFace=myInstance.faces.findAt((topFaceCenter,))
#在模型的上表面施加压力荷载
topSurface=((topFace,SIDE1),)
mm.Pressure(name='Pressure',createStepName='cubeLoad',region=topSurface,magnitude=5e6)
#单位为Pa
#导入mesh模块
import mesh
#为部件实例指定单元类型
region=(myInstance.cells,)
elemType=mesh.ElemType(elemCode=C3D8I,elemLibrary=STANDARD)
myAssembly.setElementType(regions=region,elemTypes=(elemType,))
#为部件实例撒种子
myAssembly.seedPartInstance(regions=(myInstance,),size=2.5)
#为部件实例划分网格
myAssembly.generateMesh(regions=(myInstance,))
#显示划分网格后的模型
myViewport.assemblyDisplay.setValues(mesh=ON)
myViewport.assemblyDisplay.meshOptions.setValues(meshTechnique=ON)
myViewport.setValues(displayedObject=myAssembly)
#导入job模块
import job
#为模型创建并提交分析作业
jobName='cube_tutorial'
myJob=mdb.Job(name=jobName,model='Model-1',description='Cantilever
cube tutorial')
#等待分析作业完成
myJob.submit()
myJob.waitForCompletion()
print'分析已顺利完成,下面进行后处理。'
#导入Visualization模块
import visualization
#打开输出数据库,显示默认的云图
myOdb=visualization.openOdb(path=jobName+'.odb')
myViewport.setValues(displayedObject=myOdb)
myViewport.odbDisplay.display.setValues(plotState=CONTOURS_ON_DEF)
myViewport.odbDisplay.commonOptions.setValues(renderStyle=FILLED)
#设置显示背景为白色
session.graphicsOptions.setValues(backgroundColor='#FFFFFF',
#将Mises云图输出为PNG格式的文件
session.printToFile(fileName='Mises',format=PNG,canvasObjects=(myViewport,))
print'文件Mises.png保存于工作目录下,请查看!'
#Save .cae file
mdb.saveAs(pathName=FilePath + FileName + '.cae')
http://s10/mw690/001rwgmJzy6NpdMKbUd79&690 python实例003 sat文件导入并运算 立方体模型的导入(与实例" TITLE="Abaqus python实例003 sat文件导入并运算 立方体模型的导入(与实例" />
{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{
700 0 1 0
33 Autodesk Inventor 19.0.15900.0000 19 ASM 220.0.0.5900 NT 24
Tue Oct 14 13:16:44 2014
1 9.9999999999999995e-007 1e-010
body $-1 -1 $-1 $1 $-1 $-1 #
lump $-1 -1 $-1 $-1 $2 $0 #
shell $-1 -1 $-1 $-1 $-1 $3 $-1 $1 #
face $-1 -1 $-1 $4 $5 $2 $-1 $6 forward single #
face $-1 -1 $-1 $7 $8 $2 $-1 $9 forward single #
loop $-1 -1 $-1 $-1 $10 $3 #
plane-surface $-1 -1 $-1 0 -6.6613381477509392e-015 10 0 -1
-6.6613381477509392e-016 0 6.6613381477509392e-015 -10 forward_v I
I I I #
face $-1 -1 $-1 $11 $12 $2 $-1 $13 forward single #
loop $-1 -1 $-1 $-1 $14 $4 #
plane-surface $-1 -1 $-1 0 9.9999999999999929 10 0 0 1 0 -10 0
forward_v I I I I #
coedge $-1 -1 $-1 $15 $16 $17 $18 forward $5 $-1 #
face $-1 -1 $-1 $19 $20 $2 $-1 $21 forward single #
loop $-1 -1 $-1 $-1 $22 $7 #
plane-surface $-1 -1 $-1 0 10 0 0 1 6.6613381477509392e-016 0
-6.6613381477509392e-015 10 forward_v I I I I #
coedge $-1 -1 $-1 $23 $24 $25 $26 forward $8 $-1 #
coedge $-1 -1 $-1 $27 $10 $28 $29 forward $5 $-1 #
coedge $-1 -1 $-1 $10 $27 $23 $30 reversed $5 $-1 #
coedge $-1 -1 $-1 $25 $31 $10 $18 reversed $32 $-1 #
edge $-1 -1 $-1 $33 0 $34 1 $10 $35 forward @7 unknown #
face $-1 -1 $-1 $36 $37 $2 $-1 $38 forward single #
loop $-1 -1 $-1 $-1 $39 $11 #
plane-surface $-1 -1 $-1 0 0 0 0 0 -1 0 10 0 forward_v I I I I
#
coedge $-1 -1 $-1 $40 $41 $42 $43 forward $12 $-1 #
coedge $-1 -1 $-1 $44 $14 $16 $30 forward $8 $-1 #
coedge $-1 -1 $-1 $14 $44 $40 $45 reversed $8 $-1 #
coedge $-1 -1 $-1 $42 $17 $14 $26 reversed $32 $-1 #
edge $-1 -1 $-1 $46 0 $33 1 $14 $47 forward @7 unknown #
coedge $-1 -1 $-1 $16 $15 $48 $49 reversed $5 $-1 #
coedge $-1 -1 $-1 $39 $50 $15 $29 reversed $20 $-1 #
edge $-1 -1 $-1 $34 0 $51 1 $28 $52 forward @7 unknown #
edge $-1 -1 $-1 $33 0 $53 1 $23 $54 forward @7 unknown #
coedge $-1 -1 $-1 $17 $42 $39 $55 reversed $32 $-1 #
loop $-1 -1 $-1 $-1 $31 $36 #
vertex $-1 -1 $-1 $26 $56 #
vertex $-1 -1 $-1 $55 $57 #
straight-curve $-1 -1 $-1 0 -6.6613381477509392e-015 10 0
6.6613381477509392e-015 -10 F 0 F 1 #
face $-1 -1 $-1 $-1 $32 $2 $-1 $58 reversed single #
loop $-1 -1 $-1 $-1 $59 $19 #
plane-surface $-1 -1 $-1 10 4.9999999999999964 5 1 0 0 0 0 -10
forward_v I I I I #
coedge $-1 -1 $-1 $60 $28 $31 $55 forward $20 $-1 #
coedge $-1 -1 $-1 $61 $22 $24 $45 forward $12 $-1 #
coedge $-1 -1 $-1 $22 $61 $60 $62 reversed $12 $-1 #
coedge $-1 -1 $-1 $31 $25 $22 $43 reversed $32 $-1 #
edge $-1 -1 $-1 $63 0 $46 1 $22 $64 forward @7 unknown #
coedge $-1 -1 $-1 $24 $23 $65 $66 reversed $8 $-1 #
edge $-1 -1 $-1 $46 0 $67 1 $40 $68 forward @7 unknown #
vertex $-1 -1 $-1 $43 $69 #
straight-curve $-1 -1 $-1 0 9.9999999999999929 10 0 -10 0 F 0
F 1 #
coedge $-1 -1 $-1 $59 $65 $27 $49 forward $37 $-1 #
edge $-1 -1 $-1 $53 0 $51 1 $27 $70 forward @7 unknown #
coedge $-1 -1 $-1 $28 $60 $59 $71 reversed $20 $-1 #
vertex $-1 -1 $-1 $29 $72 #
straight-curve $-1 -1 $-1 0 0 0 10 0 0 I I #
vertex $-1 -1 $-1 $30 $73 #
straight-curve $-1 -1 $-1 0 -6.6613381477509392e-015 10 10 0 0
I I #
edge $-1 -1 $-1 $34 0 $63 1 $39 $74 forward @7 unknown #
point $-1 -1 $-1 0 -6.6613381477509392e-015 10 #
point $-1 -1 $-1 0 0 0 #
plane-surface $-1 -1 $-1 0 4.9999999999999964 5 1 0 0 0 0 -10
forward_v I I I I #
coedge $-1 -1 $-1 $75 $48 $50 $71 forward $37 $-1 #
coedge $-1 -1 $-1 $50 $39 $41 $62 forward $20 $-1 #
coedge $-1 -1 $-1 $41 $40 $75 $76 reversed $12 $-1 #
edge $-1 -1 $-1 $63 0 $77 1 $60 $78 forward @7 unknown #
vertex $-1 -1 $-1 $55 $79 #
straight-curve $-1 -1 $-1 0 10 0 0 -6.6613381477509392e-015 10
F 0 F 1 #
coedge $-1 -1 $-1 $48 $75 $44 $66 forward $37 $-1 #
edge $-1 -1 $-1 $67 0 $53 1 $44 $80 forward @7 unknown #
vertex $-1 -1 $-1 $45 $81 #
straight-curve $-1 -1 $-1 0 9.9999999999999929 10 10 0 0 I I
#
point $-1 -1 $-1 0 9.9999999999999929 10 #
straight-curve $-1 -1 $-1 10 -6.6613381477509392e-015 10 0
6.6613381477509392e-015 -10 I I #
edge $-1 -1 $-1 $51 0 $77 1 $50 $82 forward @7 unknown #
point $-1 -1 $-1 10 0 0 #
point $-1 -1 $-1 10 -6.6613381477509392e-015 10 #
straight-curve $-1 -1 $-1 0 0 0 0 10 0 F 0 F 1 #
coedge $-1 -1 $-1 $65 $59 $61 $76 forward $37 $-1 #
edge $-1 -1 $-1 $77 0 $67 1 $61 $83 forward @7 unknown #
vertex $-1 -1 $-1 $62 $84 #
straight-curve $-1 -1 $-1 0 10 0 10 0 0 I I #
point $-1 -1 $-1 0 10 0 #
straight-curve $-1 -1 $-1 10 9.9999999999999929 10 0 -10 0 I I
#
point $-1 -1 $-1 10 9.9999999999999929 10 #
straight-curve $-1 -1 $-1 10 0 0 0 10 0 I I #
straight-curve $-1 -1 $-1 10 10 0 0 -6.6613381477509392e-015
10 I I #
point $-1 -1 $-1 10 10 0 #
End-of-ACIS-data
}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}