surf154加载问题

16443

想通过surf154单元，对一表面施加3个方向的面压，即法向和两个切向。
commands出入如下命令，为什么总是加载不上那，请各位给指点以下错误

cmsel,s,surf1                               
esel,s,ename,,154  !select elem with Ename
*GET, maxelem, ELEM, 0, NUM,MAX                       
*GET, Par_type, ELEM, maxelem, ATTR, type
*dim,p1,,maxelem
*dim,p2,,maxelem
*dim,p3,,maxelem

/prep7
keyopt,Par_type,2,1           !Pressure loads according to the local coordinate system
FINISH
!*
/solu
*afun,deg                               
csys,12                   !select local coordinate system
esys,12
elnum=0
*DO,i,1,maxelem
  elnum=elnext(elnum)
  *GET, locx, ELEM, elnum, CENT, X                 !get elem centroid location
  *GET, locy, ELEM, elnum, CENT, Y
  *GET, locz, ELEM, elnum, CENT, Z
        lambda1=lambdap1-u1+(locy*u6)-(locz*u5)   
        lambda2=lambdap2-u2-(locz*u4)+(locx*u6)  
        lambda3=lambdap3-u3+(locx*u5)-(locy*u4)  
        *if,lambda2,GE,0,THEN
                sign=1
        *else
                sign=0
        *endif
        p1(i)=alphat*(lambda2**betat)*lambda1*sign       
        p2(i)=alphan*(lambda2**betan)*sign           
        p3(i)=alphat*(lambda2**betat)*lambda3*sign   
       
            SFE,elnum,1,PRESS,,-p1,-p2,-p3   ！是不是这里有问题呢 ？
  *ENDDO

cumtmbqq

surf154加载你说的那个力应该加三次吧！
我对命令流不了解，是个新手。
surf154  的 LEKY 1：-Z方向
                2：tanx
                            3：tany

16443

2# cumtmbqq
谢谢提醒。

SFE, ELEM, LKEY, Lab, KVAL, VAL1, VAL2, VAL3, VAL4中，lkey表示154（153）单元的面号（如face1，face2，。。。face5）；val1~4 含义各不相同，总的意思就是载荷的大小。

SFE,elnum,1,PRESS,,-p1,-p2,-p3   ！是不是这里有问题呢 ？

这个应该是:
sfe,elnum,1,PRESS,,-p1
sfe,elnum,2,PRESS,,-p2
sfe,elnum,3,PRESS,,-p3

cumtmbqq

3# 16443
呵呵，楼主，你学命令流的过程能和我分享下不？
谢谢！！