声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1819|回复: 2

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

[复制链接]
发表于 2006-4-14 11:21 | 显示全部楼层 |阅读模式

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

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

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

使用道具 举报

发表于 2006-4-16 07:29 | 显示全部楼层

回复:(cxg1223)[求助]混凝土水化热施加上面出了问题...

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

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-15 05:07 , Processed in 0.061580 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表