空心板跨中截面作用单位集中力弯矩不对称

这是命令
finish
/clear
/prep7
et,1,82 ! 划分网路单元类型
et,2,44 ! 空心板
et,3,44 ! 刚性连杆
et,4,44 ! 铰缝
keyopt，4,8,111    !释放44单元的J节点的z,y,x的转动释放   
mp,dens,2,2600
mp,ex,2,3.3e10
mp,prxy,2,0.1667   !定义空心板的材料属性
mp,ex,3,3e16
mp,dens,3,0
mp,prxy,3,0.1667   !定义刚性连杆
mp,ex,4,3e16
mp,dens,4,0
mp,prxy,4,0.1667   
r,3,1e-3,1,1,0.03,0.03
!生成空心板截面并储存文件
rectng,0,0.99,0,0.6
cyl4,0.27,0.26,0.19,-180
cyl4,0.27,0.34,0.19,180
rectng,0.08,0.46,0.26,0.34
/pnum,area,1
/replot    !重新显示图形
agen,2,2,4,1,0.45,0,0,100    !
/replot
asel,s,,,2,7,1
cm,ar,area   !将几何元素分组形成组，以便剪除
allsel
asba,1,ar   !生成空心板
asel,all
esize,0.1
amesh,all
secwrite,kxb
sectype,1,beam,mesh
secoffset,cent   !定义梁节点在形心上
secread,kxb,,,mesh        
                     
asel,all
aclear,all  !将面上的实体部分删除
adele,all,,,1   !清除有所面上的附着的线，关键点
/replot   !这一步之前可能要实现/graphics,full,,/graphics,power-
n,1
n,21,,,12.6
fill,1,21  !填充节点
n,22,0,1,0  !定义参考点
ngen,9,100,1,22,1,1,0,0
type,2     !单元代号
secnum,1
mat,2
*do,i,1,20,1   !循环生成单元
*do,k,1,9,1
e,(k-1)*100+i,(k-1)*100+i+1,22+
*enddo
*enddo
esel,s,,,1,180,1
cm,e2,elem
allsel
type,4
mat,3
r,3
*do,i,1,21,1  !循环生成连杆单元
*do,k,1,8,1
e,(k-1)*100+i,k*100+i
*enddo
*enddo
esel,s,,,181,348,1
cm,e4,elem
!  选择梁的端点进行约束
*do,i,0,8,1
allsel
nsel,s,loc,z,0
nsel,r,loc,y,0
nsel,r,loc,x,i
d,all,ux,,,,,uy,uz,rotz,roty
allsel
nsel,s,loc,z,12.6
nsel,r,loc,y,0
nsel,r,loc,x,i
d,all,ux,,,,,uy,rotz,roty
allsel
*enddo
!!!开始加载单位集中力
/solu
*do,i,0,8,1
time,i+1
fdele,all,all
nsel,s,loc,z,6.3
nsel,r,loc,y,0
nsel,r,loc,x,i
f,all,fy,-1
allsel
solve
*enddo
/post1
set,1
etable,z1,smisc,12   !J端得弯矩
plls,z1,z1

请高手指教，应该不是约束的问题吧！

9块空心板，跨径12。6m。我查看单元的I,J两端弯矩是对称的！为什么查看J端得时候就不对称啊，这是桥工（范立础）上的例题，请高手帮忙！

不同单元在共同节点处存在应力差，因为形函数的假定是人为的。所以求得的弯矩也也是不同的，因此常常需要应力磨平，所以选弯矩表时也要选择I，J两端，然后程序自动将弯矩平均得到。

加载有问题，你看你的*ENDDO位置放错了，还在你只在梁边上加了一个偏载，这是你的图

还有就是你的虚拟横梁刚度也不用那么大吧？
都比空心板还大了？你看那个载面尺寸

看看，不过我是不想运行这陀命令流