先对问题说明如下:
偶要用温度场模拟渗流场,做了一个四层的问题,第三层有一个点源模拟抽水或注水,想要看它每个时间步的变化情况,但每次都是只算了一次,且每次都出现警告。偶把程序贴上,哪位大虾帮忙看看问题出在哪里咯,叩谢了!
/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. |