- 积分
- 0
- 注册时间
- 2013-11-5
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2014-2-25 20:39:46
|
显示全部楼层
来自 吉林长春
总结一下:
1、白金泽书上第六章前两个例子分别对应:各向同性线弹性材料(等同于相同参数的1号材料)和各向同性强化双线性弹塑性材料(等同于相同参数的3号材料);
赵海鸥书上第十五章的例子是:自定义多线性弹塑性等向硬化材料(仅对体单元)(等同于相同参数的24号材料);
2、其中各向同性线弹性材料同1号材料的对比结果几乎一致;而各向同性强化双线性弹塑性材料同3号材料的对比结果几乎一致;而自定义多线性弹塑性等向硬化材料同3号及24号材料的对比结果均有一定偏差;(对von mises应力、x向位移进行的对比)第三种情况可能是由于多段线情况的复杂性造成(暂时认为)
3、其中白金泽书上各向同性强化双线性弹塑性材料的k文件中有一处错误,即材料中NHV=1应改为NHV=2因为子程序里用到了两个历史变量,否则会出现计算错误!
4、还有就是很多高手提到的ls-dyna版本的问题,ls-971的语法与之前版本不同。如:以前的子程序(标量41号中的参数):
subroutine umat41(cm,eps,sig,hisv,dt,capa,etype,time,tmep,i,ixs,x,k,j)
971中的子程序(标量41号中的参数):
subroutine umat41(cm,eps,sig,epsp,hsv,dt1,capa,etype,tt,temper,failel,crv)
需要参见所用dyna版本的帮助进行对应的更改,否则计算出来的结果没有应力输出,并且结构的变形也不正常。
5、对于程序本身很多理论性知识还不是很懂,可能是力学知识的欠缺,应该多参考一些相关方面的书籍。
|
|