24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1162  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

985379629

铁虫 (初入文坛)

[求助] ANSYS温度场边界条件施加的问题已有1人参与

请教各位大佬,目前在研究温度场的问题,需要把辐射换热等效成热流密度施加到结构表面节点上,但是计算方法与结构表面的温度相关,即我每计算完成一个荷载步都需要获取上一荷载步下结构的表面温度,以计算当前荷载步的热流密度大小,再进行施加。我用*get命令循环获取表面每个节点的温度值,计算节点处的热流密度,再循环施加到节点上,此种方式太过繁琐,请教各位大佬有没有更好的办法实现这个过程,以及有没有什么命令能批量获取温度值,下附部分命令流,还有一个问题,就是这个计算会报错,错误信息如下。如果去掉计算初始温度场的部分就不会报错,还不知道原因,我应该如何在计算初始温度场的情况下解决这个报错信息呢。金币数是我的全部家当了,感谢各位大佬。
*** ERROR ***                           
An unexpected error ( SIG$SEGV ) has occurred...  ANSYS internal data   
has been corrupted.  ANSYS is unable to recover and will terminate.     
Previously saved files are unaffected.  Please send the data leading   
to this operation to your technical support provider, as this will      
allow ANSYS, Inc to improve the program.  
                 
/SOLU
ANTYPE,TRANS
TIMINT,OFF
NSUBST,10
TIME,1E-6                   !计算初始温度场
ALLSEL
NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ
ESLN,R
D,ALL,TEMP,32.2           
ALLSEL
NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.05
ESLN,R
D,ALL,TEMP,34.1         
ALLSEL
NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.1
ESLN,R
D,ALL,TEMP,35            
ALLSEL
NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.15
ESLN,R
D,ALL,TEMP,35.6            
ALLSEL
NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.2
ESLN,R
D,ALL,TEMP,36.2         
ALLSEL
NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.215
ESLN,R
D,ALL,TEMP,36.3         
ALLSEL
NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.23
ESLN,R
D,ALL,TEMP,36.3           
ALLSEL
NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.28
ESLN,R
D,ALL,TEMP,36.2            
ALLSEL
SOLVE
DDELE,ALL,ALL  
TIMINT,ON
AUTOTS,ON                                                                                                                                   
NSUBST,25                                                              
NEQIT,50
CSYS,0
KBC,1                                                                                                                                                                        
TIME,1E-6+3600                                                                                                       
DELTIM,300
NROPT,FULL
OUTRES,ALL,ALL
*DO,I,1,3600/300,1
        TIME,300*I
        SFDELE,ALL,ALL
  *DO,J,1,ELEM_NUM            !遍历表面每个单元施加荷载
        ALLSEL
          *GET,TEMPP1,NODE,TOP_NODE(J,1),TEMP                      !获取单元表面四个角节点温度值
                *GET,TEMPP2,NODE,TOP_NODE(J,2),TEMP
                *GET,TEMPP3,NODE,TOP_NODE(J,3),TEMP
                *GET,TEMPP4,NODE,TOP_NODE(J,4),TEMP
                TEMP_ELEM=(TEMPP1+TEMPP2+TEMPP3+TEMPP4)/4   !平均温度当做单表面温度
                ESEL,S,,,TOP_ELEM(J)
                NSLE,R
                NSEL,R,LOC,Y,GDBZ+CAZ+DZBZ       
          SF,ALL,HFLUX,-FAILL*C*((273+TEMP_ELEM)**4)       !计算施加当前荷载步的热流密度
  *ENDDO
  ALLSEL
  SOLVE
*ENDDO
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
985379629: 金币+18, ★★★很有帮助 2022-05-17 16:47:39
你质的是哪个错误?
  高于  DELTIM,300, 30,300, 15, 1 命令,是我写错了,多了一个参数 300.应该改为:
       DELTIM, 30,300, 15, 1
  其中 命令之后的 4 个数据分别是:
    初始的子步长,最大子步长,最小子步长;第 4 个参数 1 表示新的载荷步的初始子步长采用上一个载荷步结束时的子步长。
  具体可见帮助文件中高于 deltim 的内容。
  我的失误造成你的困扰,特表示抱歉。
5楼2022-05-15 08:28:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
看了以下你的命令流,可能有几个问题,给你参考一下:
1.  命令流中多处使用了类似如下的命令流:
    NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.05
    ESLN,R
    D,ALL,TEMP,34.1         
    ALLSEL
因为你是在节点上施加温度载荷,其中的 ESLN 没有用,可以去掉。
2.  命令 NSUBST,25最好写完整,比如改为:
      NSUBST,25, 100,2
    同样,DELTIM,300 最好改为,如:
      DELTIM,300, 30,300, 15, 1
  这列出修改是为了和打开自动时间步相配套。
3.  在计算初始温度前,增加命令:
      KBC, 0
  有利于收敛。
后面的 KBC, 1 最好也改为  KBC,0
4.  命令流中使用了数组 TOP_NODE,但是这段命令流中没有定义过该数组,检查以下,前面是否定义过。
先说这些吧。

» 本帖已获得的红花(最新10朵)

2楼2022-05-11 08:35:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

985379629

铁虫 (初入文坛)

送红花一朵
引用回帖:
2楼: Originally posted by htbbzzg at 2022-05-11 08:35:21
看了以下你的命令流,可能有几个问题,给你参考一下:
1.  命令流中多处使用了类似如下的命令流:
    NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.05
    ESLN,R
    D,ALL,TEMP,34.1         
    ALLSEL
因为你是在节 ...

感谢您的帮助,我按照您的建议修改了,第四个问题我前面是定义过这里没有放上来,请问下这个报错信息一般是什么情况导致的呢?
3楼2022-05-13 11:57:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

985379629

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by htbbzzg at 2022-05-11 08:35:21
看了以下你的命令流,可能有几个问题,给你参考一下:
1.  命令流中多处使用了类似如下的命令流:
    NSEL,S,LOC,Y,DZBZ+CAZ+GDBZ-0.05
    ESLN,R
    D,ALL,TEMP,34.1         
    ALLSEL
因为你是在节 ...

您好,关于您的第二点建议,这个命令流我没看懂,我查了下帮助文档,DELTIM, DTIME, DTMIN, DTMAX, Carry,DELTIME后面是四个数值,您写的这个是什么含义呢?
4楼2022-05-13 16:41:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见