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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 怎么才能吧ansys中每个节点的温度提取提取出来导入matlab处理
333/4返回列表上一页1234下一页
查看: 2115  |  回复: 32
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

晓嘉嘉

新虫 (小有名气)
引用回帖:
20楼: Originally posted by mygt_hit at 2017-03-20 16:41:15
1. 试试*get,Tem_Node(i,j),NODE,nodes(i,j),TEMP
2. 试试nsol命令
3. 采用PRNSOL列出结果，拷到TXT，导入MATLAB变形为151*26矩阵...

FINI
/clea
/filname,thermal_net_condANDconv
/title,wafer and water
!!!!!!!!!前处理!!!!!!!!!!
/prep7
et,1,link33
r,1,1e-12
MPTEMP,1,300,577,854,1131,1408,1685,
MPTEMP,7,1687, 1688,1919 ,2150, 2381,2612,
MPDATA,DENS,1,1,2.3475E3,2.3402E3,2.3329E3,2.3256E3,2.3183E3,2.3111E3,
MPDATA,DENS,1,7,2.3110E3,2.5798E3,2.5317E3,2.4663E3,2.3838E3,2.2841E3,
MPDATA,KXX,1,1,148.7562, 62.5523,39.6026,28.9729,22.8419,18.8525,
MPDATA,KXX,1,7,18.8287,18.8169,16.4286,14.5783,13.1025,11.898,
MPDATA,C,1,1,0.7007e3,0.7702e3,0.8402e3,0.9106e3,0.9814e3,1.0527e3,
MPDATA,C,1,7,1.0532e3,0.9427e3,0.9606e3,0.9861e3,1.0202e3,1.0648e3,
MPTEMP,1,300,577,854,1131,1408,1685
MPTEMP,7, 1688,1919 ,2150, 2381,2612,2843,
MPDATA,ENTH,1,1,0,4.775e8,9.986e8,1.5634e9,2.1718e9,2.8239e9,
MPDATA,ENTH,1,7,3.3774e9,3.9392e9,4.501e9,5.063e9,5.6246e9,6.1864e9,
et,2,link34
r,2,1e-12
MPTEMP,1,300,577,854,1131,1408,1685
MPDATA,hf,2,1,120,1000,8000,15000,18000,20000
et,3,link31
r,3,1e-12,1,1,5.67e-8
mp,emis,3,0.66
et,4,link34
r,4,1e-12
MPTEMP,1,300,577,854,1131,1408,1685
MPDATA,hf,4,1,120,1000,8000,15000,18000,20000
*DO,I,1,11
*DO,J,1,6,1
N,100000+6*(I-1)+J,2.00000e-06*J,2.00000e-06*(I-1)
*ENDDO
*ENDDO
mat,1
type,1
real,1
*DO,I,1,11
*DO,J,1,5,1
E,100000+6*(I-1)+J,100001+6*(I-1)+J
*ENDDO
*ENDDO
*DO,I,1,6
*DO,J,1,10,1
E,100000+I+6*(J-1),100000+I+6*J
*ENDDO
*ENDDO
nsel,all
save
physics,write,thermal,,,
/solu
physics,read,thermal,,,
antype,trans
ic,all,temp,300
nsel,s,loc,x,6*2e-06
d,all,temp,300
nsel,s,loc,y,10*2e-06
d,all,temp,300
allsel
eplot
timint,off
nropt,full
solcontrol,on
pred,on,,on
NCNV,0
autots,on
outres,basic,LAST
toffst,273
tintp,,,,1,
time,1e-11
solve
*DO,I,1,2,1
*do,jj,1,5,1
*do,ii,1,11,1
nsel,s,node,,100000+(ii-1)*6+jj,100001+(ii-1)*6+jj
esln,s,1
esel,r,type,,1
ex=2e-6*jj
ey=2e-6*(ii-1)
bfe,all,hgen,1,5.4126e+016*(1-0.380)*exp(-4*(ex-0)**2/(20e-6)**2)*exp(-2792.79*ey)
*enddo
*enddo
*do,jj,1,6,1
*do,ii,1,9,1
nsel,s,node,,100000+(ii-1)*6+jj,100000+ii*6+jj,6
esln,s,1
esel,r,type,,1
ex=2e-6*(jj-1)
ey=2e-6*ii
bfe,all,hgen,1,5.4126e+016*(1-0.380)*exp(-4*(ex-0)**2/(20e-6)**2)*exp(-2792.79*ey)
*enddo
*enddo
timint,on
lnsrch,auto
time,I*.05e-9
allsel
solve
esel,all
bfedele,all,hgen
*ENDDO
/post1
allsel,all
!set,list
SET,LAST
/EFACET,1
PLNSOL,TEMP
/auto，1
file,alijiajia,rst
*DIM,Tem_Node,ARRAY,11,6
*DIM,nodes,ARRAY,11,6
*do,I,1,11,1
*do,J,1,6,1
nodes(I,J)=100000+6*(I-1)+J
*GET,Tem_Node(I,J),NODE,nodes(I,J),TEMP
*enddo
*enddo
*CFOPEN,alijiajia,txt
*DO,I,1,11,1
*Vwrite,(Tem_Node(I,1)),(Tem_Node(I,2)),(Tem_Node(I,3)),(Tem_Node(I,4)),(Tem_Node(I,5))(Tem_Node(I,6))
(26f8.5)
*ENDDO
*CFCLOSE
形式如上程序
麻烦老师帮我看看
*get,Tem_Node(i,j),NODE,nodes(i,j),TEMP这个提示出错
赞一下 回复此楼
21楼2017-03-20 20:15:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黑太阳

