为什么我看到的梁单元两个节点的弯矩会是一样的？

chenyanru

我作了一个小模型，发现提取弯矩时梁单元的两个不同的点，弯矩为一样的。不知道是计算问题还是数据提取的问题，附命令流，请大家帮我看看。
/prep7
/TITLE,10.5
!********截面1-框架截面**********
ET,1,82
RECTNG,-0.2,0.2,-0.1,0.1
ASEL,ALL
SMRTSIZE,5
AMESH,ALL
SECWRITE,kuangjia,SECT,,1
SECTYPE,1,BEAM,MESH
SECOFFSET,CENT,,,
SECREAD,'kuangjia','SECT',,MESH
ASEL,ALL
ACLEAR,ALL
ADELE,ALL,,,1

!********截面2-柱截面**********

RECTNG,-0.25,0.25,-0.25,0.25
ASEL,ALL
SMRTSIZE,5
AMESH,ALL
SECWRITE,zhu,SECT,,1
SECTYPE,2,BEAM,MESH
SECOFFSET,CENT,,,
SECREAD,'zhu','SECT',,MESH
ASEL,ALL
ACLEAR,ALL
ADELE,ALL,,,1

!******************************
!***定义节点************
n,1,0,0,0
n,2,10,0,0
n,3,0,5,0
n,4,0,10,0
n,5,0,15,0
n,6,5,5,0
n,7,5,10,0
n,8,5,15,0
n,9,10,5,0
n,10,10,10,0
n,11,10,15,0


etdele,1
et,1,beam188
mp,ex,1,3.0e7
mp,dens,1,26
mp,prxy,1,0.1667

/type1
r,1
real,1
mat,1

!****第一阶段
secnum,2
e,1,3
e,2,9
secnum,1
e,3,6
e,6,9
allsel
cm,diyi,elem
esel,none
!****第二阶段
secnum,2
e,3,4
e,9,10
e,6,7
secnum,1
e,4,7
e,7,10
cm,dier,elem
esel,none


!****第三阶段

secnum,2
e,4,5
e,10,11
e,7,8
secnum,1
e,5,8
e,8,11
cm,disan,elem
esel,none


allsel

d,1,all
d,2,all

/solu
acel,,9.8
eqslv,peg

!第一个施工阶段
time,1
nlgeom,on
nropt,full
estif,1e-5
allsel
ekill,all
esel,s,,,diyi
ealive,all
esel,s,live
nsle,s
nsel,inve
d,all,all,0
allsel
nsel,s,loc,y,0
d,all,all
allsel
solve
save
!*********第二阶段******
time,2
allsel
esel,s,,,dier
ealive,all
allsel
ddele,all,all
esel,s,live
nsle,s
nsel,inve
d,all,all,0
allsel
nsel,s,loc,y,0
d,all,all
allsel
solve
save

!*********第三阶段******
time,3
allsel
esel,s,,,disan
ealive,all
allsel
ddele,all,all
esel,s,live
nsle,s
nsel,inve
d,all,all,0
allsel
nsel,s,loc,y,0
d,all,all
allsel
solve
save

/post1
set,1
etable,z11,smisc,3
etable,z12,smisc,16
pretab,z11
pretab,z12

!*dim,diyi1,array,1,11
!*vget,diyi1,elem,1,etab,z11
!*dim,diyi2,array,1,11
!*vget,diyi2,elem,1,etab,z12
!*status,diyi1


set,2
etable,z21,smisc,3
etable,z22,smisc,16
pretab,z21
pretab,z22


set,3
etable,z31,smisc,3
etable,z32,smisc,16
pretab,z31
pretab,z32

chenyanru

是我划分的单元不够，但是现在新的问题出现了！

xinyuxf

新问题是什么？