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

24小时热门版块排行榜    

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

shgao20

专家顾问 (职业作家)
1  先回答 30 楼的问题吧：
    你的循环中，time,k 命令应该放到 *if 命令之前，否则，当 k > 120 时，time 一直是 120，就无法继续计算了。
    修改后，240 秒时，中心点温度是 16.43。

  2  对于 29 楼的问题，在计算完成 240 秒的载荷步之后，可以试试如下命令流 (大致)：
      time, 2*3600+240      !  新的求解时间，在原来时间上加上 2 小时。如不对，自己改一下
      ddele,temp, all           !  删除前面的载荷
      d, 18521, temp, -15   !  定义新的载荷 - 中心点温度
         然后可以自己添加其它载荷条件
      nsubst,  100, 10000,10  !  设置子步数
         其它求解设置同前，如有需要可以自己修改
      solve

      time, 7440+???         !  新的求解时间，自己确定，必须大于 7440
      ddele,temp, all           !  删除前面的载荷，如有其它载荷，应一起删除
      d, 18521, temp, 8     !  定义新的载荷 - 中心点温度
         然后可以自己添加其它载荷条件
      nsubst,  100, 10000,10  !  设置子步数，如果新的时间间隔很大或很小，可以自己修改为合理的子步数
         其它求解设置同前，如有需要可以自己修改
      solve

  没有验证过，自己试试，完善一下吧。
赞一下(1人) 回复此楼
31楼2016-03-24 09:05:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)
一点建议
29楼和30楼第一问：尝试把瞬态温度分析变为稳态多荷载步分析，看看效果。
30楼第二问：问题出在*ELSEIF,K,GT,120 这条命令之后了，进入这个分支后TIME没有进行更新，还是一直停留在120。因此改进方法有二：1. 如楼上所说，将TIME, K命令放到IF之前；2. *IF,K,LE,120,THEN 和 *ELSEIF,K,GT,120之后都加TIME, K，更新时间步。
赞一下(1人) 回复此楼

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

德群 德群
知其然，知其所以然。
32楼2016-03-24 11:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)
引用回帖:
32楼: Originally posted by mygt_hit at 2016-03-24 11:55:28
一点建议
29楼和30楼第一问：尝试把瞬态温度分析变为稳态多荷载步分析，看看效果。
30楼第二问：问题出在*ELSEIF,K,GT,120 这条命令之后了，进入这个分支后TIME没有进行更新，还是一直停留在120。因此改进方法有二 ...

以混凝土温度为目标，需要先预估一下周围温度需要达到多少，然后往上加。或者迭代实现，但还需要考虑时间、温度扩散速度等。

另，周围温度-50，中心混凝土还是19度，模型和材料性质是否有问题，需要检查一下。
赞一下(1人) 回复此楼

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

德群
知其然，知其所以然。
33楼2016-03-24 12:01:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
引用回帖:
31楼: Originally posted by shgao20 at 2016-03-24 09:05:05
1  先回答 30 楼的问题吧：
    你的循环中，time,k 命令应该放到 *if 命令之前，否则，当 k > 120 时，time 一直是 120，就无法继续计算了。
    修改后，240 秒时，中心点温度是 16.43。

  2  对于 29 楼 ...

谢谢大神，
第一个问题把time换到上面去就好用了，昨天弄了好久怎么也弄不出来这部分。
第二个问题，我把时间延长之后就可以使用了。

但是现在有个小问题，就是我的最终设定的time是14400，但是在read result 里面得不到14400这个时间，在12984的时候就停止了，这个是怎么回事啊？
下面是我的命令流
*DO,k,60,14400,60
TIME,k
*IF,K,LE,7200,THEN
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,20*SIN(K/2293+1.57)
ALLSEL,ALL
D,18643,TEMP,20*SIN(K/2293+1.57)
D,18561,TEMP,20*SIN(K/2293+1.57)
D,18724,TEMP,20*SIN(K/2293+1.57)
D,18643,TEMP,20*SIN(K/2293+1.57)
D,18805,TEMP,20*SIN(K/2293+1.57)
OUTRESM,ALL,1
SOLVE
*ELSEIF,K,GT,7200
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,5*SIN(K/2293+1.57)
ALLSEL,ALL
D,18643,TEMP,5*SIN(K/2293+1.57)
D,18561,TEMP,5*SIN(K/2293+1.57)
D,18724,TEMP,5*SIN(K/2293+1.57)
D,18643,TEMP,5*SIN(K/2293+1.57)
D,18805,TEMP,5*SIN(K/2293+1.57)
OUTRESM,ALL,1
SOLVE
*ENDIF
*ENDDO
麻烦大神再帮我看一下
ANSYS 热分析中外界环境温度随时间变化如何求解
荷载步.jpg
赞一下 回复此楼
34楼2016-03-24 17:18:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
引用回帖:
33楼: Originally posted by mygt_hit at 2016-03-24 12:01:24
以混凝土温度为目标，需要先预估一下周围温度需要达到多少，然后往上加。或者迭代实现，但还需要考虑时间、温度扩散速度等。

