<<to gus1977 >>:你好！感谢你上次的指导！有新问题

jingyuanchun

*** WARNING ***                         CP=      4.031   TIME= 11:02:53
The program chosen initial timestep/load-factor is arbitrary.  It is   
necessary for the user to supply a suitable initial                     
timestep/load-factor through the NSUB or DELTIM command for            
convergence and overall efficiency.
计算过程中出现了上述警告，怎么处理！
再者，我想循环计算80次，请问怎么用程序实现这样的循环计算！
还有阿，
/post26
nsol,2,node(x,y,z),temp,,temp_time
plvar,2
得到节点(x,y,z)温度时间变化曲线,如何准确找到最高温度对应的时间点！
小妹谢谢版主

gus1977

哦，这个问题好解决，是我上次给你写代码的时候，没写全
*SET,tm_start,1e-16
*SET,tm_end,4.86
tm_inc=(tm_end-tm_start)/80
*do,i, 1, 80,
tm=tm_start+tm_inc*i
time,tm
deltim,tm_inc
NSEL,S,LOC,X,0.200,
sfdele,all,hflux
sf,all,hflux,519+(66-519)/4.86*tm
nsel,all
solve
*enddo 看看这回怎么样

jingyuanchun

我想分析的问题是这样子的,汽车循环制动80次,,包括制动和冷却过程,制动时间4.86s,冷却时间60s,得到制动80次后的温度场分布!
你给写的这个程序是制动程序,
anty,trans
outres,all,all
tunif,24  
*SET,tm_start,1e-16
*SET,tm_end,4.86
*SET,tm_inc,0. 01
*do,tm, tm_start, tm_end, tm_inc
time,tm
deltim,tm_inc
NSEL,S,LOC,X,0.200,
sf,all,hflux,519600+(66-519600)/4.86*tm
nsel,all
solve
*enddo
冷却阶段只有对流换热边界条件
NSEL,S,LOC,X,0.215
sf,all,conv,62.3,24
你看冷却阶段程序怎么写,
以及怎么实现制动冷却,循环80次????/
多谢!!!!!!

jingyuanchun

版主,很忙吗?讨论一下吧,呵呵,谢谢啊!

gus1977

呵呵，着急了，这两天不是周末吗，我休息。另外，请你注意，为了你问题的系统性，和避免重复性灌水，请不要重复发帖，你后来发的那个帖子我会给你删掉。
其实你的问题很简单，主要是变一下你的时间就好了

anty,trans
outres,all,all
tunif,24
*set,totaltm,0
*SET,tm_start,1e-16
*SET,tm_end,4.86
*SET,tm_inc,0. 01
*do,i,0,79
*do,tm, tm_inc, tm_end, tm_inc
totaltm=tm_end*i+tm
time,totaltm
deltim,tm_inc
NSEL,S,LOC,X,0.200,
sf,all,hflux,519600+(66-519600)/4.86*tm
nsel,all
solve
*enddo
*enddo

time,totaltm+nnnn
deltim,xxxx
NSEL,S,LOC,X,0.200,
sfdele,all,all
nsel,all
NSEL,S,LOC,X,0.215
sf,all,conv,62.3,24
solve

jingyuanchun

恩,我知道了,一着急就重复发帖子了,呵呵,以后不会了,我怕你看不到,我又比较着急做课题!!小妹谢谢版主,真的不胜感激!!!!!!

jingyuanchun

我想还是要有劳版主帮忙,版主上次的程序,我看过了;我表达的不清楚,所以让你费解了;我修改了一下,有错误,运行不久就不计算了!循环80次,不光是制动部分,还要冷却部分!不知道 这次你能不能明白我的意思,我也不知道怎么清楚的表达我的意思!汽车循环制动过程是这样子的:制动,冷却,再制动,冷却......如此循环下去,到某一时间停止.我附加 流程图给你,麻烦你了,小妹多谢!

/solu
anty,trans
outres,all,all
tunif,24
*set,totaltm,0
*SET,tm_start,1e-16
*SET,tm_end,4.86
*SET,tm_inc,0. 1
*do,i,0,79
*do,tm, tm_inc, tm_end, tm_inc
totaltm=tm_end*i+tm
time,totaltm
deltim,tm_inc
NSEL,S,LOC,X,0.200,
sf,all,hflux,519600+(66-519600)/4.86*tm
nsel,all
solve
*enddo

NSEL,S,LOC,X,0.200,
sfdele,all,all
nsel,all

time,totaltm+ 60*(i+1)
deltim,3
NSEL,S,LOC,X,0.215
sf,all,conv,62.3,24
solve
*enddo

gus1977

呵呵，你在考我啊，出现什么错误，是你修改之后出现错误吗，那没修改之前那

jingyuanchun

当然不是在考你了,我哪里有那个能力啊!我修改后的有问题!你看了我发的附件了吗?我觉得你写的程序 ,不是我想表达的意思,所以我修改了一下!却不成功!!!我的理解是这样的!我的qq号:37834822!版主,您有电话吗?我打给你,如果你方便的话!不方便就算了,我们就这样交流,只是麻烦您了!!多谢多谢!!!!

gus1977

/solu
anty,trans
outres,all,all
tunif,24
*set,totaltm,0
*SET,tm_start,1e-16
*SET,tm_end,4.86
*SET,tm_inc,0. 1
!cooltm 冷却时间
*do,i,0,79
NSEL,S,LOC,X,0.215
sfdele,all
nsel,all
*do,tm, tm_inc, tm_end, tm_inc
totaltm=(tm_end+cooltm)*i+tm
time,totaltm
deltim,tm_inc
NSEL,S,LOC,X,0.200,
sf,all,hflux,519600+(66-519600)/4.86*tm
nsel,all
solve
*enddo

NSEL,S,LOC,X,0.200,
sfdele,all,all
nsel,all

time,totaltm+ cooltm
deltim,3
NSEL,S,LOC,X,0.215
sf,all,conv,62.3,24
solve
*enddo

jingyuanchun

谢谢版主!(y)(y)

jingyuanchun

版主,计算过程中出现了很多警告,程序不再运行!警告如下:
** WARNING ***   SUPPRESSED MESSAGE    CP=     80.984   TIME= 10:51:17
Node 3428 on element 497 is unselected.                                 
*** WARNING ***   SUPPRESSED MESSAGE    CP=     80.984   TIME= 10:51:18
Node 4697 on element 499 is unselected.   
怎么解决啊?什么原因啊?
*** WARNING ***                         CP=      4.219   TIME= 10:31:40
Missing surface load label on the SFDE command.  The SFDE command is     ignored.   
再有劳版主!!!!多谢!!!!!

jingyuanchun

再者,cooltm应该设置60吧,
*set,totaltm,0
*SET,tm_start,1e-16
*SET,tm_end,4.86
*SET,tm_inc,0. 1

*SET, cooltm,60

!cooltm 冷却时间
*do,i,0,79
对吧??
为什么会出现这么多的警告,不运行 了呢?
多谢!

jingyuanchun

问题解决了，谢谢版主！！！原因在于忘记 nsel,all ;;;
在我做热结构耦合的时候，如有问题，还是得麻烦版主费神！！！！
谢谢版主，送你(g)