蠕变、重启动、重新设置边界、在计算过程中改变参数

水幽寒 发表于 2006-4-1 14:02:41 | 显示全部楼层 | 阅读模式
这是我编的一段命令流,还很不完善,其中涉及蠕变、重启动、重新设置边界、在计算过程中改变参数等问题。
fini
/cle
/fil,exam0610
/COM, Structural
/config,NRES,6000
/prep7
et,1,beam23
r,1,10*9,10*9**3/12,9,,,,
mp,prxy,1,0.3
mp,dens,1,7.8e3
mp,ex,1
*dim,day1,array,50!定义时间数组
*dim,eex1,array,50!定义弹模数组
*dim,epcr1,array,50
*dim,epel1,array,50
*do,i,1,50,1
day1(i)=2*i!时间
eex1(i)=50*i!随着时间而增长的弹性模量
*enddo
eex1(50)=10000000!永久弹性模量
k,1
k,2,10
l,1,2
lsel,all
latt,1,1,1
esize,1
lmesh,all
f,2,fx,-22932000
nsel,s,loc,x,0
d,all,all
*do,i,1,50,1
time,day1(i)!时间
/soluokok.org
RESCONTRL,DEFINE,ALL,-1,1
eqslv,sparse
NLGEON,ON
NROPT,FULL
bfunif,temp,300 !温度
*if,i,gt,2,then
*if,abs(epcr1(i)/epel1(i)-1),gt,0.001,then
eex1(i)=eex1(50)
ANTYPE,,REST,
PARRES, CHANGE ,param, TXT,
*endif
*endif
uimp,1,ex,,,eex1(i)
tb,creep,1!使用显式蠕变方程
tbdata,1,0.0032,1,,0,,1
solve
finish
/POST26
ESOL,2,1,,LEPCR,1,EPCR ! STORE CREEP STRAIN
PRVAR,2 ! PRINT STRAIN VARIATION WITH TIME
esol,3,1,,lepel,1,epel
prvar,3
*get,epcr1(i+1),VARI,2,rtime,day1(i)
*get,epel1(i+1),VARI,3,rtime,day1(i)
parsav,all,param,txt
*enddo
在各个不同的时期改变试件的弹性摸量,但是到一定的条件下则采用永久的摸量。
我想在计算程序中加入一段命令流,提取结果中的蠕变和弹变值。
如果(蠕变值)/(弹变值)-1大于0.0001的话,则调用eex1(50)(永久弹性模量)。
我想在在下一步计算中判断后,调用上一步的边界条件,老是出问题,好象提示除数是0,另外还有很多的问题,大家可以讨论下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

 
 
  • QQ:56984982
  • 点击这里给我发消息
    电话:13527553862
    站务咨询群桥头堡站务咨询桥梁专业交流群:
    中国桥梁专业领袖群
    工作时间
    8:00-18:00