ABAQUS读取节点集坐标

ABAQUS读取节点集坐标
ABAQUS读取节点集坐标

在abaqus中使用python实现的功能

(2010-02-06 15:04:10)

转载▼

标签:

杂谈

功能一:实行提交多个job的功能。

对象:Job object

使用:在源文件开始写上import job,源程序用mdb.jobs[name] 使用名字为name的job对象。

建立一个job对象的方法:

λ利用已有的inp文件中建立job:mdb.JobFromInputFile()

λ利用已有的cae中建立job: Job(...)

建议用第一种方法。

设定参数的方法:

?利用第一种方法建立job的时候,可以设定很多的参数,比如type,queue,userSubroutine等。格式:mdb.JobFromInputFile(name=,inputFile=,type=,queue=,userSubroutine=,…….)。

?也可以先建立一个job,然后利用job对象的setValues来设定参数,格式:

job.setValues(type=,queue=,userSubroutine=,…….)。

一个简单的例子:

文件:job.py

from abaqusConstants import *

import job

mdb.JobFromInputFile(name='job-1-1',inputFileName='Job-1.inp')

#基于inp文件Job-1.inp建立名称为job-1-1的job

mdb.jobs['job-1-1'].setValues(waitMinutes=1)

#设定参数

mdb.jobs['job-1-1'].submit()

#提交任务

mdb.jobs['job-1-1'].waitForCompletion()

运行:

在cmd下面运行:Abaqus cae nogui=job.py

如果是多个job,同样道理了,不多说了。

功能二:

后处理,提取需要的数据,形成可以用其他软件处理的文件

功能:提取odb文件中某个set中的数据(可以是应力、应变和位移、坐标等),建立一个外部文件,把提取的结果写到这个文件中,利用tecplot处理。

1)在Odb对象中提取场变量:

odb-》steps-》frams-》fieldoutputs【变量名称】

具体odb对象中的各个成员如下图

2)在场变量中选取所需要set的变量:

Odb-》rootAssembly-》Sets

或者odb-》rootAssembly-》instances-》Sets

上面两种方法取决了你在inp文件种是在assembly定义了Set还是在Instance中定义了Set。具体的如下图

文件:plot.py

from odbAccess import *

from abaqusConstants import *

import string

print 'begin abaqus python'

print 'today is 10-29'

print 'this code is for 3node_sin.odb'

odb = openOdb(path='3node_sin.odb')

myAssembly = odb.rootAssembly

f=open('plot/3node_sin1.dat','w')

f.write('TITLE = Example: Simple XY Plot\n')

f.write('VARIABLES = "X-Coordinate", "Y-Coordinate"\n')

#建立一个tecplot的dat文件

i=1

while i<100:

Frame = odb.steps['Step-1'].frames

print odb.steps['Step-1'].frames

coordinate=Frame.fieldOutputs['COORD']

#提取节点坐标数据

center = odb.rootAssembly.nodeSets['Set-1']

centerCoordinate = coordinate.getSubset(region=center)

centerValues = centerCoordinate.values

#提取Set-1集合中的节点坐标

count=len(centerValues)

s='ZONE T=" '+str(i)+' Zone", I='+str(count)+', F=POINT\n'

f.write(s)

for v in centerValues:

s2=str(v.data[0])

s3=str(v.data[1])

f.write(s2)

f.write('\t')

f.write('\t')

f.write(s3)

f.write('\t')

f.write('\n')

print i;i=i+50

#写入tecplot的dat文件中

运行:

Cmd中运行:abaqus scrip=plot.py

之所以用python因为abaqus的界面是基于python的,用它来进行批处理或者后处理是比较不错的,希望能和大家分享abaqus的更多功能。

相关主题
相关文档
最新文档