找回密码
 注册
Simdroid-非首页
查看: 1787|回复: 28

[动态分析] abaqus中生死单元的应用

[复制链接]
发表于 2013-5-13 21:26:49 | 显示全部楼层 |阅读模式 来自 贵州贵阳
给位大神,有谁知道在abaqus中Explicit怎么设置生死单元?
发表于 2013-7-2 23:33:04 | 显示全部楼层 来自 英国
Simdroid开发平台
直接删除就是个简单的fortran小程序
        subroutine vusdfld(
c Read only variables -
     1 nblock, nstatev, nfieldv, nprops, ndir, nshr,
     2 jElem, kIntPt, kLayer, kSecPt,
     3 stepTime, totalTime, dt, cmname,
     4 coordMp, direct, T, charLength, props,
     5 stateOld,
c Write only variables -
     6 stateNew, field )
        include 'vaba_param.inc'
C
      dimension jElem(nblock), coordMp(nblock,*),
     * direct(nblock,3,3), T(nblock,3,3),
     * charLength(nblock), props(nprops),
     * stateOld(nblock,nstatev),
     * stateNew(nblock,nstatev),
     * field(nblock,nfieldv),
     * StatusOld(NBLOCK,1),
     * StatusNew(NBLOCK,1)
      character*80 cmname
C
C Local arrays from vgetvrm are dimensioned to
C maximum block size (maxblk)
C
      parameter( nrData=6 )
      character*3 cData(maxblk*nrData)
      dimension rData(maxblk*nrData), jData(maxblk*nrData),
     * strain(nblock,nrData),stress(nblock,nrData)
DO k = 1, nblock
!获取应力
CALL VGETVRM('S',RDATA,JDATA,CDATA,JSTATUS)
          SS11=rData(k*1)
          SS22=rData(k*2)
          SS33=rData(k*3)
          SS12=rData(k*4)
          SS23=rData(k*5)
          SS13=rData(k*6)
!判断是否删除
          IF (SS11.LT.100)THEN ! 1方向应力大于小于100时不删除单元, 第一个sdev=1
          stateNew(k,1)=1
          ELSEIF(SS11.GT.100)THEN ! 1方向应力大于大于100时不删除单元, 第一个sdev=0
          stateNew(k,1)=0
       END IF
       END DO
c
      RETURN
      END
具体什么条件删除自己编
在inp的material中加几行:
*Depvar, delete=1 (第1个sdev 即sdev(k,1)负责控制删除单元)
      2, (总共需要的state dependent variable数量)
*User Defined Field (usdfld子程序)
要用explicit
简单来说就这样啦
回复 1 不支持 0

使用道具 举报

发表于 2013-5-13 22:19:13 | 显示全部楼层 来自 江苏
Interfect 模块下。
Mode change. 可以指定某部分材料在某步失效或激活。

你如果是做冲裁、子弹穿透的话,就要用定义失效判断标准了。另外方法了。
回复 1 不支持 0

使用道具 举报

 楼主| 发表于 2013-5-13 22:24:36 | 显示全部楼层 来自 贵州贵阳
你这个model change 只能在standard里使用,在Explicit用不了。
回复 不支持

使用道具 举报

发表于 2013-5-13 22:47:34 | 显示全部楼层 来自 爱尔兰
需要你定义失效,可以用ABAQUS自带功能,也可以用子程序,Abaqus mannul里都有。用类似VUSDFLD跟VUMAT的话,设置depvar,默认0删除,看看帮助文档,都有
回复 不支持

使用道具 举报

发表于 2013-5-13 23:17:57 | 显示全部楼层 来自 英国
4楼说的很对 小补充下
用VUSDFLD 定义一个status变量就好 有多个status变量的话 需要指定第几个变量控制单元删除
需要在section control 里面加上element deletion(CAE 在mesh type 里面)
回复 不支持

使用道具 举报

 楼主| 发表于 2013-5-13 23:20:53 | 显示全部楼层 来自 贵州贵阳
abaqus自带功能怎么用,因为我要在动力分析中去掉一个构件,所以要用到生死单元?
回复 不支持

使用道具 举报

 楼主| 发表于 2013-5-14 13:35:54 | 显示全部楼层 来自 贵州贵阳
拜托上楼的大神,能讲的详细点吗,我不是很懂,现在要在显示动力中去除一些单元,所以我要在显示动力中设置生死单元,希望您能讲解的详细点,谢谢!
回复 不支持

