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

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

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

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

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

所以我们需要使用Python。

下面直接给出相关Python程序。

1. #一般在ABAQUS中使用Python必须加上这三行
2. from abaqus import *
3. from abaqusConstants import *
4. from caeModules import *
5. session.journalOptions.setValues(replayGeometry=COORDINATE,recoverGeometry= COORDINATE)
6. blist = []
7. myAssembly = mdb.models['Model-1'].rootAssembly
8. #将力的时程存在txt文件里。存储格式是第一列为时间，其余每列是荷载时程。每列之间逗号分隔
9. with open('G:/zhanfang/Faster_focce.txt', 'r') as bf:
10. bRVlist = [b.strip().split(',') for b in bf]
11. with open('G:/zhanfang/Faster_focce.txt', 'r') as bf:
12. bRHlist = [b.strip().split(',') for b in bf]
13. #建立要加集中力的节点集合。因为节点排列没有顺利，所以按顺序依次建立。在我的模型里节点号为2,1,15,24,21。。。即为[x:x+1]里面的后一个值。
14. n1 = myAssembly.instances['GD-slab-1'].nodes
15. 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]
16. jiedianshu=9
17. for i in range(1,jiedianshu+1):
18. nodesi = nodes1[i-1:i]
19. F1nodelist = ''.join(['line7load',str(i)])
20. myAssembly.Set(nodes=nodesi, name=F1nodelist)
21. nodes2i = nodes2[i-1:i]
22. F2nodelist = ''.join(['line8load',str(i)])
23. myAssembly.Set(nodes=nodes2i, name=F2nodelist)
24. LoadnamelistV = ''.join(['Load8V-',str(i)])
25. LoadnamelistH = ''.join(['Load8H-',str(i)])
26. #下面这个程序块是建立垂向力的幅值
27. train250V = ()
28. strlist = ['AmpRV-',str(i)]
29. cstrlistV = ''.join(strlist)
30. #这个for循环是将荷载数据转换成ABAQUS建立幅值的形式
31. for bl in bRVlist:
32. train250V = train250V+((float(bl[0]),float(bl[i])),)
33. mdb.models['Model-1'].TabularAmplitude(name=cstrlistV, timeSpan=STEP, smooth=SOLVER_DEFAULT, data=train250V)
34. #下面这个程序块是建立横向力的幅值
35. train250H = ()
36. strlist = ['AmpRH-',str(i)]
37. cstrlistH = ''.join(strlist)
38. for bl in bRHlist:
39. train250H = train250H+((float(bl[0]),float(bl[i])),)
40. mdb.models['Model-1'].TabularAmplitude(name=cstrlistH, timeSpan=STEP,smooth=SOLVER_DEFAULT, data=train250H)
41. #建立集中力
42. region1 = myAssembly.sets[F2nodelist]
43. mdb.models['Model-1'].ConcentratedForce(name=LoadnamelistV, createStepName='Step-2',
44. region=region1, cf3=1.0, amplitude=cstrlistV, distributionType=UNIFORM,
45. field='', localCsys=None)
46. region2 = myAssembly.sets[F2nodelist]
47. mdb.models['Model-1'].ConcentratedForce(name=LoadnamelistH, createStepName='Step-2',
48. region=region1, cf1=1.0, amplitude=cstrlistH, distributionType=UNIFORM,
49. field='', localCsys=None)

tips：

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

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

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