应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » ANSYS 热分析中外界环境温度随时间变化如何求解
533/6返回列表上一页123456下一页
查看: 3258  |  回复: 52
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

德群

新虫 (小有名气)
引用回帖:
19楼: Originally posted by mygt_hit at 2016-03-22 19:00:43
还是table用法上的问题，但是具体怎么用我也不太清楚。
BF加温度时，应该引用table的名，而不是把table当数组来用。而且*DIM定义talbe也有问题，应该是*DIM,T_WENDU,TABLE,240,1
可以试试BF直接加具体数值的温度 ...

我就把时间设为time 180，得到的结果温度还是稳定在那个值，大神，是不是我删除有问题，他前面提示我一个warning
*** WARNING ***                     
When using the DDELE command with a first order system in a transient  (especially following a ramp change), the backward Euler method (TINTP,,,,1.0) usually gives more stable and accurate results.
不知道和这个有没有关系。
麻烦大神了！
ANSYS 热分析中外界环境温度随时间变化如何求解
时间143的温度.png
赞一下 回复此楼
21楼2016-03-22 19:16:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)
内容已删除
回复此楼
知其然，知其所以然。
24楼2016-03-22 22:29:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

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

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

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

FINISH
/FILN ...

嗯嗯，问题解决了，不知道说什么感谢您。谢谢您大神，真的帮了我大忙了！！！！
赞一下 回复此楼
28楼2016-03-23 09:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
引用回帖:
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 生热 ...

大神，还要麻烦您个事，我老师想要，2个小时的降温冻结过程中混凝土的温度能下降到最低-15度，然后融化升温的过程中，混凝土的最高温度能够在8度左右，这个过程怎么实现啊
赞一下 回复此楼
29楼2016-03-23 22:09:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
引用回帖:
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, ...

我从新改了下命令流，有两个问题
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
ANSYS 热分析中外界环境温度随时间变化如何求解-1
-50度.png
赞一下 回复此楼
30楼2016-03-23 23:52:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 德群 的主题更新
533/6返回列表上一页123456下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