节点耦合自由度释放的实现

本人对节点耦合情况不是很清楚，从王新敏的ansys书上找到的例子ex6.4，ex6.5（P352）能得到一些节点耦合的命令流分析。但是本人对ex6.5中第二个耦合节点的分析有一些疑问（节点编号见图）

怎么实现第二个耦合节点的相应图示自由度释放？
!王新敏ansys教材中的例子
!EX6.5 BEAM189自由度释放
finish
/clear
/prep7
et,1,beam189
mp,ex,1,2.1e11
mp,prxy,1,0.3
sectype,1,beam,csolid
secdata,0.05
k,1
k,2,0,4
k,3,3,0
k,4,3,4
k,5,3,7
k,6,6,4
k,7,6,7
k,8,9,7
k,100,-10,110
l,1,2
l,3,4
l,4,5
l,2,4
l,6,7
l,5,7
l,7,8
dk,1,all
dk,3,all
dk,6,all
dk,8,ux,,,,uy
latt,1,,1,,100,,1
lesize,all,,,9
lmesh,all
lsel,s,loc,y,4
esll
sfbeam,all,1,pres,2000
lsel,s,loc,y,7
esll
sfbeam,all,1,pres,3000
allsel

lsel,s,,,1,4,3
esll
endrelease,,,ux

lsel,s,,,2,4,2
esll
endrelease,,,rotz

lsel,s,,,5,7
esll
endrelease,,,rotz
allsel
finish
/solu
solve
finish
/post1
pldisp,1
etable,mi,smisc,2
etable,mj,smisc,15
plls,mi,mj

这个好像不难吧，铰节点属于那个单元就把那个单元相应的转角自由度释放就好。

ls不妨试试，贴出命令流吧。谢谢啦~