24小时热门版块排行榜

返回列表

supertgy

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 143.1
散金: 5
帖子: 38
在线: 6.8小时
虫号: 791447
注册: 2009-06-10
性别: GG
专业: 机械工程

[交流] 【求助】请大侠帮忙看看，为何已经写入了场变量，但是后处理里面却看不到创建的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 ]

回复此楼

» 猜你喜欢

1楼 2010-07-26 22:02:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 supertgy 的主题更新

返回列表