|
楼主 |
发表于 2016-12-14 17:56:56
|
显示全部楼层
本帖最后由 淋雨后的晴朗 于 2016-12-14 17:59 编辑
我能再问你个问题吗?又多写了几行,写到下面的时候又出错误了
from abaqus import *
import math
filePath='E:/Temp/'
Emod = [9000 for i in range(20)]
EdgeEltList = map(int, open(r'E:\Temp\fanyan-20ceng-force-model-EdgeEltList.txt').readline().split(','))
b = zip(EdgeEltList,range(20))
DICx = map(float, open(r'E:\Temp\fanyan-20ceng-force-model-DIC-le11.txt').readline().split(',')) #####################
DICy = map(float, open(r'E:\Temp\fanyan-20ceng-force-model-DIC-le22.txt').readline().split(',')) #####################
DICxy= map(float, open(r'E:\Temp\fanyan-20ceng-force-model-DIC-le12.txt').readline().split(',')) #####################
d=list(DICx)
e=list(DICy)
f=list(DICxy)
DIC = [math.sqrt(d**2+e**2+f**2) for i in range(20)]
inputTxt = open(filePath + '20ceng-force-model.inp').readlines()
fileNameNew='fanyan-20ceng-0'
for i in range(20):
inputTxt[525 + i*3] = str(Emod)+', 0.3\n' ##modify some lines in inp file
open(filePath + fileNameNew + '.inp','w').writelines(inputTxt)
def mdb():
session.Viewport(name='Viewport: 1', origin=(0.0, 0.0), width=100, height=100)
sv1 = session.viewports['Viewport: 1']
sv1.makeCurrent()
sv1.maximize()
executeOnCaeStartup()
for iLoad in range(10):
fileNameNew='fanyan-20ceng-' + str(iLoad)
mdb.ModelFromInputFile(name='Model-1',
inputFileName=filePath + fileNameNew + '.inp')
mdb.Job(name=fileNameNew, model='Model-1', description='', type=ANALYSIS,
atTime=None, waitMinutes=0, waitHours=0, queue=None, memory=60,
memoryUnits=PERCENTAGE, getMemoryFromAnalysis=True,
explicitPrecision=SINGLE, nodalOutputPrecision=SINGLE, echoPrint=OFF,
modelPrint=OFF, contactPrint=OFF, historyPrint=OFF, userSubroutine='',
scratch='', multiprocessingMode=DEFAULT, numCpus=1)
mdb.jobs[fileNameNew].waitForCompletion()
|
|