ANSYS 热分析中外界环境温度随时间变化如何求解 - 仿真模拟 - ANSYS & ABAQUS

21楼2016-03-22 19:16:53

mygt_hit

专家顾问 (职业作家)

专家经验: +362
应助: 438 (硕士)
贵宾: 0.019
金币: 19571.1
散金: 5130
红花: 135
沙发: 3
帖子: 4938
在线: 990.5小时
虫号: 1489764
注册: 2011-11-13
性别: GG
专业: 结构工程
管辖: 土木建筑

引用回帖:

21楼: Originally posted by 德群 at 2016-03-22 19:16:53
我就把时间设为time 180，得到的结果温度还是稳定在那个值，大神，是不是我删除有问题，他前面提示我一个warning
*** WARNING ***
When using the DDELE command with a first order syst ...

这样试一下
把温度还是定义为ARRAY （不是TABLE），将BF命令改为D命令，类似于D,ALL,TEMP,T_WENDU(k)，看看结果是不是你预期的。

SOLID70单元的节点自由度是温度，节点上可加的体荷载BF是Heat Generations 生热率，请看单元帮助。
另外可以参考 http://wenku.baidu.com/link?url= ... 18wfa0zr0_prXX5K###

知其然，知其所以然。

22楼2016-03-22 21:31:28

德群

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 60.3
红花: 1
帖子: 74
在线: 254小时
虫号: 3524726
注册: 2014-11-07
专业: 结构工程

引用回帖:

22楼: Originally posted by mygt_hit at 2016-03-22 21:31:28
这样试一下
把温度还是定义为ARRAY （不是TABLE），将BF命令改为D命令，类似于D,ALL,TEMP,T_WENDU(k)，看看结果是不是你预期的。

SOLID70单元的节点自由度是温度，节点上可加的体荷载BF是Heat Generations 生热 ...

大神，您说的是这个意思么？
FINISH
/FILNAME,YICIXUNHUAN
/CLEAR
/PREP7
ET,1,SOLID70
MP,KXX,1,1.8
MP,KYY,1,1.8
MP,KZZ,1,1.8
MP,DENS,1,2400
MP,C,1,950
MP,KXX,2,0.61
MP,KYY,2,0.61
MP,KZZ,2,0.61
MP,DENS,2,996
MP,C,2,4185
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
BLOCK,0,200E-3,0,200E-3,0,600E-3
WPAVE,0,0,400E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

WPROTA,0,-90,0
WPAVE,0,0,50E-3
CSYS,4
VSBW,ALL
WPAVE,0,0,100E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

WPROTA,0,0,90
WPAVE,0,0,50E-3
CSYS,4
VSBW,ALL
WPAVE,0,0,100E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

ALLSEL,ALL
VSEL,S,LOC,Z,0,400E-3
VSEL,R,LOC,X,50E-3,150E-3
VSEL,R,LOC,Y,50E-3,150E-3
VATT,1,,1
ALLSEL,ALL
VSEL,U,MAT,,1,1,1
VATT,2,,1

ALLSEL,ALL
lesize,all,,,10
ALLSEL,ALL
MSHKEY,1
VMESH,ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/SOLU
ALLSEL,ALL
VSEL,S,MAT,,1
NSLV,S,1
D,ALL,TEMP,20                                              !ic不能用？
ALLSEL,ALL
VSEL,S,MAT,,2
NSLV,S,1
D,ALL,TEMP,15
ALLSEL,ALL
ANTYPE,TRANS
TIMINT,OFF
TIME,0.01
DELTIM,0.01
SOLVE
DDELE,ALL,TEMP

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*DIM,T_WENDU,ARRAY,240,1
*DO,k,1,240,1
*SET,T_WENDU(k),10*SIN(0.0262*k)
*ENDDO

*DO,k,1,240,1
TIME,k
TIMINT,ON
AUTOTS,ON
KBC,1
DELTIM,13,2,50
ALLSEL,ALL
DDELE,ALL,TEMP
ALLSEL,ALL
VSEL,S,MAT,,2
NSLV,S,1
D,ALL,TEMP,%T_WENDU(k)%
ALLSEL,ALL
D,18643,TEMP,%T_WENDU(k)%
D,18561,TEMP,%T_WENDU(k)%
D,18724,TEMP,%T_WENDU(k)%
D,18643,TEMP,%T_WENDU(k)%
D,18805,TEMP,%T_WENDU(k)%
OUTRESM,ALL,1
SOLVE
*ENDDO

23楼2016-03-22 21:36:14

mygt_hit

专家顾问 (职业作家)

专家经验: +362
应助: 438 (硕士)
贵宾: 0.019
金币: 19571.1
散金: 5130
红花: 135
沙发: 3
帖子: 4938
在线: 990.5小时
虫号: 1489764
注册: 2011-11-13
性别: GG
专业: 结构工程
管辖: 土木建筑

内容已删除

知其然，知其所以然。

24楼2016-03-22 22:29:30

