本帖最后由 daytime123 于 2010-1-15 21:02 编辑
二、本人使用的手动输入节点温度法
Finish !(从此处到“——”处命令无变化)
/clear
/filname,example2,1
/prep7
et,1,plane55,,,1
mp,kxx,1,1.2
rectng,0.2,0.8,0,0.05
type,1
lsel,s,line,,1,3,2
lesize,all,,,10
lsel,s,line,,2,4,2
lesize,all,,,2
amesh,1
finish
/solu
antype,static
lsel,s,line,,4
nsll,s,1
d,all,temp,80
lsel,s,line,,2
nsll,s,1
d,all,temp,10
allsel
outpr,basic,all
solve
finish
/post1
plnsol,temp
—————————————————
allsel
nsel,all
*get,n1,node,,num,min !最小节点值
*get,n2,node,,num,max !最大节点值
n3=n2-n1+1 !总节点数
*dim,abc,array,n3 !定义一个记录所有节点温度的数组
*do,i,1,n3
*get,abc(i),node,n1+i-1,temp !给数组赋温度值
*enddo
——————————————————
/prep7 !(从此处到“——”处命令与example1中的相同)
etchg,tts
keyopt,1,3,1
keyopt,1,6,1
mp,ex,1,120e9
mp,alpx,,1.3e-6
mp,nuxy,1,.3
lsel,s,line,,3
nsll,s,1
cp,1,uy,all
allsel
cp,2,ux,1,14,24
cp,3,ux,2,12,13
finish
/solu
antype,static
d,1,uy,,,11
——————————————————
*do,i,1,n3
bf,n1+i-1,temp,abc(i) !给结构单元上的所有节点赋温度值
*enddo
——————————————————
Solve !(从此处到结束的命令与example1中的相同)
finish
/post1
plnsol,s,x
nsel,s,loc,x,0.2
nsel,r,loc,y,0.05
*get,node1,node,,num,max
nsel,s,loc,x,0.8
nsel,r,loc,y,0.05
*get,node2,node,,num,max
allsel
path,spath1,2,,30
ppath,1,node1
ppath,2,node2
pdef,sx,s,x
plpath,sx
大家有兴趣的话可以分别运行一下这两段程序,会发现两种解法得出来的应力分布图是一样的。在这里我只张贴出沿矩形最上面那条边的温度分布曲线:
1、
2、
我们发现它们的数值是完全一样的。这说明我使用的这种方法可行。
这仅仅是个简单的例子,用来说明我的思路,有不足或错误的地方还望大家不吝赐教! |