另，周围温度-50，中心混凝土还是19度，模型和材料性质是否有问题，需要检查一下。...

大神，谢谢您一直帮我解决问题！
第一点，经过您说的那种修改已经没问题了，昨天弄它弄了半夜，都快放弃了，没想到就换一下位置。果然自己不太行
第二点，我又重新过了下模型发现我把时间弄错了，因为导热系数是w/(m*c)所以时间是秒，当时没注意到这点，所以那时候计算的是240秒的温度变化情况。

现在还有一个小问题，想请教下大神
我把time最终设为14400 秒，但是在read result里面取结果的时候出现了些问题，就是我得不到14400秒的结果。只能到12984秒就停止了，不知道怎么回事。
下面是我从新修改的命令流，不知道是不是我的那个deltim有问题。麻烦大神了！
*DO,k,60,14400,60
TIME,k
*IF,K,LE,7200,THEN
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,20*SIN(K/2293+1.57)
ALLSEL,ALL
D,18643,TEMP,20*SIN(K/2293+1.57)
D,18561,TEMP,20*SIN(K/2293+1.57)
D,18724,TEMP,20*SIN(K/2293+1.57)
D,18643,TEMP,20*SIN(K/2293+1.57)
D,18805,TEMP,20*SIN(K/2293+1.57)
OUTRESM,ALL,1
SOLVE
*ELSEIF,K,GT,7200
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,5*SIN(K/2293+1.57)
ALLSEL,ALL
D,18643,TEMP,5*SIN(K/2293+1.57)
D,18561,TEMP,5*SIN(K/2293+1.57)
D,18724,TEMP,5*SIN(K/2293+1.57)
D,18643,TEMP,5*SIN(K/2293+1.57)
D,18805,TEMP,5*SIN(K/2293+1.57)
OUTRESM,ALL,1
SOLVE
*ENDIF
*ENDDO
ANSYS 热分析中外界环境温度随时间变化如何求解-1
荷载步情况.jpg
赞一下 回复此楼
35楼2016-03-24 17:26:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
送红花一朵
引用回帖:
32楼: Originally posted by mygt_hit at 2016-03-24 11:55:28
一点建议
29楼和30楼第一问：尝试把瞬态温度分析变为稳态多荷载步分析，看看效果。
30楼第二问：问题出在*ELSEIF,K,GT,120 这条命令之后了，进入这个分支后TIME没有进行更新，还是一直停留在120。因此改进方法有二 ...

大神， 能不能麻烦您再帮我看一下这个部分，这两天都在改这个，可是怎么改都改不好。麻烦了大神！
赞一下 回复此楼
36楼2016-03-25 13:50:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
送红花一朵
引用回帖:
33楼: Originally posted by mygt_hit at 2016-03-24 12:01:24
以混凝土温度为目标，需要先预估一下周围温度需要达到多少，然后往上加。或者迭代实现，但还需要考虑时间、温度扩散速度等。

另，周围温度-50，中心混凝土还是19度，模型和材料性质是否有问题，需要检查一下。...

大神，我在进行温度和力学耦合的模拟，就是第一个荷载步的模拟，但是我现在有些问题。
首先，我定义的材料2是水，但是水是没有弹性模量的，但是不输入又不行，这个该怎么办？
其次，我把水的弹性模量设为1，运算出现一个warning， Coefficient ratio exceeds 1.0e8 - Check results. 是不是和我那个弹性模量设置的有关系。
最后，当我点开read result 的时候想读取荷载步的时候显示的是 An error occurred while attempting to open the results file YICIXUNHUAN.rst.  节点的应变图也得不到，这又是什么原因？
下面是我的命令流
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
!  DELTIM,0.01,    0.02,0.005     !***  设置为一个子步即可，不用 DELTIM 命令，改为下一句
SOLVE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

/PREP7
ETCHG,TTS
MP,EX,1,3E10
MP,ALPX,1,10E-6
MP,NUXY,1,0.2

