| 查看: 3179 | 回复: 5 | ||
| 本帖产生 1 个 EPI ，点击这里进行查看 | ||
ppzzyy金虫 (正式写手)
|
[求助]
Abaqus的子程序问题
|
|
| 我在用Abaqus的子程序vumat，有人知道怎样得到积分点编号吗？ |
回复此楼» 猜你喜欢
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然评审意见
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
-
有没有仍没收到信息的
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
-
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ABAQUS中的用户子程序fortran编写问题
已经有4人回复
- 求助！！！！在线等！！！怎样从应力-应变曲线上找出断裂点
已经有8人回复
- 求助}abaqus循环计算问题
已经有4人回复
- 一个abaqus装了两天还没搞定，郁闷了
已经有11人回复
- ABAQUS的UMAT子程序学习心得
已经有243人回复
- 【讨论】有同仁用abaqus模拟材料的热力耦合吗？
已经有8人回复
- 【求助】abaqus里面能用梁单元做复合材料的分析吗
已经有4人回复
- 【讨论】子程序Umat讨论帖
已经有12人回复
- 【资源】浅谈ABAQUS用户子程序【已搜索无重复】
已经有146人回复
ppzzyy
金虫 (正式写手)
- EPI: 1
- 应助: 5 (幼儿园)
- 金币: 1230.5
- 散金: 615
- 红花: 1
- 帖子: 316
- 在线: 142.4小时
- 虫号: 1041469
- 注册: 2010-06-14
- 性别: MM
- 专业: 结构工程
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kingslcq: 金币+10, EPI+1, 鼓励讨论! 这个议题甚好! 2013-02-16 17:41:56
kingslcq: 金币+10, EPI+1, 鼓励讨论! 这个议题甚好! 2013-02-16 17:41:56
二楼找到的程序，我也搜到了。只是看不太明白，刚开始学ABAQUS 我搜到的完整版本如下： c User subroutine VUMAT subroutine vumat ( c Read only - * jblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal, * stepTime, totalTime, dt, cmname, coordMp, charLength, * props, density, strainInc, relSpinInc, * tempOld, stretchOld, defgradOld, fieldOld, * stressOld, stateOld, enerInternOld, enerInelasOld, * tempNew, stretchNew, defgradNew, fieldNew, c Write only - * stressNew, stateNew, enerInternNew, enerInelasNew ) c include 'vabua_param.inc' c dimension jblock(*), props(nprops),density(*), coordMp(*), 1 charLength(*), strainInc(*), 2 relSpinInc(*), tempOld(*), 3 stretchOld(*), 4 defgradOld(*), 5 fieldOld(*), stressOld(*), 6 stateOld(*), enerInternOld(*), 7 enerInelasOld(*), tempNew(*), 8 stretchNew(*), 9 defgradNew(*), 1 fieldNew(*), 2 stressNew(*), stateNew(*), 3 enerInternNew(*), enerInelasNew(*) c character*80 cmname parameter ( * i_umt_nblock = 1, * i_umt_npt = 2, * i_umt_layer = 3, * i_umt_kspt = 4, * i_umt_noel = 5 ) call vumatXtrArg ( jblock(i_umt_nblock), * ndir, nshr, nstatev, nfieldv, nprops, lanneal, * stepTime, totalTime, dt, cmname, coordMp, charLength, * props, density, strainInc, relSpinInc, * tempOld, stretchOld, defgradOld, fieldOld, * stressOld, stateOld, enerInternOld, enerInelasOld, * tempNew, stretchNew, defgradNew, fieldNew, * stressNew, stateNew, enerInternNew, enerInelasNew, * jblock(i_umt_noel), jblock(i_umt_npt), * jblock(i_umt_layer), jblock(i_umt_kspt)) return end c subroutine vumatXtrArg ( c Read only - * nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal, * stepTime, totalTime, dt, cmname, coordMp, charLength, * props, density, strainInc, relSpinInc, * tempOld, stretchOld, defgradOld, fieldOld, * stressOld, stateOld, enerInternOld, enerInelasOld, * tempNew, stretchNew, defgradNew, fieldNew, c Write only - * stressNew, stateNew, enerInternNew, enerInelasNew, c Read only extra arguments - * nElement, nMatPoint, nLayer, nSecPoint ) include 'vaba_param.inc' c dimension props(nprops), density(nblock), coordMp(nblock,*), 1 charLength(nblock), strainInc(nblock,ndir+nshr), 2 relSpinInc(nblock,nshr), tempOld(nblock), 3 stretchOld(nblock,ndir+nshr), 4 defgradOld(nblock,ndir+nshr+nshr), 5 fieldOld(nblock,nfieldv), stressOld(nblock,ndir+nshr), 6 stateOld(nblock,nstatev), enerInternOld(nblock), 7 enerInelasOld(nblock), tempNew(nblock), 8 stretchNew(nblock,ndir+nshr), 9 defgradNew(nblock,ndir+nshr+nshr), 1 fieldNew(nblock,nfieldv), 2 stressNew(nblock,ndir+nshr), stateNew(nblock,nstatev), 3 enerInternNew(nblock), enerInelasNew(nblock) c c Documentation of extra arguments: c nElement: Array of internal element numbers dimension nElement(nblock) c nMatPoint: Integration point number c nLayer : Layer number for composite shells and layered solids c nSecPoint: Section point number within the current layer c character*80 cmname c do km = 1, nblock c c User coding ... c end do return end 为什么VMAT里面是jblock,不是nblock？尤其是这段 parameter ( * i_umt_nblock = 1, * i_umt_npt = 2, * i_umt_layer = 3, * i_umt_kspt = 4, * i_umt_noel = 5 ) call vumatXtrArg ( jblock(i_umt_nblock),... 不明白怎么就获取积分点的单元信息了。 |
4楼2013-02-12 22:18:37
ehe333333
木虫 (著名写手)
- 应助: 22 (小学生)
- 金币: 1804.8
- 散金: 1800
- 红花: 11
- 帖子: 1561
- 在线: 184小时
- 虫号: 2198330
- 注册: 2012-12-21
- 专业: 结构工程
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
ppzzyy: 金币+20, ★★★很有帮助 2013-03-22 03:28:31
感谢参与，应助指数 +1
ppzzyy: 金币+20, ★★★很有帮助 2013-03-22 03:28:31
|
找到一段子程序，应该就是来取得积分点编号的吧，6.4以上版本才适用： c subroutine vumatXtrArg ( c Read only - * nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal, * stepTime, totalTime, dt, cmname, coordMp, charLength, * props, density, strainInc, relSpinInc, * tempOld, stretchOld, defgradOld, fieldOld, * stressOld, stateOld, enerInternOld, enerInelasOld, * tempNew, stretchNew, defgradNew, fieldNew, c Write only - * stressNew, stateNew, enerInternNew, enerInelasNew, c Read only extra arguments - * nElement, nMatPoint, nLayer, nSecPoint ) include 'vaba_param.inc' |
2楼2013-02-12 18:35:32
ppzzyy
金虫 (正式写手)
- EPI: 1
- 应助: 5 (幼儿园)
- 金币: 1230.5
- 散金: 615
- 红花: 1
- 帖子: 316
- 在线: 142.4小时
- 虫号: 1041469
- 注册: 2010-06-14
- 性别: MM
- 专业: 结构工程
3楼2013-02-12 22:04:00
5楼2015-04-03 10:22:56
6楼2017-06-20 21:26:55
投票: