如何选取节点

wangwenna218

我的模型是矩形，我想在最下面的一条y=0的线上加载。
我要把这条线分为5段来移动载荷，在每一段的两端的载荷是中间节点载荷的二分之一。
我的问题在于选取这些节点。总是选不对，有时候会选上好多个。请大家看看是为什么？
/com, units: SI, m - degree C - W - kg -sec - J
/title,bs
/filname,the analysis of a crane
/prep7
et,1,shell57
R,1,0.022  !盖板厚度
R,2,0.006  !腹板厚度
!*****************定义热分析物理参数**********************
mptemp,1,20,250,500,750  !定义材料温度区间(C)
mptemp,,1000,1500, 1700,2500
mpdata,dens,1,1,7820,7700,7610,7550    !定义材料密(kg/m3)
mpdata,dens,1,,7490,7350,7300,7090
mpdata,kxx,1,1,50,47,40,27               !热传导系数(W/(m*K))
mpdata,kxx,1,,30,35,14,14.2
mpdata,c,1,1,460,480,530,675                  !比热容(J/kg/K)
mpdata,c,1,,670,660,780,820
!**********定义应力分析参数****************************                          
mpdata,ex,1,1,2.05e11,1.87e11,1.5e11,0.7e11   !定义弹性模量
mpdata,ex,1,,0.2e11,0.19e11,0.18e11,0.12e11
tb,bkin,1,6                                !定义屈服强度和切变模量
tbtemp,20
tbdata,,2.3e8,2.05e10
tbtemp,250
tbdata,,1.75e8,1.87e10
tbtemp,500
tbdata,,1.3e8,1.5e10
tbtemp,1000
tbdata,,0.25e8,0.2e10
tbtemp,1500
tbdata,,2e6,0.19e10
tbtemp,2000
tbdata,,0.7e6,0.15e10
mpdata,alpx,1,1,1.1e-5,1.22e-5,1.39e-5,1.48e-5            !定义线膨胀系数
mpdata,alpx,1,,1.34e-5,1.33e-5,1.32e-5,1.32e-5
mpdata,prxy,1,1,0.28,0.29,0.31,0.35              !定义泊松比
mpdata,prxy,1,,0.4,0.45,0.48,0.5
!***********建立模型*************
wprot,,-90  
BLC5,0,0,11.53,0.6/2          !上盖板和下盖板
wpoff,1.14,,-1.450  
BLC5,0,0,18.5/2,0.6/2
wpoff,-6.615,,0.65
BLC5,0,0,0.58,0.6/2

WPCSYS,-1,0
k,13,-5.185,0,0.12         !腹板
k,14,-5.185,-0.8,0.12   
k,15,-3.485,-1.450,0.12
k,16,11.53/2,-1.450,0.12
k,17,11.53/2,0,0.12
a,13,14,15,16,17
a,10,14,15,5
a,14,15,8,11
wpoff,,,0.12
ASBW,       1  
wpoff,,-1.45
ASBW,       2   
WPCSYS,-1,0
aglue,all
LCCAT,17,32
LCCAT,18,23
LCCAT,6,13
!*******************************划分网格********************
lsize=0.08
!************上盖板和下盖板
mshkey,0
esize,lsize
amesh,11,12
amesh,1,2
amesh,9,10
amesh,13
!**************腹板
mshkey,0
real,2
esize，lsize
amesh,14
!****************细化分焊缝网格1/3
LREF,17,29,12,2,1,1,1   
LREF,6,,,2,1,1,1
nummrg,all
!************************************非线性求解选项
/solu
antype,4
trnopt,full            !完全瞬态积分法
pred,on                !打开预测校正
solcontrol,on          !自动控制
nropt,full,,on         !完全牛顿-拉普森方法
!cnvtol,u,,0.5          !设置收敛值，位移
!cnvtol,f,,0.05         !设置收敛值，力
cutcontrol,plslimit,10,!设置在一个非线性求解中时间步的减少
!eqslv,sparse,,2        !求解器类型
lnsrch,on              !线性搜索
autots,on              !自动步长        
tintp,,,,1,,           !向后积分        
timint,on              !瞬态效果  
tref,25                !参考温度
kbc,1                 !阶跃性载荷
neqit,50               !规定每个子步中最大迭代次数,默认为25
outres,all
!***********************热源
lsize=lsize/3
u=30
i=300
v=0.006
rb=0.005
k=1/rb**2
qm=0.75*u*i/(rb*rb)
Q=qm*lsize*sqrt(3.1415926/k)
t=sqrt(3.1415926/k)/v

!********************焊接上盖板和腹板的双面焊缝**********************
sf,all,conv,15,25
NWPAVE,147 !***********在节点147处建立局部坐标
wprot,180  
CSYS,4
!***********杀死下盖板
ASEL,S, , ,1,2
ASEL,a, , ,9,10
ASEL,a, , ,13
esla,s,all
eplot
ekill,all
alls
eplot
!***********焊接第一道
d=10.95
n=5
xmin=0                     
xmax=d/n
*do,i,1,n
time,i*t
AUTOTS,0
nsubst,5
nsel,s,loc,y,0
nsel,r,loc,z,0
nsel,r,loc,x,xmin+1e-5,xmax-lsize,lsize
f,all,heat,Q        !施加中间热流率
alls
nsel,s,loc,y,0
nsel,r,loc,z,0
nsel,r,loc,x,xmin
f,all,heat,Q/2    !施加两端热流率
alls
nsel,s,loc,y,0
nsel,r,loc,z,0
nsel,r,loc,x,xmax
f,all,heat,Q/2
alls
solve
fdele,all,heat  !撤掉热流率
xmin=xmin+d/n
xmax=xmax+d/n
*enddo

wangwenna218

我的QQ是1023400904
王文娜

wangwenna218

这么冷清？呵呵
不过好像我解决掉了一点

spiderman216

顶啊，哈哈