【讨论】如何用subroutine在不同的载荷步间实现材料性质的改变？

小百

可以用子程序UPOWDR给出相对密度和应力、应变等的关系，详细内容在帮助文档的d卷。

huhu

多谢小百的回复！
可是，这里的密度不是相对密度，是普通的材料的密度，不是粉末金属。

小百

密度和相对密度的关系只是相差一个基准，真正影响应力分布的也是模量的相对关系而不是模量的绝对大小。

huhu

谢了，我试试看。

huhu

谢谢小百，ych 和liuguangtiany的热心帮助！
这里的密度和模量之间的关系是已知的，难点在于：在载荷步之间传递质量密度，这一密度是由上一载荷步得出的应力和应变计算出来的。即：
密度（n+1）=应力（n）*应变（n）/密度（n）
这里n是载荷步（或时间步）
密度是质量密度，每个单元都不相同。

huhu

尊敬的ych斑竹,关于我提的这个问题,已经解决了.我用了和marc捆绑在一起的python.非常方便.它可以直接读取后处理数据.用了大概两周的时间,包括学习和编程.
用python从后处理中读数据,然后按照我所说的算法计算,再用python把结果(单元性质)放回去,整个过程加一个循环就行了.
感觉python的优点是:可以把所有的操作过程用语句表达出来,而且,可以比procedure文件多一些控制语句,象for while只类.而且,由于它可以和parameter传递参数,并可以直接操作marc的数据库,所以,对骨再造这类反复调用模型,反复读取后处理结果的题目,它比子程序灵活,比直接用fortran简单.
感谢各位高手的热心帮助和指点,让我学到了很多东西.谢谢!

liulangtiany

我想读取单元的性质是不一定要读后处理文件的，用elmvar子程序就是很方便了。我注意到marC网站上有个子程序hypla.f的用法（包括它的功能和语句），和你说的功能很相似，我想也许可以更加轻松的解决这个问题。
但是我有一点疑惑的是，你如何将结果又放入到计算中的呢
？

ych

huhu wrote:
尊敬的ych斑竹,关于我提的这个问题,已经解决了.我用了和marc捆绑在一起的python.非常方便.它可以直接读取后处理数据.用了大概两周的时间,包括学习和编程.
用python从后处理中读数据,然后按照我所说的算法计算,再用python把结果(单元性质)放回去,整个过程加一个循环就行了.
感觉python的优点是:可以把所有的操作过程用语句表达出来,而且,可以比procedure文件多一些控制语句,象for while只类.而且,由于它可以和parameter传递参数,并可以直接操作marc的数据库,所以,对骨再造这类反复调用模型,反复读取后处理结果的题目,它比子程序灵活,比直接用fortran简单.
感谢各位高手的热心帮助和指点,让我学到了很多东西.谢谢!

  
很感謝你將你的心得分享給大家知道,在不牽涉到學術機密的情況下,可否幫我們
以小程式展示python進階的用法,謝謝你.

liulangtiany

woding