【讨论】怎样访问ABAQUS的原代码？

ilovenili

请大虾指教，ABAQUS有没有对用户开放原代码,
怎样才能看到ABAQUS的原代码.
举一个简单的例子,我想知道ABAQUS中关于画圆
的PYTHON代码该怎么办?
请指教!!

blast

from abaqus import *
from abaqusConstants import *
session.Viewport(name='Viewport: 1', origin=(0.0, 0.0), width=246.796875,  
    height=172.265625)
session.viewports['Viewport: 1'].makeCurrent()
session.viewports['Viewport: 1'].maximize()
from driverUtils import executeOnCaeStartup
executeOnCaeStartup()
Mdb()
#: A new model database has been created.
#: The model "Model-1" has been created.
session.viewports['Viewport: 1'].setValues(displayedObject=None)
import part
import regionToolset
import displayGroupMdbToolset as dgm
import sketch
session.viewports['Viewport: 1'].setValues(displayedObject=None)
s = mdb.models['Model-1'].Sketch(name='__profile__', sheetSize=200.0)
g, v, d = s.geometry, s.vertices, s.dimensions
s.sketchOptions.setValues(sheetSize=200.0, gridSpacing=5.0, grid=ON,  
    gridFrequency=2, constructionGeometry=ON, dimensionTextHeight=5.0,  
    decimalPlaces=2)
s.setPrimaryObject(option=STANDALONE)
s.CircleByCenterPerimeter(center=(0.0, 0.0), point1=(30.0, 0.0))
#
#在Sketch中，画中心为（0，0），圆周上有一点位于（30，0）的圆
#
mdb.models['Model-1'].sketches.changeKey(fromName='__profile__',  
    toName='Sketch-1')
s.unsetPrimaryObject()

ilovenili

遇到牛人了，这是自己编写的，还是在哪里找的？
清指教！

ilovenili

不过看起来好象不是画圆的ABAQUS原代码！！
呵呵，见笑！

blast

是PYTHON的代码
可以满足画圆要求:I
我没有说是fortran的源代码:I
感觉不是很清楚你想要什么:I，从源代码、PYTHON级水平建模吗
我肯定不是牛仁了，刚自学几天^^b

ilovenili

因为ABAQUS是用PYTHON编写的，
那么我想要PYTHON 的原代码，因为我想
在其中加一个椭圆程序，不知可行么？

ilovenili

其实准确的说，我是想知道画圆函数
CircleByCenterPerimeter（）的源代码！
呵呵，不知道我现在表述清楚了没有：（

blast

要获得初步的代码很简单，
每次abaqus启动都会修改当前目录下的abaqus.rpy文件
这个文件实际就是PYTHON语言写成的
abaqus.rpy纪录每次启动计算机进行的所有步骤
我上面地代码就是直接从abaqus.rpy中获取的
  
椭圆的函数不知道有没有
不过有苯办法
虽然方法有点苯
用直线连出来的
不过，只要能分网，应该不会影响分网后的计算吧
你可以先找到画短线的方法，然后再一条一条连出个椭圆（不知道有没有直接画椭圆的函数:I)
  
我以前做的椭圆函数找不到了……不过那个没几行相信你自己可以做一个
  
PS:提醒一下PYTHON的计算功能不是很好，注意除错，这点也许画椭圆用得着

ilovenili

能不能够再深入一点，知道里面的函数是怎么编写的？

blast

象下面的方式就可以调用\写函数

1. a=f(30)
2. #a=f(30)是调用函数f，令a=60
3. def f(x):
4.   return x*2
5.   
6. def fax1(x):
7.   return x*x
8.   
9. c=fax1(20)
10. #c=fax1(20)是调用函数f，令c=400
11.   
12. #调用函数和函数位置没有先后分别
13. #python是完全靠缩进判断块的，注意缩进定义函数内容

复制代码

ilovenili

看来大虾是一位PYTHON高手
可以取取经么？看哪本书较好

blast

上面代码缩进不对
上传附件……

blast

高手^6^…………:I:I
  
市面上没见过针对abaqus的python教程
其他的随便看了。真没什么可推荐的= =

cafe0519

ilovenili wrote:
因为ABAQUS是用PYTHON编写的，
那么我想要PYTHON 的原代码，因为我想
在其中加一个椭圆程序，不知可行么？

  
;):(
所有的PY都是差不多一样的,把ABAQUS的PY拿到MARC玩玩看?!

ilovenili

可是我不懂MARC哦
呵呵，失望吧！