mygt_hit

专家顾问 (职业作家)

专家经验: +362
应助: 438 (硕士)
贵宾: 0.019
金币: 19571.1
散金: 5130
红花: 135
沙发: 3
帖子: 4938
在线: 990.5小时
虫号: 1489764
注册: 2011-11-13
性别: GG
专业: 结构工程
管辖: 土木建筑

引用回帖:

23楼: Originally posted by 德群 at 2016-03-22 21:36:14
大神，您说的是这个意思么？
FINISH
/FILNAME,YICIXUNHUAN
/CLEAR
/PREP7
ET,1,SOLID70
MP,KXX,1,1.8
MP,KYY,1,1.8
MP,KZZ,1,1.8
MP,DENS,1,2400
MP,C,1,950
MP ...

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*DIM,T_WENDU,ARRAY,240,1
*DO,k,1,240,1
*SET,T_WENDU(k),10*SIN(0.0262*k)
*ENDDO

*DO,k,1,240,1
TIME,k
TIMINT,ON
AUTOTS,ON
KBC,1
DELTIM,13,2,50
ALLSEL,ALL
DDELE,ALL,TEMP
ALLSEL,ALL
VSEL,S,MAT,,2
NSLV,S,1
D,ALL,TEMP,T_WENDU(k)
ALLSEL,ALL
D,18643,TEMP,T_WENDU(k)
D,18561,TEMP,T_WENDU(k)
D,18724,TEMP,T_WENDU(k)
D,18643,TEMP,T_WENDU(k)
D,18805,TEMP,T_WENDU(k)
OUTRESM,ALL,1
SOLVE
*ENDDO

赞一下(1人)

知其然，知其所以然。

25楼2016-03-22 22:33:17

德群

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 60.3
红花: 1
帖子: 74
在线: 254小时
虫号: 3524726
注册: 2014-11-07
专业: 结构工程

引用回帖:

25楼: Originally posted by mygt_hit at 2016-03-22 22:33:17
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*DIM,T_WENDU,ARRAY,240,1
*DO,k,1,240,1
*SET,T_WENDU(k),10*SIN(0.0262*k)
*ENDDO

*DO,k,1,240,1
TIME,k
TIMINT,ON
AUTOTS,ON
KBC,1
DELTIM,13,2,50
ALLSEL, ...

谢谢大神，我运行下试试

26楼2016-03-22 22:38:05

shgao20

专家顾问 (职业作家)

专家经验: +809
仿真EPI: 1
应助: 1311 (讲师)
金币: 13381.4
红花: 267
帖子: 3531
在线: 473.4小时
虫号: 3644138
注册: 2015-01-14
专业: 动力学与控制
管辖: 仿真模拟

看你们讨论很热烈，就没有参与。不过，看来没有讨论到点子上。而且，以前我提的建议，似乎你也没有考虑。这里不想多说了，给你一个经过修改的命令流试试吧，修改的地方并不多，以 !*** 作为标记：

FINISH
/FILNAME,YICIXUNHUAN, 1
/CLEAR
/PREP7
ET,1,SOLID70
MP,KXX,1,1.8
MP,KYY,1,1.8
MP,KZZ,1,1.8
MP,DENS,1,2400
MP,C,1,950
MP,KXX,2,0.61
MP,KYY,2,0.61
MP,KZZ,2,0.61
MP,DENS,2,996
MP,C,2,4185
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
BLOCK,0,200E-3,0,200E-3,0,600E-3
WPAVE,0,0,400E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

WPROTA,0,-90,0
WPAVE,0,0,50E-3
CSYS,4
VSBW,ALL
WPAVE,0,0,100E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

WPROTA,0,0,90
WPAVE,0,0,50E-3
CSYS,4
VSBW,ALL
WPAVE,0,0,100E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

ALLSEL,ALL
VSEL,S,LOC,Z,0,400E-3
VSEL,R,LOC,X,50E-3,150E-3
VSEL,R,LOC,Y,50E-3,150E-3
VATT,1,,1
ALLSEL,ALL
VSEL,U,MAT,,1,1,1
VATT,2,,1

ALLSEL,ALL
lesize,all,,,10
ALLSEL,ALL
MSHKEY,1
VMESH,ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
save
fini

/SOLU
ALLSEL
VSEL,S,MAT,,1
NSLV,S,1
D,ALL,TEMP,20
ALLSEL,ALL
VSEL,S,MAT,,2
NSLV,S,1
D,ALL,TEMP,10
ALLSEL,ALL
ANTYPE, TRANS
TIMINT,OFF
TIME,0.01
!  DELTIM,0.01, 0.02,0.005    !***  设置为一个子步即可，不用 DELTIM 命令，改为下一句
nsubst, 1,1,1
SOLVE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

*DIM,T_WENDU,ARRAY,240,1
*DO,k,1,240,1
*SET,T_WENDU(k),10*SIN(0.0262*k)
*ENDDO

