popohopo 发表于 2012-3-17 16:49:08

关于热分析中提取所属某面上节点的热流密度分量值的问题

关于热分析中提取所属某面上的所有节点的热流密度分量值或者热梯度分量值的问题
背景:
当对某一器件进行完热分析之后,我想提取某一面上所有节点的热流密度或者是热梯度分量值,然后保存到一个.TXT格式的文本文件中,已方便后期用Matlab继续做处理。
问题:
我先将需要提取节点热流密度分量值的面单独选中,然后在用NSEL命令将所属该面上的所有节点选中,之后我使用下面这段命令流进行提取,保存数据

*GET,NUMBER_NODE,NODE,0,COUNT                   !获取所选择的节点总数
*DIM,QW_TF,ARRAY,NUMBER_NODE,3               !定义1个数组存放数据
*DO,I,1,NUMBER_NODE,1
         *GET,MIN_NODE,NODE,0,NUM,MIN            !获取最小的节点编号
         QW_TG(I,1)=NX(MIN_NODE)                         !节点的X坐标放数组第1列
         QW_TG(I,2)=NZ(MIN_NODE)                         !节点的Z坐标放数组第2列
         *GET,QW_TG(I,3),NODE,MIN_NODE,TF,X       !节点的热流密度X分量放数组第3列
         NSEL,U,NODE,,MIN_NODE                              !从选择集中剔除已提取坐标与热流量的节点
*ENDDO
*CFOPEN,MAX_TF_NODE_INFORMATION,TXT,,
*VWRITE,QW_TF(1,1),QW_TF(1,2),QW_TF(1,3)
(F15.5' 'F15.5' 'F15.5)
*CFCLOS

运行时,我发现运行的速度相当的慢且没有一次成功的将热流密度分量值提取,而我对上述命令流稍作修改(把 *GET,QW_TG(I,3),NODE,MIN_NODE,TF,X改为了 *GET,QW_TG(I,3),NODE,MIN_NODE,TEMP),用来提取该面上所有节点的温度值时,却提取速度很快!不知道为什么?
是不是因为对于热分析单元而言,温度是其基本解,而热流密度以及热梯度是其导出解,在提取的时候只能提取节点温度值,而不可以将节点的热流密度分量与热梯度分量值提取出来?
若是哪位懂的,或者做过类似处理的人能给予解答,先谢了!!!
页: [1]
查看完整版本: 关于热分析中提取所属某面上节点的热流密度分量值的问题