使用道具 举报

发表于 2013-5-14 14:56:05 | 显示全部楼层 来自 吉林长春
想问一下楼主,利用生死单元除去的一部分会对其他部分(应力情况)造成影响么?还是只是不显示了,但是实际上还是存在的,类似于被遮盖了,对其他部分无影响?
回复 不支持

使用道具 举报

发表于 2013-5-14 22:21:24 | 显示全部楼层 来自 英国
zxg123789 发表于 2013-5-14 14:56
想问一下楼主,利用生死单元除去的一部分会对其他部分(应力情况)造成影响么?还是只是不显示了,但是实际 ...

选择element deletion的话 删除了就是彻底从模型里删除了 再后续计算中也无法恢复
回复 不支持

使用道具 举报

发表于 2013-5-14 22:29:30 | 显示全部楼层 来自 英国
guidingfa 发表于 2013-5-14 13:35
拜托上楼的大神,能讲的详细点吗,我不是很懂,现在要在显示动力中去除一些单元,所以我要在显示动力中设置 ...

需要用到fortran子程序 先简单学一下子程序的用法吧 这个很难一下讲明白。。。。
如果是指单纯删除默写单元的话 使用VUSDFLD这一个subroutine 就可以完成 具体程序自己看一下abaqus USERSUB 手册里面的VUSDFLD这一部分
我做过的是当单元应力满足一定条件就删除 这个条件就是通过子程序VUSDFLD里设置的
使用VUSDFLD需要在material里面加入depvar 也就是state变量 数量自定义
设置depvar后 每个计算节点都会
通过status变量=0或=1来控制单元是否删除

回复 不支持

使用道具 举报

 楼主| 发表于 2013-5-14 23:49:12 | 显示全部楼层 来自 贵州贵阳
如果只单纯去除一个构件中的某些单元那应该怎样设置?希望您能相信讲讲,谢谢!
回复 不支持

使用道具 举报

 楼主| 发表于 2013-5-15 22:23:28 | 显示全部楼层 来自 贵州贵阳
您有那个删除的子程序吗,能给我一份吗,我看看学习学习,行不?
回复 不支持

使用道具 举报

发表于 2013-7-2 10:38:58 | 显示全部楼层 来自 陕西西安
diesure 发表于 2013-5-14 22:29
需要用到fortran子程序 先简单学一下子程序的用法吧 这个很难一下讲明白。。。。
如果是指单纯删除默写单 ...

层主,我是新学Abaqus的,做砌体结构的时候要用到生死单元定义砂浆和砌体的接触,想问一下,可不可以具体的说一下生死单元的设置,不甚感激
回复 不支持

使用道具 举报

发表于 2013-8-25 10:19:02 | 显示全部楼层 来自 福建福州
江冰韵 发表于 2013-5-13 22:19
Interfect 模块下。
Mode change. 可以指定某部分材料在某步失效或激活。

斑竹  问下混凝土塑性损伤本构关系下如何定义混凝土失效判断标准 生死单元中  谢谢斑竹的指点...
回复 不支持

使用道具 举报

发表于 2013-10-15 15:22:13 | 显示全部楼层 来自 浙江杭州
问下混凝土塑性损伤本构关系下如何定义混凝土失效判断标准 生死单元中
回复 不支持

使用道具 举报

 楼主| 发表于 2014-12-30 14:07:10 | 显示全部楼层 来自 广东
上面的楼主,之前给的子程序是控制应力来删除单元,如果控制应变来删除单元,这个怎么编译,望各位大神解答?
回复 不支持

使用道具 举报

发表于 2016-4-28 16:07:27 | 显示全部楼层 来自 甘肃兰州
谢谢楼主大一
回复 不支持

使用道具 举报

发表于 2016-5-3 21:56:44 | 显示全部楼层 来自 浙江台州
学习了,谢谢!
回复 不支持

使用道具 举报

发表于 2016-6-3 15:32:25 | 显示全部楼层 来自 大连理工大学
wawa_0077 发表于 2013-10-15 15:22
问下混凝土塑性损伤本构关系下如何定义混凝土失效判断标准 生死单元中

层主解决了么?在混凝土塑性损伤模型下使用生死单元
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 21:32 , Processed in 0.051325 second(s), 12 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表