MP,EX,2,1
MP,ALPX,2,0.0002
MP,NUXY,0.5

ASEL,S,AREA,,90
NSLA,S,1
D,ALL,UZ,0
NSEL,ALL

FINISH
/SOLU
TREF,20
LDREAD,TEMP,,,,,YICIXUNHUAN,RTH

/PBF,TEMP,,1
SOLVE
FINISH
麻烦大神帮我看一下
赞一下 回复此楼
37楼2016-03-26 00:15:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)
对 35 楼的问题：
    先打开输出窗口 (output 窗口)，看一下有什么警告或错误信息，在看看 .err 文件，看一下有什么警告或错误信息，然后找出问题进行修改。
    或者，试试看：在命令流的 /prep7 之前，加一句：
        /config, nres, 10000

  对 37 楼的问题：
  在温度场分析完成后，转换到结构分析，这是实验 Ekill 命令，杀死作用的水单元，再进行求解，水就不起作用了。当然，看结果时，也需要把水单元去掉 (Unselected)。
赞一下 回复此楼

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

德群
38楼2016-03-26 09:14:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
送红花一朵
引用回帖:
38楼: Originally posted by shgao20 at 2016-03-26 09:14:29
对 35 楼的问题：
    先打开输出窗口 (output 窗口)，看一下有什么警告或错误信息，在看看 .err 文件，看一下有什么警告或错误信息，然后找出问题进行修改。
    或者，试试看：在命令流的 /prep7 之前，加一句： ...

大神，第一个问题解决了！
第二个问题就是我想用水的体积弹性模量这样靠谱么？还是只能把水删去？
赞一下 回复此楼

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

德群
39楼2016-03-26 11:01:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

德群

新虫 (小有名气)
送红花一朵
引用回帖:
39楼: Originally posted by 德群 at 2016-03-26 11:01:20
大神，第一个问题解决了！
第二个问题就是我想用水的体积弹性模量这样靠谱么？还是只能把水删去？...

大神，你好，还得麻烦您一下，我按照您的说法做了把水的单元用ekill命令杀死，得到的结果好像也差不多。但是有两个问题，麻烦大神帮我分析下
1.我想在混凝土的下表面施加ASEL,S,LOC,Z,0    NSLA,S,1   D,ALL,UZ,0就是在Z=0的平面固定Z方向的位移， 但是弹出warning   An unknown degree of freedom label (5) was found on the D command.和 No valid degree of freedom labels were input. 这两个。然后D命令被忽略

2，杀死单元之后，我最后得到的结果永远都是最后一个时间步的位移，用outres命令也找不到之前时间步的数据，这又是怎么回事？

下面是我的命令流
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
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

/SOLU
TIME,1  
TIMINT,ON
AUTOTS,ON
KBC,1
NSUBST,1,1,1
ALLSEL,ALL  
DDELE,ALL,TEMP                     
ALLSEL,ALL
VSEL,S,MAT,,2
NSLV,S,1
D,ALL,TEMP,5
ALLSEL,ALL
D,18643,5
D,18561,5
D,18724,5
D,18643,5
D,18805,5
OUTRESM,ALL,1
SOLVE

/PREP7
ETCHG,TTS
MP,EX,1,3E10
MP,ALPX,1,10E-6
MP,NUXY,1,0.2

MP,EX,2,2.18E9
MP,ALPX,2,0.0002
MP,NUXY,0.2

ALLSEL,ALL
ASEL,S,LOC,Z,0
NSLA,S,1
D,ALL,UZ,0                          !此处命令被忽略了

/SOLU
TIME,0.01
NLGEOM,ON
NROPT,FULL
ESTIF,,
ESEL,S,MAT,,2
EKILL,ALL
ALLSEL,ALL
ESEL,S,LIVE
NSLE,S
LDREAD,TEMP,,,0.01,,MONI,RTH
SOLVE
FINISH

/SOLU
TIME,1
NLGEOM,ON
NROPT,FULL
ESTIF,,
ESEL,S,MAT,,2
EKILL,ALL
ALLSEL,ALL
ESEL,S,LIVE
NSLE,S
LDREAD,TEMP,,,1,,MONI,RTH      ！我有两个荷载步，但是结果只是      显示最后的一步荷载，应该用什么命令？

SOLVE
FINISH
麻烦大神再指导下，谢谢您！
赞一下 回复此楼
40楼2016-03-27 23:34:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 德群 的主题更新
534/6返回列表上一页123456下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