- 积分
- 41
- 注册时间
- 2003-5-4
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2009-4-25 21:40:20
|
显示全部楼层
来自 日本
本帖最后由 tozhou 于 2009-4-25 21:50 编辑
我还没解决啊~~~~
cypemail 发表于 2009-4-25 17:08
两种方法解决:
1, 后处理时在 tools-creat field output- -s.press/s.mises
适合于只是某几步用。
2, 把下面这段贴过去建立一个.for子程序。适合于每步都需输出用。
而且别忘了在材料特性里把用户自定义输出变量给选上。
************************************************************************
* The Subroutine for Stress Triaxiality
* for ABAQUS/Standard
* by tozhou@simwe
* 2009/4/25
************************************************************************
subroutine uvarm(uvar,direct,t,time,dtime,cmname,orname,
1 nuvarm,noel,npt,layer,kspt,kstep,kinc,ndi,nshr,
2 coord,jmac,jmatyp,matlayo,laccfla)
include 'aba_param.inc'
character*80 cmname,orname
character*3 flgray(15)
dimension uvar(nuvarm),direct(3,3),t(3,3),time(2)
dimension array(15),jarray(15),jmac(*),jmatyp(*),coord(*)
call getvrm('SINV',array,jarray,flgray,jrcd,jmac,jmatyp,
1 matlayo,laccfla)
***** Average Stress
***** -array(3)
***** Equivalent Stress
***** array(1)
***** Stress Triaxiality
if(array(1).eq.0.0) then
uvar(1)=0.0
else
uvar(1)=-array(3)/array(1)
endif
return
end |
评分
-
1
查看全部评分
-
|