函数加载面载荷的问题？

8hao

在一个圆环内表面的250°-343°之间加载随节点角度变化的面载荷,下面是我的命令流，请高手帮忙啊！
!进行物料载荷的面加载
p=2300    !物料的密度
r=2.0       !窑衬的内径
w=1.9       !筒体的转速
a=0.462     !物料安息角
b=0.813     !物料填充角的一半
pi=3.14
g=9.8
D=p*g*r**2
F=0.5*p*g*w**2*r**3
csys,1　　　　　　             ! 进入柱坐标
nsel,all
nsel,s,loc,x,2.0                   ! 选中窑衬内径
nsel,r,loc,y,250,343   
*get,nmax,node,,num,max        !得到所选节点的最大编号，并赋给nmax
*get,nmin,node,,num,min
*dim,t1,array,nmax,1,1
*do,k,nmin,nmax
*if,nsel(k),eq,1,then
c=(ny(k)-270)/180*pi              !取得节点的角度并转换为弧度单位
fn=D*(cos(c-a)-cos(b))*cos(c)/cos(a)+F*((cos(c-a))**2-(cos(b))**2)/((cos(c-a))**2)
t1(k)=fn                         !函数加力
*else
t1(k)=0
*endif
*enddo
sffun,pres,t1(1)
sf,all,pres,0

pfsr

一个类似的例子
http://www.simwe.com/forum/thread-809713-1-1.html

jinpingzhu

money is so less for me to take some time to do this :-)

8hao

原帖由 jinpingzhu 于 2008-3-25 10:56 发表
money is so less for me to take some time to do this :-)

见钱眼开的家伙！

wxwieandy

http://www.simwe.com/forum/thread-809713-1-1.html
很想看看，却打不来，站长能维护一下以前的老帖子么