用vusdfld子程序实现材料属性随纵坐标变化
本帖最后由 taozi136161055 于 2013-5-18 21:02 编辑最近做个实例，需要实现材料随着纵坐标变化，于是考虑子程序vusdfld，学了一段时间的Fortran，着手编了一个小程序，经验证可以运行，代码如下。
*******************************************
c
c User subroutine VUSDFLD for user-defined fields
c
subroutine vusdfld(
c Read only -
* nblock, nstatev, nfieldv, nprops, ndir, nshr,
* jElemUid, kIntPt, kLayer, kSecPt,
* stepTime, totalTime, dt, cmname,
* coordMp, direct, T, charLength, props,
* stateOld,
c Write only -
* stateNew, field )
c
include 'vaba_param.inc'
c
dimension props(nprops),
* jElemUid(nblock), coordMp(nblock, *),
* direct(nblock, 3, 3), T(nblock,3,3),
* charLength(nblock),
* stateOld(nblock, nstatev),
* stateNew(nblock, nstatev),
* field(nblock, nfieldv)
character*80 cmname
parameter( nrData=6 )
c
character*3 cData(maxblk*6)
dimension jData(maxblk*nrData)
dimension eqps(maxblk*6),stress(maxblk*6)
c
if(stepTime.eq.0) then
do k = 1, nblock
y = coordMp(k,2)
field(k,1) = 0 - y
c
end do
end if
c
return
end
************************************分割线
希望对初学vusdfld的大家有所帮助。
关于abaqus、vusdfld与usdfld的问题可以多多交流，共同提高！！
楼主这个子程序怎么调用和USDFLD一样吗，在property中选中‘user defined field’?? 好贴正需要顶起 好帖子，希望顶起来，作者写的更详细些
有没有inp文件啊！急需啊！ props(pprops)是user matiaral中的那些吗 楼主还在吗,如果在usdfld里面怎么用呢 coordmp提出的是物质点的坐标值么？为什么我提取出来的全是零啊 感谢楼主分享的资料，对我很有用。我在想是不是初始时刻判定的时候可以改为 if stepTime.le.dt，相当于第一个增量步的时候，开始赋予材料属性。因为之前用 stepTime.eq.0，发现并没有赋值成功。 本帖最后由 jacklee_94 于 2022-1-4 13:07 编辑
18354223911 发表于 2021-9-15 10:04
感谢楼主分享的资料，对我很有用。我在想是不是初始时刻判定的时候可以改为 if stepTime.le.dt，相当于第一 ...
"steptime.le.dt"这个语句赋值成功 西兰花佛系读博 发表于 2018-5-4 21:39
coordmp提出的是物质点的坐标值么？为什么我提取出来的全是零啊
同问，这个coordmp提取出来的是积分点（单元中心点）的坐标，还是节点的坐标？
页:
[1]