lxg613 发表于 2009-9-12 16:42:45

求助:关于时间步控制的问题

先对问题说明如下:
偶要用温度场模拟渗流场,做了一个四层的问题,第三层有一个点源模拟抽水或注水,想要看它每个时间步的变化情况,但每次都是只算了一次,且每次都出现警告。偶把程序贴上,哪位大虾帮忙看看问题出在哪里咯,叩谢了!
/PREP7
/TITLE, deformation in thermal load(four layers)
nrow=41!行的点数
nlie=41!列的点数
xdiv=0.25   !x-方向的间距
ydiv=0.25   !y-方向的间距
!建两个点,插入点
N,1,0,0,0,,,,   
N,nrow,(nrow-1)*xdiv,0,0,,,,
!*
FILL,1,nrow,nrow-2, , ,1,1,1,
!进行COPY
nsel,all
NGEN,nlie,nrow,all, , ,0,ydiv,0,1,
ET,1,PLANE13,4,,2
!材料1属性
MP,EX,1,5.0e9             !弹性模量
MP,ALPX,1,-0.1627E-10*98000*44.6   !热膨胀系数
MP,KXX,1,5.92e-5/0.743          !热传导率
MP,DENS,1,2647            !密度
MP,C,1,1.496E-17          !比热容
MP,PRXY, 1,0.25
!材料2属性
MP,EX,2,2.0e9             !弹性模量
MP,ALPX,2,-1.1627E-10*98000*44.6   !热膨胀系数
MP,KXX,2,1.92e-6/0.743         !热传导率
MP,DENS,2,2647            !密度
MP,C,2,1.496E-17          !比热容
MP,PRXY, 2,0.25
!four layers 2:1:1:1,mat:2,1,2,1
*do,j,1,0.4*(nlie-1),1
*do,i,1,nrow-1,1
type,1
mat,2
e,(j-1)*nrow+i,(j-1)*nrow+i+1,j*nrow+i+1,j*nrow+i
*enddo
*enddo
*do,j,0.4*(nlie-1)+1,0.6*(nlie-1),1
*do,i,1,nrow-1,1
type,1
mat,1
e,(j-1)*nrow+i,(j-1)*nrow+i+1,j*nrow+i+1,j*nrow+i
*enddo
*enddo
*do,j,0.6*(nlie-1)+1,0.7*(nlie-1),1
*do,i,1,nrow-1,1
type,1
mat,2
e,(j-1)*nrow+i,(j-1)*nrow+i+1,j*nrow+i+1,j*nrow+i
*enddo
*enddo
*do,j,0.7*(nlie-1)+1,nlie-1,1
*do,i,1,nrow-1,1
type,1
mat,1
e,(j-1)*nrow+i,(j-1)*nrow+i+1,j*nrow+i+1,j*nrow+i
*enddo
*enddo
!边界条件和初始温度
NSEL,S,LOC,Y,0
D,ALL,UY
D,ALL,UX
D,ALL,TEMP,0
NSEL,S,LOC,X,0
D,ALL,UY
D,ALL,UX
D,ALL,TEMP,0
NSEL,S,LOC,X,(nrow-1)*xdiv
D,ALL,UY
D,ALL,UX
D,ALL,TEMP,0
NSEL,S,LOC,Y,(nlie-1)*ydiv
D,ALL,UY
D,ALL,UX
D,ALL,TEMP,0
NSEL,S,LOC,Y,5
NSEL,R,LOC,X,5
BF,ALL,HGEN,3.125e-1
NSEL,ALL
FINISH
!求解
/SOLU
ANTYPE,TRANS       !瞬态分析
TIME,1
TIMINT,ON          !打开时间积分
DELTIM,0.05,0.1,0.2
OUTRES,BASIC,ALL
TUNIF,0         !initial uniform temperature
SOLVE
!*** WARNING ***
! The maximum time step (DTMAX = 0.1) must be greater than the starting   
! time step (DTIME = 1).                                                
!The DELTIM command is ignored.

lxg613 发表于 2009-9-12 16:56:55

大虾们,帮忙看看咯

lxg613 发表于 2009-9-13 18:12:28

解决了
把后面那句DELTIM,0.05,0.1,0.2
改成DELTIM,0.05就可以了

呵呵呵,给自己鼓励一下

sigh08 发表于 2009-9-14 15:05:16

3# lxg613
学习了,:)
页: [1]
查看完整版本: 求助:关于时间步控制的问题