- 积分
- 0
- 注册时间
- 2009-4-11
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2015-3-19 05:30:56
|
显示全部楼层
来自 江苏南京
不知道这两列数组(第一列存放外边界的节点号,第二列准备存放相应的温度。)如何定义和提取数值,以及后续如何读取外边界的温度。
熊少 5:26:53
!模拟华新丽华能源桩在竖直方向上的热演化情况,3D,以换热管内壁热流密度为控制条件,如何是热力场研究应该采用SOLID186
FINISH !退出当前处理器
/CLEAR,START !重新开始一个新的分析
/FILNAME,Try !定义工作名
jpgprf,500,100,1 !将背景调成白色
/prep7 !进入前处理器
et,1,plane55 !定义体单元3-D 20-Node Thermal Solid
!定义换热管热物参数
MP,KXX,1,0.42 !定义pe管的导热系数P180
MP,C,1,1465 !定义pe管的比热容
MP,DENS,1,1100 !定义pe管的密度
blc4,0,0,4,20
mshkey,0 !自由网格划分,不选择映射网格划分是因为通过命令nsla从面上选择点Valid only if the nodes were generated by an area meshing operation [AMESH, VMESH]
allsel
lesize,all,0.5 !网格尺寸为0.5
mat,1 !桩材料
amesh,all
/solu
ANTYPE,TRANS !瞬态分析类型
TIMINT,ON,THERM !打开瞬态效应P245
AUTOTS,ON !打开自动时间步长p240
OUTRES,,ALL !写入每个子步
timint,off !先进行稳态分析,确定初始地温场
time,0.01
!根据监测数据,按层施加温度荷载计算初始温度
nsel,s,loc,y,20
d,all,temp,5
nsel,s,loc,y,10
d,all,temp,10
nsel,s,loc,y,0
d,all,temp,20
allsel
solve
/post1
nsel,s,loc,x,4 !选择所有节点
*get,ncount,node,0,count !获得已选节点的总数Ncount
*dim,narray,array,ncount,2 !定义ncount*2的数组,一列存储节点编号,另一列存储节点温度值
*get,nd,node,0,num,min !得到最小的节点编号
*do,i,1,ncount,1
narray(i,1)=nd
*get,md,node,nd,temp
narray(i,2)=md
nd=ndnext(nd) !读出下一个节点编号
*enddo
finish
/solu
antype,,restart
timint,on !打开瞬态
time,60 !计算600s
deltim,5,5,60
ddele,all,temp !删除温度荷载
NSEL,S,LOC,X,4
D,ALL,TEMP,%narray%
sfl,4,hflux,156.1
allsel
solve
save,,
FINISH
但是不对,请大神帮我看看吧!我对数组不太了解(请问这种情况是用数值数组or表格数组?能简单给出相应的代码,那就nice了!)
我突然想到一种最简单的办法:将边界上的节点反选,然后删除剩余选中的节点温度:
allsel
nsel,u,loc,x,4
ddele,all,temp
...
|
|