Abaqus 在命令提示符中中止脚本 - IndexError: -1
我遇到了一个错误,提示 IndexError: -1
我在命令提示符下运行一个Python脚本时出现这个问题。在Abaqus CAE中运行时,一切正常。但是当我在没有图形界面的命令提示符下运行时,就会出现这个错误。
我能完成这个任务,但我的脚本还会打开输出数据库。错误就是在这个步骤出现的。
我试着在网上搜索这个问题,但没能找到有用的信息。这在Python中一般意味着什么呢?为什么我的脚本在Abaqus中可以正常工作,但在命令提示符下却只完成了一半?
谢谢!
from time import gmtime, strftime
import time, os
from abaqus import *
from part import *
from material import *
from section import *
from assembly import *
from step import *
from interaction import *
from load import *
from mesh import *
from job import *
from sketch import *
from visualization import *
from connectorBehavior import *
from odbAccess import *
mdb = openMdb(pathName='X:\MinPreload-VaryingBC')
modelName = 'minPreload'
cMdb = mdb.models[modelName]
date = strftime("%Y-%d-%b", gmtime())
logPath_bu = 'C:\Users\yyt\maxMin-%s.txt' %date
logObj_bu = open(logPath_bu, 'a')
TopOfThePage = 'JobName\tCFNM at Preload\tTime Finished\tPreload\n\n'
logObj_bu.write(TopOfThePage)
logObj_bu.close()
output = open('C:\yyt\Output-%s.txt' %date, 'a')
output.write('Job \t\t\t'+'PS723\t\t'+'PS721\t\t'+'PL\t\t'+'FRIC\n\n')
output.close()
jobName = ""
maxCf = 0.0
#Temperature
minTemp = -49
maxTemp = 101
temperature = [minTemp, maxTemp]
bcTemp = 'Temperature'
#Friction
friction = [0.1, 0.1, 0.2, 0.2, 0.3, 0,3, 0.4, 0,4, 0.5, 0.5]
interactionProperty = 'Contact Seal Coupler'
#PressureSwitcher
Pressure = True
internalPressure = 76
cMdb.loads['IP'].resume('Pressure')
internalPressBC = cMdb.loads['IP']
setInternalPressureBC = internalPressBC.setValues(magnitude=internalPressure)
externalPressure = 33
cMdb.loads['ExtP'].resume('Pressure')
externalPressBC = cMdb.loads['ExtP']
setExternalPressureBC = externalPressBC.setValues(magnitude=externalPressure)
# Preload
adjustedPreload = 0
PLipht = 10
PLiplt = 10
PLextpht = 10
PLextplt = 10
for temp in temperature:
output = open('C:\yyt\Output-%s.txt' %date, 'a')
output.write('\n')
output.close()
tempBC = cMdb.boundaryConditions[bcTemp]
setTemp = tempBC.setValues(magnitude=temp)
if temp == minTemp:
adjustedPreload = 0.10315
else:
adjustedPreload = -0.21288
for fric in friction:
tangProp = cMdb.interactionProperties[interactionProperty].tangentialBehavior
tangProp.setValues(table=((fric,),))
if Pressure:
cMdb.loads['IP'].resume('Pressure')
cMdb.loads['ExtP'].suppress('Pressure')
press = 'IP'
else:
cMdb.loads['ExtP'].resume('Pressure')
cMdb.loads['IP'].suppress('Pressure')
press = 'ExtP'
if temp == maxTemp and Pressure:
PLi = PLipht
elif temp == minTemp and Pressure:
PLi = PLiplt
elif temp == maxTemp and Pressure:
PLi = PLextpht
elif temp == minTemp and Pressure:
PLi = PLextplt
for PL in range (PLi, 100):
axPL = float(PL)/100
preloadBC = cMdb.boundaryConditions['axial_move']
setPreloadBC = preloadBC.setValuesInStep('Pre', u2=axPL)
clampBC = cMdb.boundaryConditions['clamp']
modifyBottomCoupler = clampBC.setValuesInStep('Temp', -adjustedPreload )
setPreloadBC = preloadBC.setValuesInStep('Temp', u2=(axPL+adjustedPreload))
#Create a job
fricInt = fric*10
jobName = date + '_%dC-0_%df-%s-0_%dPLl' %(temp, fricInt, press, PL)
mdb.Job(model=modelName, name=jobName)
#submit Job
mdb.jobs[jobName].submit(consistencyChecking=OFF)
#Wait for job to complete
status = str(mdb.jobs[jobName].messages[-1].type)
mdb.jobs[jobName].waitForCompletion()
#Collect data from job
lastStepName = 'Pressure'
odb = session.openOdb(name=jobName+'.odb')
region = odb.steps[lastStepName].historyRegions['NodeSet Z000002']
region2 = odb.steps[lastStepName].historyRegions['NodeSet Z000001']
CFNM1 = region.historyOutputs['CFNM ASSEMBLY_SLAVEBOT/ASSEMBLY_MASTERBOT'].data
CFNM2 = region2.historyOutputs['CFNM ASSEMBLY_SLAVESURFTOP/ASSEMBLY_MASTERTOP'].data
for i, var in enumerate(CFNM1):
if i == len(CFNM1) - 1:
maxCf1 = var[-1]
print str(maxCf1) + 'CFNM1 724/723'
for i, var in enumerate(CFNM2):
if i == len(CFNM2) - 1:
maxCf2 = var[-1]
print str(maxCf2) + 'CFNM2 722/721'
output = open('C:\yyt\Output-%s.txt' %date, 'a')
output.write(jobName+'\t\t'+str(maxCf1)+'\t\t'+str(maxCf2)+'\t\t'+'0.'+str(PL)+'\t\t'+str(fric)+'\n')
output.close()
if maxCf1 >= 10000 and maxCf2 >= 10000:
time_done = strftime("%H:%M:%S", gmtime())
logObj_bu = open(logPath_bu, 'a')
logObj_bu.write(jobName+'\t'+' CFNM1 = ' + str(maxCf1)+' CFNM2 = ' + str(maxCf2)+'\t' + time_done + '%d\n\n' %PL)
logObj_bu.close()
Pressure = not Pressure
if temp == maxTemp and Pressure:
PLipht = axPL-0.02
elif temp == minTemp and Pressure:
PLiplt = axPL-0.02
elif temp == maxTemp and Pressure:
PLextpht = axPL-0.02
elif temp == minTemp and Pressure:
PLextplt = axPL-0.02
break
我得到的结果
End Abaqus/Standard Analysis
Abaqus JOB 2014-06-Mar_-49C-0_1f-IP-0_10PLl COMPLETED
IndexError: -1
Abaqus Error: cae exited with an error.
1 个回答
1
你可能会遇到三个关于 IndexError: -1 的常见问题。
[...]
status = str(mdb.jobs[jobName].messages[-1].type)
[...]
for i, var in enumerate(CFNM1):
if i == len(CFNM1) - 1:
maxCf1 = var[-1]
print str(maxCf1) + 'CFNM1 724/723'
for i, var in enumerate(CFNM2):
if i == len(CFNM2) - 1:
maxCf2 = var[-1]
print str(maxCf2) + 'CFNM2 722/721'
[...]
IndexError 这个错误说明你在尝试访问一个列表或序列中不存在的元素。因为你使用的是 -1 这个索引,这意味着你想访问列表或序列中的最后一个元素。所以在访问任何元素之前,你应该先检查一下列表是否为空,比如:
[...]
for i, var in enumerate(CFNM2):
if i == len(CFNM2) - 1 and var:
maxCf2 = var[-1]
print str(maxCf2) + 'CFNM2 722/721'
[...]
注意 if 语句末尾的 and var
。