马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我看了很久也没找出原因,烦请各位前辈费心指点:
加荷载的时候有两段程序,第一段可以运行:
csys,0
asel,s,area,,2
esla,s
nsla,s,1
*get,enmax,elem,,num,max
*get,enmin,elem,,num,min
fcum,add
*do,i,enmin,enmax,1
*if,esel(i),eq,1,then
*get,ae,elem,i,area
ze=centrz(i)
p_e=-(ka*gama2*(depth-0.05)+ka*gama1*0.050)
f_tot=p_e*ae
esel,s,elem,,i
nsle,s,corner
*get,nn,node,,count
f_n=f_tot/nn
*do,j,1,nn
f,nelem(i,j),fz,f_n
*enddo
*endif
esla,s
*enddo
但是把第一段中加一个*if循环,成为第二段就运行不了了:
csys,11
asel,s,area,,1
esla,s
nsla,s,1
*get,enmax,elem,,num,max
*get,enmin,elem,,num,min
fcum,add
*do,i,enmin,enmax,1
*if,esel(i),eq,1,then
*get,ae,elem,i,area
ze=centrz(i)
*if,ze,ge,0,and,ze,lt,depth-0.05,then
p_e=-(ka*gama2*(depth-ze-0.050)+ka*gama1*0.05
*elseif,ze,ge,depth-0.05,and,ze,le,depth,then
p_e=-(ka*gama1*(depth-ze))
*else
p_e=0
*endif
f_tot=p_e*ae
esel,s,elem,,i
nsle,s,corner
*get,nn,node,,count
f_n=f_tot/nn
*do,j,1,nn
f,nelem(i,j),fx,f_n
*enddo
*endif
esla,s
*enddo
加上这段,就会在file.err 中出现 Specified forces from the previous F command were applied only to the
selected NODAL LOAD set.
[ 本帖最后由 p1gx 于 2007-7-17 20:49 编辑 ] |