金虫 (正式写手)
节点有它自己的节点号，好像不能重新定义，

发自小木虫Android客户端
赞一下 回复此楼
22楼2017-03-20 20:44:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晓嘉嘉

新虫 (小有名气)
你说的不能重新定义是什么意思，是提取的时候还是在计算的过程？在list results中是存在结果的，前面的图片可以说明 @mygt_hit

发自小木虫Android客户端
赞一下 回复此楼
23楼2017-03-20 20:51:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晓嘉嘉

新虫 (小有名气)
引用回帖:
22楼: Originally posted by 黑太阳 at 2017-03-20 20:44:18
节点有它自己的节点号，好像不能重新定义，

你说的不能重新定义是什么意思，是提取的时候还是在计算的过程？在list results中是存在结果的，前面的图片可以说明

发自小木虫Android客户端
赞一下 回复此楼
24楼2017-03-20 20:52:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

引用回帖:
21楼: Originally posted by 晓嘉嘉 at 2017-03-20 20:15:46
FINI
/clea
/filname,thermal_net_condANDconv
/title,wafer and water
!!!!!!!!!前处理!!!!!!!!!!
/prep7
et,1,link33
r,1,1e-12
MPTEMP,1,300,577,854,1131,1408,1685,
MPTEMP,7,1687, 1688,1919 ,2150 ...

问题出在file,alijiajia,rst，计算结果中并没有这个.rst文件。注释掉即可。另外最后写入文件命令也有问题，给你调整过来了。


FINI
/clea
/filname,thermal_net_condANDconv
/title,wafer and water
!!!!!!!!!前处理!!!!!!!!!!
/prep7
et,1,link33
r,1,1e-12
MPTEMP,1,300,577,854,1131,1408,1685,
MPTEMP,7,1687, 1688,1919 ,2150, 2381,2612,
MPDATA,DENS,1,1,2.3475E3,2.3402E3,2.3329E3,2.3256E3,2.3183E3,2.3111E3,
MPDATA,DENS,1,7,2.3110E3,2.5798E3,2.5317E3,2.4663E3,2.3838E3,2.2841E3,
MPDATA,KXX,1,1,148.7562, 62.5523,39.6026,28.9729,22.8419,18.8525,
MPDATA,KXX,1,7,18.8287,18.8169,16.4286,14.5783,13.1025,11.898,
MPDATA,C,1,1,0.7007e3,0.7702e3,0.8402e3,0.9106e3,0.9814e3,1.0527e3,
MPDATA,C,1,7,1.0532e3,0.9427e3,0.9606e3,0.9861e3,1.0202e3,1.0648e3,
MPTEMP,1,300,577,854,1131,1408,1685
MPTEMP,7, 1688,1919 ,2150, 2381,2612,2843,
MPDATA,ENTH,1,1,0,4.775e8,9.986e8,1.5634e9,2.1718e9,2.8239e9,
MPDATA,ENTH,1,7,3.3774e9,3.9392e9,4.501e9,5.063e9,5.6246e9,6.1864e9,
et,2,link34
r,2,1e-12
MPTEMP,1,300,577,854,1131,1408,1685
MPDATA,hf,2,1,120,1000,8000,15000,18000,20000
et,3,link31
r,3,1e-12,1,1,5.67e-8
mp,emis,3,0.66
et,4,link34
r,4,1e-12
MPTEMP,1,300,577,854,1131,1408,1685
MPDATA,hf,4,1,120,1000,8000,15000,18000,20000
*DO,I,1,11
*DO,J,1,6,1
N,100000+6*(I-1)+J,2.00000e-06*J,2.00000e-06*(I-1)
*ENDDO
*ENDDO
mat,1
type,1
real,1
*DO,I,1,11
*DO,J,1,5,1
E,100000+6*(I-1)+J,100001+6*(I-1)+J
*ENDDO
*ENDDO
*DO,I,1,6
*DO,J,1,10,1
E,100000+I+6*(J-1),100000+I+6*J
*ENDDO
*ENDDO
nsel,all
save
physics,write,thermal,,,
/solu
physics,read,thermal,,,
antype,trans
ic,all,temp,300
nsel,s,loc,x,6*2e-06
d,all,temp,300
nsel,s,loc,y,10*2e-06
d,all,temp,300
allsel
eplot
timint,off
nropt,full
solcontrol,on
pred,on,,on
NCNV,0
autots,on
outres,basic,LAST
toffst,273
tintp,,,,1,
time,1e-11
solve

