对于建立过大型的ABAQUS分析模型的同学们，肯定会面对繁琐的重复性建模过程。比如对上百个点施加集中力，比如建立几百个耦合约束。这些重复性操作，如果使用cae操作的话，不仅费时，而且极其容易出错。

利用Python程序来处理这些问题吧，把繁琐的工作留给程序吧，我们留下的只有创造力！！！

如下图所示的一块板，我想要在图示的节点处加上水平向和垂向的集中力。每个集中力是时间的函数，所以利用Amplitude中的Tabular来导入荷载时程。

你准备怎么做呢？是直接在cae里面来建立吗，幅值通过复制粘贴来建立？？

对于这个只有9个点的情况，当然可以这样，但是假如是有九百个呢？

所以我们需要使用Python。

下面直接给出相关Python程序。#一般在ABAQUS中使用Python必须加上这三行

from abaqus import *

from abaqusConstants import *

from caeModules import *

session.journalOptions.setValues(replayGeometry=COORDINATE,recoverGeometry= COORDINATE)

blist = []

myAssembly = mdb.models['Model-1'].rootAssembly

#将力的时程存在txt文件里。存储格式是第一列为时间，其余每列是荷载时程。每列之间逗号分隔

with open('G:/zhanfang/Faster_focce.txt', 'r') as bf:

bRVlist = [b.strip().split(',') for b in bf]

with open('G:/zhanfang/Faster_focce.txt', 'r') as bf:

bRHlist = [b.strip().split(',') for b in bf]

#建立要加集中力的节点集合。因为节点排列没有顺利，所以按顺序依次建立。在我的模型里节点号为2,1,15,24,21。。。即为[x:x+1]里面的后一个值。

n1 = myAssembly.instances['GD-slab-1'].nodes

nodes1=n1[1:2]+n1[0:1]+n1[14:15]+n1[23:24]+n1[20:21]+n1[31:32]+n1[28:29]+n1[39:40]+n1[36:37]

jiedianshu=9

for i in range(1,jiedianshu+1):

nodesi = nodes1[i-1:i]

F1nodelist = ''.join(['line7load',str(i)])

myAssembly.Set(nodes=nodesi, name=F1nodelist)

nodes2i = nodes2[i-1:i]

F2nodelist = ''.join(['line8load',str(i)])

myAssembly.Set(nodes=nodes2i, name=F2nodelist)

LoadnamelistV = ''.join(['Load8V-',str(i)])

LoadnamelistH = ''.join(['Load8H-',str(i)])

#下面这个程序块是建立垂向力的幅值

train250V = ()

strlist = ['AmpRV-',str(i)]

cstrlistV = ''.join(strlist)

#这个for循环是将荷载数据转换成ABAQUS建立幅值的形式

for bl in bRVlist:

train250V = train250V+((float(bl[0]),float(bl[i])),)

mdb.models['Model-1'].TabularAmplitude(name=cstrlistV, timeSpan=STEP, smooth=SOLVER_DEFAULT, data=train250V)

#下面这个程序块是建立横向力的幅值

train250H = ()

strlist = ['AmpRH-',str(i)]

cstrlistH = ''.join(strlist)

for bl in bRHlist:

train250H = train250H+((float(bl[0]),float(bl[i])),)

mdb.models['Model-1'].TabularAmplitude(name=cstrlistH, timeSpan=STEP,smooth=SOLVER_DEFAULT, data=train250H)

#建立集中力

region1 = myAssembly.sets[F2nodelist]

mdb.models['Model-1'].ConcentratedForce(name=LoadnamelistV, createStepName='Step-2',

region=region1, cf3=1.0, amplitude=cstrlistV, distributionType=UNIFORM,

field='', localCsys=None)

region2 = myAssembly.sets[F2nodelist]

mdb.models['Model-1'].ConcentratedForce(name=LoadnamelistH, createStepName='Step-2',

region=region1, cf1=1.0, amplitude=cstrlistH, distributionType=UNIFORM,

field='', localCsys=None)添加荷载后的板

tips：

1.对于完成某个操作的程序命令具体是怎样的，可以从打开每个模型就会自动生成的abaqus.py文件查看。这个文件记录了每步cae操作对应的命令流。

2.由于Python自身的强制缩进的机制，上面给出的程序在缩进方面需要修改。建立熟悉基本的python语言规则，在来学习将Python和abaqus结合使用。

因为上一篇的关于abaqus的专栏文章多了好多关注呢，谢谢各位的关注啊。本小硕会持续把学习过程中的点滴收获和大家分享哦。