[求助]混凝土水化热施加上面出了问题

cxg1223 · 发表于 2006-4-14 11:21

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

我最近在混凝土水化热施加上面出了问题了：我用绝热温升公式＝24.9＊（1－exp(-0.351*time))导出水化热公式。 但是施加时如用函数： HE00=59867*(exp(-0.351*(time-1))-exp(-0.351*time))施加和用table施加： *DIM,B,TABLE,10,1,1,TIME,TEMP, !* BFV,P51X,HGEN, %B% *SET,B(1,0,1) , 1 *SET,B(1,1,1) , 12021.72 *SET,B(2,0,1) , 2 *SET,B(2,1,1) , 18032.58 *SET,B(3,0,1) , 3 *SET,B(3,1,1) , 7213.032 *SET,B(4,0,1) , 4 *SET,B(4,1,1) , 6010.86 *SET,B(5,0,1) , 5 *SET,B(5,1,1) , 4808.688 *SET,B(6,0,1) , 6 *SET,B(6,1,1) , 3606.516 *SET,B(7,0,1) , 7 *SET,B(7,1,1) , 2404.344 *SET,B(8,0,1) , 8 *SET,B(8,1,1) , 2404.344 *SET,B(9,0,1) , 9 *SET,B(9,1,1) , 2163.91 *SET,B(10,0,1) , 10 *SET,B(10,1,1) , 1202.172 模拟分析（就是设体初温为0，无对流）结果大不一样，函数的结果是14度多，而table的是28度。这与预想的24.9度都不吻合，这是为什么？是我那里出错了？ [em06]

AaronSpark · 发表于 2006-4-16 07:29