*DO,k,1,240,1
TIME,k
TIMINT,ON
AUTOTS,ON
!  DELTIM,13,2,50          !***  时间间隔为 1，设置的时间子步长大于 1，难以执行，改为下一句
nsubst, 1,1,1
DDELET,ALL,TEMP
VSEL,S,MAT,,2
NSLV,S,1
!  BF,ALL,TEMP,T_WENDU(k)
D,ALL,TEMP,T_WENDU(k)          !***  还是使用 D 命令
OUTRESM,ALL,1
allsel,all
SOLVE       !***  最关键的，求解必须放到循环内部
*ENDDO

求解后，最终结果，材料 2 (不是 TYPE,2) 的最小温度只有  0.048 左右，基本与 T_WENDU(240) 的值相同。

» 本帖已获得的红花（最新10朵）

德群

27楼2016-03-23 07:17:36

德群

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 60.3
红花: 1
帖子: 74
在线: 254小时
虫号: 3524726
注册: 2014-11-07
专业: 结构工程

送红花一朵

引用回帖:

27楼: Originally posted by shgao20 at 2016-03-23 07:17:36
看你们讨论很热烈，就没有参与。不过，看来没有讨论到点子上。而且，以前我提的建议，似乎你也没有考虑。这里不想多说了，给你一个经过修改的命令流试试吧，修改的地方并不多，以 !*** 作为标记：

FINISH
/FILN ...

嗯嗯，问题解决了，不知道说什么感谢您。谢谢您大神，真的帮了我大忙了！！！！

28楼2016-03-23 09:14:43

德群

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 60.3
红花: 1
帖子: 74
在线: 254小时
虫号: 3524726
注册: 2014-11-07
专业: 结构工程

引用回帖:

大神，还要麻烦您个事，我老师想要，2个小时的降温冻结过程中混凝土的温度能下降到最低-15度，然后融化升温的过程中，混凝土的最高温度能够在8度左右，这个过程怎么实现啊

29楼2016-03-23 22:09:19

德群

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 60.3
红花: 1
帖子: 74
在线: 254小时
虫号: 3524726
注册: 2014-11-07
专业: 结构工程

引用回帖:

我从新改了下命令流，有两个问题
1 无论我温度设置的多低，混凝土内部的温度永远都是19度，我最后设置的温度为-50

2，算到120步的时候自动停止了，是不是我的if语句有问题
命令流如下FINISH
/FILNAME,YICIXUNHUAN, 1
/CLEAR
/PREP7
ET,1,SOLID70
MP,KXX,1,1.8
MP,KYY,1,1.8
MP,KZZ,1,1.8
MP,DENS,1,2400
MP,C,1,950
MP,KXX,2,0.61
MP,KYY,2,0.61
MP,KZZ,2,0.61
MP,DENS,2,996
MP,C,2,4185
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
BLOCK,0,200E-3,0,200E-3,0,600E-3
WPAVE,0,0,400E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

WPROTA,0,-90,0
WPAVE,0,0,50E-3
CSYS,4
VSBW,ALL
WPAVE,0,0,100E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

WPROTA,0,0,90
WPAVE,0,0,50E-3
CSYS,4
VSBW,ALL
WPAVE,0,0,100E-3
CSYS,4
VSBW,ALL
WPCSYS,-1,0

ALLSEL,ALL
VSEL,S,LOC,Z,0,400E-3
VSEL,R,LOC,X,50E-3,150E-3
VSEL,R,LOC,Y,50E-3,150E-3
VATT,1,,1
ALLSEL,ALL
VSEL,U,MAT,,1,1,1
VATT,2,,1

ALLSEL,ALL
lesize,all,,,10
ALLSEL,ALL
MSHKEY,1
VMESH,ALL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
save
finiSH

/SOLU
ALLSEL
VSEL,S,MAT,,1
NSLV,S,1
D,ALL,TEMP,20
ALLSEL,ALL
VSEL,S,MAT,,2
NSLV,S,1
D,ALL,TEMP,10
ALLSEL,ALL
ANTYPE, TRANS
TIMINT,OFF
TIME,0.01

nsubst, 1,1,1
SOLVE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

*DO,k,1,240,1
*IF,K,LE,120,THEN
TIME,k
TIMINT,ON
AUTOTS,ON
nsubst, 1,1,1
DDELET,ALL,TEMP
VSEL,S,MAT,,2
NSLV,S,1
D,ALL,TEMP,-50          !***  还是使用 D 命令
OUTRESM,ALL,1
allsel,all
SOLVE
*ELSEIF,K,GT,120
TIMINT,ON
AUTOTS,ON
nsubst, 1,1,1
DDELET,ALL,TEMP
VSEL,S,MAT,,2
NSLV,S,1
D,ALL,TEMP,10       !***  还是使用 D 命令
OUTRESM,ALL,1
allsel,all
SOLVE
*ENDIF
*ENDDO

-50度.png

30楼2016-03-23 23:52:20