*DO,I,1,2,1
*do,jj,1,5,1
*do,ii,1,11,1
nsel,s,node,,100000+(ii-1)*6+jj,100001+(ii-1)*6+jj
esln,s,1
esel,r,type,,1
ex=2e-6*jj
ey=2e-6*(ii-1)
bfe,all,hgen,1,5.4126e+016*(1-0.380)*exp(-4*(ex-0)**2/(20e-6)**2)*exp(-2792.79*ey)
*enddo
*enddo
*do,jj,1,6,1
*do,ii,1,9,1
nsel,s,node,,100000+(ii-1)*6+jj,100000+ii*6+jj,6
esln,s,1
esel,r,type,,1
ex=2e-6*(jj-1)
ey=2e-6*ii
bfe,all,hgen,1,5.4126e+016*(1-0.380)*exp(-4*(ex-0)**2/(20e-6)**2)*exp(-2792.79*ey)
*enddo
*enddo
timint,on
lnsrch,auto
time,I*.05e-9
allsel
solve
esel,all
bfedele,all,hgen
*ENDDO

/post1
allsel,all
!set,list
SET,LAST
/EFACET,1
PLNSOL,TEMP
/auto,1

! file,alijiajia,rst
*DIM,Tem_Node,ARRAY,11,6
*DIM,nodes,ARRAY,11,6
*do,I,1,11,1
*do,J,1,6,1
nodes(I,J)=100000+6*(I-1)+J
*GET,Tem_Node(I,J),NODE,nodes(I,J),TEMP
*enddo
*enddo
*CFOPEN,alijiajia,txt
I=1
*vwrite,Tem_Node(I,1),Tem_Node(I,2),Tem_Node(I,3),Tem_Node(I,4),Tem_Node(I,5),Tem_Node(I,6)
(6F10.3)
*CFCLOSE
赞一下(3人) 回复此楼
知其然，知其所以然。
25楼2017-03-20 21:02:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黑太阳

金虫 (正式写手)
引用回帖:
24楼: Originally posted by 晓嘉嘉 at 2017-03-20 20:52:40
你说的不能重新定义是什么意思，是提取的时候还是在计算的过程？在list results中是存在结果的，前面的图片可以说明
...

提取的时候。ansys提示的是找不到对应的节点，说明你定义的节点号数组与实际存在的节点不对应。

发自小木虫Android客户端
赞一下 回复此楼
26楼2017-03-20 21:16:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晓嘉嘉

新虫 (小有名气)
引用回帖:
25楼: Originally posted by mygt_hit at 2017-03-20 21:02:48
问题出在file,alijiajia,rst，计算结果中并没有这个.rst文件。注释掉即可。另外最后写入文件命令也有问题，给你调整过来了。


FINI
/clea
/filname,thermal_net_condANDconv
/title,wafer and water
!!!!! ...

你好老师，运行之后，在txt中并没有得到结果如图，提示警告如图
怎么才能吧ansys中每个节点的温度提取提取出来导入matlab处理


怎么才能吧ansys中每个节点的温度提取提取出来导入matlab处理-1



发自小木虫Android客户端
赞一下 回复此楼
27楼2017-03-20 21:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晓嘉嘉

新虫 (小有名气)
引用回帖:
26楼: Originally posted by 黑太阳 at 2017-03-20 21:16:23
提取的时候。ansys提示的是找不到对应的节点，说明你定义的节点号数组与实际存在的节点不对应。
...

上面有我回复的程序，你可以帮忙看看么，我感觉没有不对应啊，节点号都是对的

发自小木虫Android客户端
赞一下 回复此楼
28楼2017-03-20 21:27:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
晓嘉嘉: 金币+5, ★★★★★最佳答案, 得到想要的结果 2017-03-20 21:55:58
引用回帖:
27楼: Originally posted by 晓嘉嘉 at 2017-03-20 21:25:18
你好老师，运行之后，在txt中并没有得到结果如图，提示警告如图


...

从菜单直接读入，不要拷到命令窗口运行。

发自小木虫Android客户端
赞一下(1人) 回复此楼

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

晓嘉嘉
知其然，知其所以然。
29楼2017-03-20 21:39:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晓嘉嘉

新虫 (小有名气)
送红花一朵
引用回帖:
29楼: Originally posted by mygt_hit at 2017-03-20 21:39:57
从菜单直接读入，不要拷到命令窗口运行。
...

谢谢，麻烦老师了！！出来了，就是还想请教你一个问题，这两种输入方式为啥会有如此大的区别？！

发自小木虫Android客户端
赞一下 回复此楼
30楼2017-03-20 21:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 晓嘉嘉 的主题更新
333/4返回列表上一页1234下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