大体积混凝土温控讨论

jacobi

不知wang doc 的公式FLUX（1）=C*ρ*ΔT/Δt
是如何体现冷却的，
如果是使用冷却水管的话，应该是冷却水管通过的nodes处的节点考虑负热源，
但是单从wang doc 的公式看，
不知是如何体现的，

不知lxcad版主是如何理解的？
谢谢。

Lxcad

如无水管冷却，混凝土绝热温升为θ(t)=θ0(1-e-mt)，有水管冷却时，混凝土的温升为θ0ψ(t)，因此，对于混凝土绝热温升，水管冷却在单位时间内所吸收的热量为图示：

Lxcad

管冷讨论总结：
1,计算机性能允许，水管不多，可以考虑用FLUENT算温度场后导入到ABAQUS算应力；
2,用朱的等效方程考虑用FLUX建立负热源；
思路清晰后就可以往下了。

Wang_Jian

（1）我前面帖子里 FLUX（1）=C*ρ*ΔT/Δt 是HETVAL子程序中定义体积生热量的，是计算生热量的通用公式；给出的程序仅是水化热部分的，没有水管等效热传导部分。
（2）水管冷却的等效热传导中的水管吸收的热量是与当前混凝土温度有关的，因此是不断变化而计算前又无法估算的，无法用*DFLUX，*AMPLITUDE定义。
（3）用HETVAL实现等效水管冷却，需仔细研究朱伯芳先生的公式，包括前面的章节。原理与水化热的模拟相同，只是水化热的放热量可认为与当前温度无关，比较简单；而水管冷却与当前温度有关，要复杂一些。

[ 本帖最后由 Wang_Jian 于 2007-4-19 21:25 编辑 ]

jacobi

感谢wang doc以及lxcad版主的回复，
从两位的回复来看，这个管冷对我来说，确是难了些。

对于水化热，
朱伯芳先生的书上有以下放热公式：
Q（t）＝Q0×（1－exp（－mt））

以及他本人提出的公式

个人以为用生热函数公式计算会好些，
那个绝热公式和此接近却又不同，

请赐教，谢谢。

Lxcad

有新的计算方法如图：
但该论文没有说该计算阶段混凝土的平均温度怎样求！

[ 本帖最后由 LXCAD 于 2007-4-19 23:52 编辑 ]

Lxcad

平均温度只有设一变量记录，每次调用更新！！

[ 本帖最后由 LXCAD 于 2007-4-19 23:28 编辑 ]

jacobi

因为个人水平有限，
我再试试，
不知版主可否给出上述计算方法的出处，
如果是论文更好，
想弄明白些。
谢谢。

另外，不知版主有没有试试用两个不同函数（即绝热函数
和直接用那个放热函数）计算的出来的温度时程是否相同。
谢谢。

[ 本帖最后由 jacobi 于 2007-4-19 23:19 编辑 ]

Lxcad

肯定一样的，朱伯芳先生的书上2公式的系数都一样，只是温度*比热*质量=热量 了。算了，还是建议用22楼的公式，可以大体把握趋势的！

aresaran

LXCAD斑竹的照片真靓啊

Lxcad

原帖由 Wang_Jian 于 2007-4-16 11:19 发表
在子程序 HETVAL中：
FLUX（1）是生热率 【Heat flux, r (thermal energy per time per volume)】
FLUX（2）是FLUX(1)对时间的导数。对于混凝土温度计算没必要考虑FLUX（2）【FLUX（2）：Rate of change of h ...

我理解该段代码时有不同意见
        ag1 = TIME(1)
        ag2 = TIME(1) + DTIME
        IF (FID .GT. 0) THEN   
            TEMP1 = TH0 * (1.0D0 -DEXP(-bm*ag1))
            TEMP2 = TH0 * (1.0D0 -DEXP(-bm*ag2))
        ELSE
            TEMP1 = TH0 * ag1/(ag1+bm)
            TEMP2 = TH0 * ag2/(ag2+bm)
        ENDIF