两个加载能对应上吗，没仔细看 下面是一个例子 !计算温度模块 fini /clea !!!!!!!!!!!!!! *set, day, 1 !计算天数，先调试1天 *set, hour, day*2 !计算小时数 *set, BEG_TEMP, 24 !初始温度 *set, air_max, 28 !最高气温 *set, air_min, 22 !最低气温 *set, l, 15 !浇注墙体长度的1/2 *set, h, 3.8 !浇注墙体的高度 *set, t, 0.5 !浇注墙体的厚度 *set, plank, 0.8 !底板厚 *set, board, 2 !底板宽计算长度 !!!!!!!!!!!!!!!!!!!!! /prep7 et,1,solid70 mp,dens,1,2400 mp,kxx,1,8.595 !混凝土的导热系数 mp,c,1,0.948 !混凝土的比热容 !墙体、底板、模板建模 block,,l,0.4,h,,-t !墙体X长L-0，Y高H-0.4，Z厚-T block,,l,0.4,0,,-t !高0.4墙体X长L-0，Y高-0.4， Z厚-T block,,l,-plank,0,0,-t !墙体下的底板X长L-0，Y高PLANK，Z厚-T block,,l,-plank,0,-t,-board !X长L-0，Y高PLANK，Z厚T-BOARD !把已建立的体连接起来，同时创建模板 vglue, 1,2 vglue, 3,5 vglue, 3,4 block,,l,0.4,h,,0.02 vglue, 1,4 block,,l,0.4,h,-t,-t-0.02 vglue, 1,4 block,l,l+0.02,0.4,h,0,-t !建立模板 vglue, 1,4 !墙体及底板分网,宽度方向分5份;长度方向分成35*3/2份;高度分16份 lsel,s,length,,t,t !选择长度为t的线 lesize,all,,,5 !划分为5等分 lsel,s,length,,l,l !选择长度为l的线 lesize,all,,,35*3/2 !划分为35*3/2等分 lsel,s,length,,h-0.4,h-0.4 !选择长度为h-0.4的线 lesize,all,,,16 !划分为16等分 lsel,s,length,,0.4,0.4 !选择长度为0.4的线 lesize,all,,,2 !划分为2等分 lsel,s,length,,plank,plank !选择长度为plank的线 lesize,all,,,2 !划分为2等分 lsel,s,length,,board-t,board-t !选择长度为board-t的线 lesize,all,,,5 !划分为5等分 vmesh,1 vmesh,2 vmesh,3 vmesh,5 !模板材料 mp,dens,2,2500 mp,kxx,2,0.837 !模板的导热系数 mp,c,2,1.8 !模板的比热容 mat, 2 !调用材料号为2的材料 type,1 !模板划分网格,大小与墙体对应 lsel,s,length,,l,l !选择长度为l的线 lesize,all,,,35*3/2 !划分为35*3/2等分 lsel,s,length,,0.02,0.02 !选择长度为0.02的线 lesize,all,,,1 !划分为1等分 lsel,s,length,,h-0.4,h-0.4 !选择长度为h-0.4的线 lesize,all,,,16 !划分为16等分 vmesh,8 vmesh,6 vmesh,7 !选中模板表面,建立WOOD(木头)组合面 asel,s,,,20,21 asel,a,,,37 cm,wood,area !选中的面组成1组,编号为wood !选中混凝土表面,建立conc1,conc2组合面 allsel asel,s,,,4 asel,a,,,8 asel,a,,,14 asel,a,,,11 asel,a,,,18 asel,a,,,25 asel,a,,,7 asel,a,,,9 asel,a,,,12 cm,conc1,area !选中的面组成1组,编号为conc1 allsel asel,s,,,1,2 asel,a,,,6 cm,conc2,area !得到墙体的节点数 allsel vsel,s,,,1,2 !选中编号为1，2的两个体 nslv,s,1 !选中体上的节点 !建立对称面约束 allsel asel,s,loc,x,0 !建立X坐标为0的面 da,all,symm DTRAN !生成GEN_HOT文件，保存每小时混凝土的生热量 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *dim, he, array,1 *cfopen,GEN_HOT,txt HE00=330*536*(1-exp(-0.69/24)) *vwrite,HE00 (f15.5) *do, i, 1, hour HE00=330*536*(1-exp(-0.69*(i/24)**0.56)) HE01=330*536*(1-exp(-0.69*(((i+1)/24))**0.56)) *vwrite,HE01-HE00 (f15.5) *enddo *cfclos !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !求解 /SOLU ANTYPE,TRANSIENT,NEW !定义分析类型 nropt, full allsel, all ! nsel,s,,,all !选所有节点 ic,all,temp,25 !初始温度25 *do, i, 1, hour time,i timint, on *vread, he(1),GEN_HOT,txt,,,,,,i-1 (f25.15) vsel,s,,,1 eslv,s BFe,ALL,HGEN,,he(1) !加水化热 *if,i,ge,48,then !时间超过48小时要拆模 vsel,s,,,6,8 eslv,s ekill,all cmsel,s,conc1 cmsel,a,conc2 !在混凝土表面施加对流系数和空气温度 sfa,all,,CONV,53,air_max/2+air_min/2+(air_max/2-air_min/2)*cos(15*(i-7)) *else cmsel,s,wood !在模板表面施加对流系数和空气温度 sfa,all,,CONV,49.40,air_max/2+air_min/2+(air_max/2-air_min/2)*cos(15*(i-7)) cmsel,s,conc1 !在混凝土表面施加对流系数和空气温度 sfa,all,,CONV,53,air_max/2+air_min/2+(air_max/2-air_min/2)*cos(15*(i-7)) *endif allsel TRNOPT,FULL deltim, 0.5,0.1,1 ! AUTOTS,ON KBC,0 outres,all,1 solve *enddo allsel !后处理 /post1 allsel !全选 vsel,s,,,1,2 nslv,s,1 !选择所选体的节点 *get,n_num,node,,count !得到的节点数 *cfopen,temp,txt,,append !创建文件 *do, i, 1, hour set, I *do, ii, 1, n_num *get,te,node,ii,temp *vwrite,te (f15.10) *enddo *enddo *cfclos !文件关闭 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *dim, tt1, array,n_num *cfopen,t_temp,txt,,append ! *vread,tt1(1),temp,txt (f15.10) *do,i,1,n_num *set,s,tt1(i) s=s-BEG_TEMP !得到温差 *vwrite,s (f15.10) *enddo *cfclos !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *dim, tt2, array,n_num *dim, tt3, array,n_num *dim, tt4, array,n_num *dim, tt5, array,n_num *cfopen,t_temp,txt,,append ! !创建t_temp.txt文件 *do,i,1,hour-1 *vread,tt2(1),temp,txt,,,,,,n_num*(i-1) !读出数据 (f15.10) *vread,tt3(1),temp,txt,,,,,,n_num*I !读出数据 (f15.10) *voper,tt4(1),tt3(1),sub,tt2(1) !数组相减 *vwrite,tt4(1) (f15.10) *enddo *cfclos save /post26 !进入时间历程后处理 nsol,2,2286,temp,,nodetemp !变量2286处的温度放入变量2中 plvar,2 !显示变量的图形 /COLOR,PBAK,OFF /replot /ui,copy,save fini

cxg1223 · 发表于 2006-4-16 09:45

谢谢回答，这个例题我看过了。我目前关注的是这两种加载方式应该没问题，抛去其他因素，验证但是没有一个得出正确结果的。我的Q(TIME)=24.9*(1-EXP(-0.351*TIME))*0.996(比热)*2400(密度）

账号		自动登录	找回密码
密码			我要加入

[综合讨论] [求助]混凝土水化热施加上面出了问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（cxg1223）[求助]混凝土水化热施加上面出了问题...