梁单元vumat的问题

saint123 · 发表于 2010-8-22 20:52:48

本帖最后由 saint123 于 2010-8-22 21:20 编辑

这是我编的一个梁单元的vumat，施加如图所示的力，vumat里面我根据材料的单轴应力应变关系只更新了轴力

这是结果，mesh为3个小单元，不知道是不是我厚度方向的应变没有更新的缘故

当我只施加轴力的话，计算结果和手算得是差不多的
不知道是什么原因，还请高人指点

lr00011 · 发表于 2010-8-22 22:27:07

应该是没考虑全所有方向的本构关系。

saint123 · 发表于 2010-8-23 09:35:40

本帖最后由 saint123 于 2010-8-25 14:07 编辑

应该是没考虑全所有方向的本构关系。
lr00011 发表于 2010-8-22 22:27

梁单元就一个轴向的力（sigma11）和一个扭转（sigma12）的力需要更新，我这个例子上页没有扭转啊，另外，vuamt要求要自己定义straininc（*，2）和straininc（*，3），不知道怎么定义，

saint123 · 发表于 2010-8-23 12:40:25

自己顶一下期待高手出现..

lr00011 · 发表于 2010-8-23 13:55:28

梁单元的本构关系还真不知道。这个最好自己找书籍看。

jerson1020 · 发表于 2010-8-23 15:53:59

将子程序本构关系段贴上，以便发现问题

saint123 · 发表于 2010-8-23 16:51:37

将子程序本构关系段贴上，以便发现问题
jerson1020 发表于3

这个就是材料的单轴本构关系，我是根据straininc（*，1）和strainnew（*，1）的关系再根据上图曲线关系来更新stressnew（*，1）的，也就是说我把strain（*，1）和stress（*，1）直接当做上图中的epsilon和sigma，不知道这么做妥不妥

saint123 · 发表于 2010-8-23 16:52:25

本帖最后由 saint123 于 2010-8-23 20:51 编辑

......................................

jerson1020 · 发表于 2010-8-24 08:47:46

你试试位移加载

jerson1020 · 发表于 2010-8-24 08:53:17

共同探讨，我的qq154442326

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[子程序] 梁单元vumat的问题

本帖子中包含更多资源

相关帖子

本帖子中包含更多资源

.....................................