求助：用PYTHON写了个小程序，不知道错在哪里

mrhenry

刚学python，不知道错在哪里了。提交job的时候貌似出错了~希望大虾们帮帮忙~~~

JingheSu

注意看错误提示:
1.
直接运行后错误:
('invalid syntax',('C:/xx.py',line160,28,"\tmdb.Job(name='scape_%s'%n,, model=modelName, description='', type=ANALYSIS,\n"))
第160行28列处有错误,而且是符号错误:仔细看应该会发现多写了一个逗号删了就好了.

2.
不过此时运行还有第二个错误:
File "C:/xx.py",line166,in <module>multiprocessingMode=DEFAULT, numDomains=1, numCpus=1)
AbaqusNameError:invalid job name
在abaqus pde中调试,再'Scape_%s'%n前一步设置断点,监视n的值发现此时
n已经不是一个简单的数字了,而是
{scape_['MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object', 'MeshNode object']}
看来是你的n在程序过程中间被重新赋值了一次.所以仿照modelName的定义,把jobName定义放在循环的最开始就行了.具体见附件.

ps.你的模型不错,算是优化设计么?

mrhenry

恩，是的，我现在只是稍微试试，还没有具体的建模，多谢了~我看看~

mrhenry

注意看错误提示:
1.
直接运行后错误:
('invalid syntax',('C:/xx.py',line160,28,"\tmdb.Job(name='scape_%s'%n,, model=modelName, description='', type=ANALYSIS,\n"))
第160行28列处有错误,而且是符号错误:仔 ...
JingheSu 发表于 2011-7-20 00:22

n确实在程序运行时发生了变化，还有2个疑问：
1、我循环了5次，但是model和job只出现最后一次的，为什么前面几步循环的modal都没有出现？看了附件里网友贴的一个小程序，感觉和我的差不多，但是他的可以产生循环几次的modal。
2、我首先在三角形截面下定义了刚体和可变形柱体是tie连接，并且在底部截面定义了约束，但是当循环变成四边形，五边形，六边形之后，tie连接的两个面发生了改变，本来应该是在底部截面的约束也跑到了侧面上，不知道该如何解决。
附件是一个网友的循环程序，与我的做个对比，希望那个网友见到莫要见怪

mrhenry

希望大家给予帮助

JingheSu

1、我循环了5次，但是model和job只出现最后一次的，为什么前面几步循环的modal都没有出现？看了附件里网友贴的一个小程序，感觉和我的差不多，但是他的可以产生循环几次的modal。
********
这个很好解决,注意看你们两个文件的区别,你的

1. Mdb()

复制代码

在循环里面所以,每次循环都新生成一个cae文件,当然把以前的覆盖掉了.
只要把这个放到循环前面就行了,

1. ns=range(3,8,1)
   
2. Mdb()
   
3. for n in ns:

复制代码

2、我首先在三角形截面下定义了刚体和可变形柱体是tie连接，并且在底部截面定义了约束，但是当循环变成四边形，五边形，六边形之后，tie连接的两个面发生了改变，本来应该是在底部截面的约束也跑到了侧面上，不知道该如何解决。
**************************
这问题很好解释,你在选取面的时候都是靠instance中面的序号来做的,那么三棱柱,五棱住和七棱柱面的数目不一样,你用同样的序号去引用,当然会出问题了:对三棱柱是底面,对于其他来说可能就是侧面了.
给你一段我查询Model_4中棱柱面编号的信息(在cae的命令行下输入,然后看到的信息)

1. >>> mdb=mdb.models['Model_4']
   
2. >>> print mdb
   
3. ({'absoluteZero': None, 'adaptiveMeshConstraints': 'Repository object', 'adaptiveMeshControls': 'Repository object', 'allFields': 'Repository object', 'amplitudes': 'Repository object', 'analyticalFields': 'Repository object', 'boundaryConditions': 'Repository object', 'calibrations': 'Repository object', 'constraints': 'Repository object', 'contactControls': 'Repository object', 'contactInitializations': 'Repository object', 'description': '', 'discreteFields': 'Repository object', 'endRestartStep': ON, 'featureOptions': 'FeatureOptions object', 'fieldOutputRequests': 'Repository object', 'filters': 'Repository object', 'globalJob': '', 'historyOutputRequests': 'Repository object', 'inputFileBlocks': 'Repository object', 'integratedOutputSections': 'Repository object', 'interactionProperties': 'Repository object', 'interactions': 'Repository object', 'keywordBlock': 'KeywordBlock object', 'lastChangedCount': 1311173432.294, 'linkedInstances': 'unknown', 'linkedParts': 'unknown', 'loads': 'Repository object', 'materials': 'Repository object', 'modelType': STANDARD_EXPLICIT, 'name': 'Model_4', 'noPartsInputFile': OFF, 'parts': 'Repository object', 'predefinedFields': 'Repository object', 'profiles': 'Repository object', 'remeshingRules': 'Repository object', 'restartIncrement': STEP_END, 'restartJob': '', 'restartStep': '', 'rootAssembly': 'Assembly object', 'sections': 'Repository object', 'shellToSolid': OFF, 'sketches': 'Repository object', 'stefanBoltzmann': None, 'steps': 'Repository object', 'surfaceProperties': 'Repository object', 'timePoints': 'Repository object', 'universalGas': None, 'waveFormulation': NOT_SET})
   