我认为不考虑分层分Step模拟
         ag1 = TIME(1)-DTIME
         ag2 = TIME(1)
因为在hetval中对TIME定义如下：
TIME(1)      Step time at the end of the increment.该增量步结束时的时间，不是增量步开始时时间



TIME(2)    Total time at the end of the increment 总的时间

是否wang博的time(2)不是上面含义！
谢谢wang博光临！

doutsh

LXCAD版主在大体积混凝土温控方面展开深层次的讨论，并且引起各位大侠的共鸣。温控分析确实有一定的难度，同时还需要了解大体积混凝土施工的工程背景，这样的讨论肯定会使ABAQUS在温控方面发挥越来越大的作用。

doutsh

LXCAD版主在大体积混凝土温控方面展开深层次的讨论，并且引起各位大侠的共鸣。温控分析确实有一定的难度，同时还需要了解大体积混凝土施工的工程背景，这样的讨论肯定会使ABAQUS在温控方面发挥越来越大的作用。

jacobi

原帖由 LXCAD 于 2007-5-5 17:46 发表


我理解该段代码时有不同意见
        ag1 = TIME(1)
        ag2 = TIME(1) + DTIME
        IF (FID .GT. 0) THEN   
            TEMP1 = TH0 * (1.0D0 -DEXP(-bm*ag1))
            TEMP2 = TH0 *  ...

估计wang doc的分析中仅一个step，这种情况下效果是一样的。
呵呵。

Lxcad

原帖由 jacobi 于 2007-5-6 21:09 发表




估计wang doc的分析中仅一个step，这种情况下效果是一样的。
呵呵。

我查了，abaqus中time(1)是该增量步结束时在该step中的时间。且在msdev中调试成功

wang doc采用该法也未尝不可，他相当与把时间往后推了，而增量步设很小的情况下，计算结果也差不多的！

jacobi

如果有多个steps时，就必须要用time（2）了，
一般需要平移时间，
分层浇筑时的单元生死配合bc变化在琢磨中。。。

Wang_Jian

我是这么实现分层浇筑模拟的。注意： Age = TIME(2) -CastT -DTIME/2.0
这是真正的龄期。这个问题前面讨论了很多。


      SUBROUTINE HETVAL(CMNAME,TEMP,TIME,DTIME,STATEV,FLUX,
     1 PREDEF,DPRED)

            ……

C     mat10201  材料名称，包含浇筑区、层、材料号。  注意材料命名规则
        ZoneNoStr= CMNAME(4:4)                                                      
        LayNoStr = CMNAME(5:6)
        MatNoStr = CMNAME(7:8)

c       字符型转化为数值
        READ (ZoneNoStr*) ZoneNo   !取区号
        READ (LayNoStr*) LayNo       !取层号
        READ (MatNoStr*) MatNo      !取材料号

  c    调用子程序GetCastTime，根据浇筑区、浇筑层和材料编号确定浇筑时间，用于确定龄期
        CastT = GetCastTime(ZoneNo,LayNo,MatNo)

        Age = TIME(2) -CastT -DTIME/2.0
                ……

        FLUX(1) = C(i)*DEN(i)*TH0(i)*b(i)/(AGE+b(i))/(AGE+b(i))  !双曲线型

        ……

[ 本帖最后由 Wang_Jian 于 2007-5-21 19:36 编辑 ]

Lxcad

欢迎wang 继续来到

[ 本帖最后由 LXCAD 于 2007-5-27 09:03 编辑 ]

Lxcad

以上我就一基础上分3次浇筑的温度场云图，其中考虑管冷和用SFILM子程序考虑外界对流，唯一没有考虑是交界面上节点温度的处理，有文献建议加一负热源，wang博和其他高手如J能发表高论吗！

wxm19820202

请问
ZoneNoStr= CMNAME(4:4)                                                   
  LayNoStr = CMNAME(5:6)|
  MatNoStr = CMNAME(7:8)
中（4：4）（5：6）（7:8）代表什么意思呢？