[求助]哪位大侠能帮我解决下图的速度梯度的问题！

johnp

[HTML]<IMG align=baseline alt="" border=0 hspace=0 src="D:\Documents and Settings\pjh1\My Documents\My Pictures\12.JPG">[/HTML]

johnp

我现在做的是螺槽是物料的流动分析，大家看到的是螺槽中物料展开的情况，我想在图中红色的面上，加有梯度（梯度沿高度方向）的速度（速度方向沿长度方向），因为螺槽中物料是旋转运动，沿半径方向有速度梯度（红色侧面）。我用的是ANSYS5.7，不知哪位大哥能帮我指点一下，我已经被它折磨了N天了，可怜！
谢了！

redhong731

你把图弄清数些，最好压缩上传，。我看不清！

johnp

版主大哥！你好！我把我问题的主要图形都贴上来吧，这是一个螺杆，螺杆旋转工作带动物料最终产生轴向运动，物料得到的速度有螺杆轴毂面传递的速度和螺棱给物料侧面传递的速度，螺棱的旋转运动，其给物料传递的切向速度沿半径方向肯定是有梯度，后面的两副图是一副是物料原始形状C形，一副是简化的图形长方体。

johnp

这副图是物料在螺槽中的原始形状，其运动状态是绕Z轴旋转，C形侧面的切向速度沿半径方向有梯度。

johnp

这副图是物料的简化图，原来的C形侧面就变为简化图的长的高度方向的侧面，原来的切向速度就转化为X方向速度，原来的径向就转化为Y方向，因此简化图中的X方向的速度沿Y向存在梯度，小弟我不知这个速度梯度如何加，请大哥帮帮我，谢谢了！

johnp

我经过一个多星期的摸索终于弄出来了，当然了也少不了流体中文网以前先驱们发过的挨到一点边的贴子啦，谢谢了!APDL传上，供大家讨论
!* 螺棱面上的速度约束
ASEL,S, , ,2
NSLA,S,1
*get,tolmax2,node,,count
*dim,nnum2,,tolmax2,3, !dim an array to restore node
*get,maxnum2,node,,num,max !get the max node number
nnum2(1,1)=maxnum2
nnum2(1,2)=(ny(nnum2(1,1))+rs)*w*cos(sa)
nnum2(1,3)=(ny(nnum2(1,1))+rs)*w*sin(sa)
D,nnum2(1,1),Vx,-nnum2(1,2)
D,nnum2(1,1),Vz,nnum2(1,3)
*do,j2,2,tolmax2
*get,nnum2(j2,1),node,nnum2(j2-1,1),nxtl,
nnum2(j2,2)=(ny(nnum2(j2,1))+rs)*w*cos(sa) !the velocity of x direction
nnum2(j2,3)=(ny(nnum2(j2,1))+rs)*w*sin(sa) !the velocity of z direction
D,nnum2(j2,1),Vx,-nnum2(j2,2)
D,nnum2(j2,1),Vz,nnum2(j2,3)
*enddo
ASEL,S, , ,4
NSLA,S,1
*get,tolmax4,node,,count
*dim,nnum4,,tolmax4,3, !dim an array to restore node
*get,maxnum4,node,,num,max!get the max node number
nnum4(1,1)=maxnum4
nnum4(1,2)=(ny(nnum4(1,1))+rs)*w*cos(sa)
nnum4(1,3)=(ny(nnum4(1,1))+rs)*w*sin(sa)
D,nnum4(1,1),Vx,-nnum4(1,2)
D,nnum4(1,1),Vz,nnum4(1,3)
*do,i4,2,tolmax4
*get,nnum4(i4,1),node,nnum4(i4-1,1),nxtl,
nnum4(i4,2)=(ny(nnum4(i4,1))+rs)*w*cos(sa) !the velocity of x direction
nnum4(i4,3)=(ny(nnum4(i4,1))+rs)*w*sin(sa) !the velocity of z direction
D,nnum4(i4,1),Vx,nnum4(i4,2)
D,nnum4(i4,1),Vz,nnum4(i4,3)
*enddo
allsel

joli

呵呵,兄弟是那个学校的,是不是搞聚合物加工的