4. >>> root=mdb.rootAssembly
   
5. >>> print root
   
6. ({'allCompositeLayups': 'Repository object', 'allInternalSets': 'Repository object', 'allInternalSurfaces': 'Repository object', 'allSets': 'Repository object', 'allSurfaces': 'Repository object', 'connectorOrientations': 'ConnectorOrientationArray object', 'datum': 'Repository object', 'datums': 'Repository object', 'edges': 'EdgeArray object', 'elements': 'MeshElementArray object', 'engineeringFeatures': 'EngineeringFeatures object', 'features': 'Repository object', 'featuresById': 'Repository object', 'instances': 'Repository object', 'isLocked': 0, 'isOutOfDate': 0, 'modelName': 'Model_4', 'name': 'rootAssembly', 'nodes': 'MeshNodeArray object', 'referencePoints': 'RPRepository object', 'regenerateConstraintsTogether': ON, 'sectionAssignments': 'SectionAssignmentArray object', 'sets': 'Repository object', 'skins': 'Repository object', 'stringers': 'Repository object', 'surfaces': 'Repository object', 'timeStamp': 1739.0, 'vertices': 'VertexArray object'})
   
7. >>> inst=root.instances
   
8. >>> print inst
   
9. {'Part-1-1': 'PartInstance object', 'Part-2-1': 'PartInstance object'}
   
10. >>> myinst=inst['Part-1-1']
    
11. >>> print myinst
    
12. ({'analysisType': 'DEFORMABLE_BODY', 'cells': 'CellArray object', 'datums': 'Repository object', 'dependent': 1, 'edges': 'EdgeArray object', 'elemEdges': 'Repository object', 'elemFaces': 'Repository object', 'elementEdges': 'MeshElemEdgeArray object', 'elementFaces': 'MeshElemFaceArray object', 'elements': 'MeshElementArray object', 'faces': 'FaceArray object', 'geometryValidity': 1, 'ignoredEdges': 'IgnoredEdgeArray object', 'ignoredVertices': 'IgnoredVertexArray object', 'ips': 'IPArray object', 'name': 'Part-1-1', 'nodes': 'MeshNodeArray object', 'part': 'Part object', 'partName': 'Part-1', 'referencePoints': 'RPRepository object', 'sets': 'Repository object', 'skins': 'Repository object', 'stringers': 'Repository object', 'surfaces': 'Repository object', 'vertices': 'VertexArray object'})
    
13. >>> faces=myinst.faces
    
14. >>> print faces
    
15. ['Face object', 'Face object', 'Face object', 'Face object', 'Face object', 'Face object']
    
16. >>> highlight(faces[0])
    
17. >>> highlight(faces[1])
    
18. >>> highlight(faces[2])

复制代码

你highlight(faces[n])就是让cae把Part-1-1的编号为n的面给你高亮显示一下.利用这样的方法可以查询你自己引用的那几个面到底是底面还是侧面了.
不过,知道了是不同模型底面的具体编号,虽说勉强可以解决问题,但是还是不完美,毕竟要建好模型自己手动弄一下才知道具体的编号,要完美解决的话建议使用

1. findat()

复制代码

建议翻阅
Abaqus Scripting Reference Manual: 6.6 FaceArray object

6.6 FaceArray object

来# W2 d/ ?. j, [* x2 `8 M

lj2009jl

6# JingheSu
有段时间没用abaqus了，highlight命令现在可以直接用了吗？我记得以前是要导入什么模块的。

JingheSu

7# lj2009jl
CLI下直接输入，没问题（我没有特异导入什么模块。）
可能打开cae的过程默认导入了一些模块吧。

highlyfly

小火鸡？该不是你吧，哈哈~

JingheSu

。。。。这个，你应该肯定是认错人了。

mrhenry

。。。。这个，你应该肯定是认错人了。
JingheSu 发表于 2011-7-21 13:20

楼上是我同学。。。

highlyfly

。。。。这个，你应该肯定是认错人了。
JingheSu 发表于 2011-7-21 13:20

哇，总版主，幸会幸会，我肯定你不是我同学，哈哈~~~

highlyfly

楼上是我同学。。。
mrhenry 发表于 2011-7-21 14:33

哈哈，还真是啊~

JingheSu

误会误会