24小时热门版块排行榜    

查看: 635  |  回复: 0

supertgy

铜虫 (初入文坛)

[交流] 【求助】请大侠帮忙看看,为何已经写入了场变量,但是后处理里面却看不到创建的Step和

请大侠帮忙看看,为何已经写入了场变量,但是后处理里面却看不到创建的Step和云图?
为了完成课题,我做了一个小程序,建了一个场变量对象,然后读取了原有ODB中的PEEQ,再将PEEQ写入这个新的场变量当中。
但是后处理里面却看不到创建的Step和云图,不知是哪一步出错了,请高手们帮忙看看,不胜感激!程序如下:
# -*- coding: cp936 -*-
from abaqus import *
from abaqusConstants import *
import visualization
import math
myViewport = session.Viewport(name='X-percent',
    origin=(10, 10), width=150, height=100)
# Open the tutorial output database.
myOdb = visualization.openOdb(path='06_700.odb',readOnly=False)
# Associate the output database with the viewport.
myViewport.setValues(displayedObject=myOdb)
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#创建新的场变量输出
newStep=myOdb.Step(name='New_Step',description='haha',domain=TIME,timePeriod=1.0)
newFrame=newStep.Frame(frameId=1,frameValue=1.0,description='haha')
newField=newFrame.FieldOutput(name='DELT',description='',type=SCALAR)
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#访问原有的等效塑性应变变量
instance1 = myOdb.rootAssembly.instances['PART-1-1']
myStep = myOdb.steps['Step-2']
myFrame = myStep.frames[-1]
myPEEQ = myFrame.fieldOutputs['PEEQ']
myPEEQV = myPEEQ.values
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#重新将等效塑性应变写入新建的场变量中
elements=instance1.elements
elementLabels=[0]*len(elements)
for i in range(len(elements)):
    elementLabels=elements.label
elementLabels.sort()#给单元号排序
elementValues=[0]*len(elements)
for v in myPEEQV:#用数组储存积分点处应变数值
    find = elementLabels.index(v.elementLabel)
    elementValues[find]=v.data#与单元号相对应
elementValuesN=a=[[0 for x in range(1)] for y in range(len(elements))]#创建二维序列储存积分点处应变数值
for u in range(len(elements)):
    elementValuesN[0]=elementValues
newField.addData(position=INTEGRATION_POINT, instance=instance1, labels=elementLabels, data=elementValuesN)#将结果写入新建场变量
myOdb.save()
print newField.values[9].data#为了检查是不是将数据真的写入了新建变量中,用这一行程序读取任意一个变量值

[ Last edited by supertgy on 2010-7-26 at 22:06 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 supertgy 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见