本帖最后由 xzyzqs 于 2013-12-27 23:57 编辑
我在做桩基换热,考虑间歇运行(运行12小时,停止12小时(热流为0))的数值模拟,就是给桩内换热管内壁一个热流密度,然后模拟桩及围岩周围的温度场时空特征。
我的做法是:由于施加了初始均匀温度场,所以第一天为正常加载,以后的时间采用*do循环求解。取了几个特征点,了解其温度随时间的变化,但是循环部分没有显示曲线,如下图(为简化,计算了两天dy=2):
inp命令流如下:
!加热工况,运行12小时停12个小时,运行一个月
FINISH !退出当前处理器
/CLEAR,START !重新开始一个新的分析
/FILNAME,heating_only !定义工作名
*SET,dy,3 !计算的天数
/prep7 !进入前处理器
ET,1,PLANE55 !定义单元类型P185
!定义换热管热物参数
MP,KXX,1,0.42 !定义pe管的导热系数P180
MP,C,1,1465 !定义pe管的比热容
MP,DENS,1,1100 !定义pe管的密度
!定义桩热物参数
MP,KXX,2,2.1 !定义桩的导热系数
MP,C,2,800 !定义桩材料的比热容
MP,DENS,2,2500 !定义桩的密度
!定义围岩的热物参数
MP,KXX,3,1.8 !定义围岩的导热系数
MP,C,3,900 !定义围岩的比热容
MP,DENS,3,1800 !定义围岩的密度
!定义几何模型
cYL4,0,0,0,0,0.5,360,0 !桩,半径0.5m,p132
cYL4,0,0,0.5,0,4,360,0 !围岩,半径4m
cyl4,0.05,0.1625,0,0,0.0125,360,0 !1#换热管
cyl4,-0.05,0.1625,0,0,0.0125,360,0 !2#换热管
wpstyl,,,,,,,,, !显示工作平面P443
wprota,120,0,0 !将工作平面绕Z轴旋转120°(逆时针旋转为正)
cyl4,0.05,0.1625,0,0,0.0125,360,0 !3#换热管
cyl4,-0.05,0.1625,0,0,0.0125,360,0 !4#换热管
wprota,120,0,0 !将工作平面绕Z轴旋转120°(逆时针旋转为正)
cyl4,0.05,0.1625,0,0,0.0125,360,0 !5#换热管
cyl4,-0.05,0.1625,0,0,0.0125,360,0 !6#换热管
/pnum,area,1 !显示面的编号p422
/pnum,line,1
asba,1,3
asba,9,4
asba,1,5
asba,3,6
asba,1,7
asba,3,8
wprota,120,0,0 !将工作平面绕Z轴旋转120°(逆时针旋转为正)
!下面创建换热管
cyl4,0.05,0.1625,0.0102,0,0.0125,360,0 !1#换热管
cyl4,-0.05,0.1625,0.0102,0,0.0125,360,0 !2#换热管
wpstyl !显示工作平面P443
wprota,120,0,0 !将工作平面绕Z轴旋转120°(逆时针旋转为正)
cyl4,0.05,0.1625,0.0102,0,0.0125,360,0 !3#换热管
cyl4,-0.05,0.1625,0.0102,0,0.0125,360,0 !4#换热管
wprota,120,0,0 !将工作平面绕Z轴旋转120°(逆时针旋转为正)
cyl4,0.05,0.1625,0.0102,0,0.0125,360,0 !5#换热管
cyl4,-0.05,0.1625,0.0102,0,0.0125,360,0 !6#换热管
aglue,all !面粘合操作P154
NUMCMP,all !压缩所有的编号P219
wpstyl !关闭工作平面
!下面划分单元
!划分换热管的单元
mshape,1,2d !指定单元的类型(2d,三角形)p198
type,1 !设置换热管单元的类型p183
mat,1 !设置换热管单元的材料p183
LSEL,S,,,9,56,1 !选择换热管上的线P402
LESIZE,ALL,,,14 !对所选择的线设置网格单元个数p191
amesh,3,8,1 !对面3-8进行网格划分p194
!划分桩体单元
mshape,1,2d
type,1
mat,2
lsel,s,,,1,4,1
LESIZE,ALL,,,20
allsel
amesh,1
!划分围岩单元
meshape,1,2d
type,1
mat,3
lsel,s,,,5,8,1
LESIZE,ALL,,,20
amesh,2
finish !Exits normally from a processor
/SOLU !进入求解器
ANTYPE,TRANS !指定分析类型p235
TIMINT,ON,THERM !打开瞬态效应P245
TIME,43200 !定义计算终止时间,半天p241
deltim,60,60,600 !定义初始荷载步p242
AUTOTS,ON !打开自动时间步长p240
OUTRES,,ALL !写入每个子步
bfunif,temp,19 !施加均匀温度荷载p287
KBC,1 !指定载荷步阶跃方式p242
LSEL,S,,,9,56,1
cm,xzy,line
sf,xzy,hflux,156.1 !施加热流荷载p273
ALLSEL
SOLVE
FINISH
/SOLU !进入求解器
ANTYPE,,rest !设定分析类型为重启动分析,重启动点为上次求解结束时
TIMINT,ON,THERM !打开瞬态效应P245
TIME,86400 !定义计算终止时间,一天p241
deltim,300,300,600 !定义初始荷载步p242
AUTOTS,ON !打开自动时间步长p240
OUTRES,,ALL !写入每个子步
KBC,1 !p242
LSEL,S,,,9,56,1
cm,xzy,line
sf,xzy,hflux,0
ALLSEL
SOLVE
finish
*DO,i,2,dy,1
time1=86400*i-43200
time2=86400*i
/SOLU
ANTYPE,,rest
TIMINT,ON,THERM
TIME,time1
deltim,300,300,600
AUTOTS,ON
OUTRES,,ALL
KBC,1
LSEL,S,,,9,56,1
cm,xzy,line
sf,xzy,hflux,156.1
ALLSEL
SOLVE
finish
/SOLU
ANTYPE,,rest
TIMINT,ON,THERM
TIME,time2
deltim,300,300,600
AUTOTS,ON
OUTRES,,ALL
KBC,1
LSEL,S,,,9,56,1
cm,xzy,line
sf,xzy,hflux,0
ALLSEL
SOLVE
finish
*ENDDO
/post26
/axlab,x,time,(sec) !指定X轴标题p433
/axlab,y,temperture,(℃) !指定y轴标题
/gthk,axis,1 !指定坐标轴的粗度
/gthk,curve,1 !指定曲线粗度
/color,curve,mred,1 !设置曲线显示颜色
/replot
/xrange,0,86400*2 !指定x坐标轴显示范围
/yrange,18,40 !指定y坐标轴显示范围
nsol,2,112,temp,,换热管内壁的温度 !定义变量2 p372
nsol,3,15,temp,,换热管外壁的温度
nsol,4,1031,temp,,桩土接触面的温度
/grid,0 !在图形显示时选择网格的类型p432
plvar,2,3,4 !绘制不同位置节点温度随时间的变化关系曲线p379
/IMAGE,SAVE,temp_2d2,BMP
finish
|