sectype,1,beam,i
secdata,16.655,16.655,36.74,1.68,1.68,0.945
secplot,1
c=1.49535
sectype,2,beam,i
secdata,c*16.655,c*16.655,c*36.74,c*1.68,c*1.68,c*0.945
secplot,2
a=400
K,1
K,2,,,a
K,3,2*a
K,4,2*a,,a
L,1,2
L,4,3
L,2,4
lesize,1,,,4
lesize,2,,,4
lesize,3,,,16
lsel,,,,1
latt,1,1,1,,,,1
lmesh,all
lsel,,,,2
latt,1,1,1,,,,2
lmesh,all
lsel,,,,3
latt,1,1,1,,,,2
lmesh,all
allsel
/solu
dk,1,all
dk,3,all
lsel,,,,3
esll
sfbeam,all,1,pres,500,500
allsel
solve
/post1
lsel,,,,1
esll !选择编号为1—4的单元,节点编号分别为1,3,4,5,2
presol,m !list elem solution(list单元解的弯矩)
prnld,m !list node solution(list节点解的弯矩)
etable,myi,smisc,2
etable,myj,smisc,15
pretab,myi,myj !先定制单元表,后list单元表的结果
三种命令得到的结果如下(MY):
1.单元解
ELEM= 1 MX MY MZ
1 -0.20200E-05-0.72770E+07-0.56976E-09
3 0.20197E-05 0.24390E+07 0.56976E-09
ELEM= 2 MX MY MZ
3 -0.20192E-05-0.24390E+07-0.56976E-09
4 0.20190E-05-0.23990E+07 0.56976E-09
ELEM= 3 MX MY MZ
4 -0.20184E-05 0.23990E+07-0.21892E-07
5 0.20182E-05-0.72371E+07 0.21892E-07
ELEM= 4 MX MY MZ
5 -0.20177E-05 0.72371E+07-0.43301E-07
2 0.20174E-05-0.12075E+08 0.43301E-07
2.节点解
NODE MX MY MZ
1 -0.20200E-05-0.72770E+07-0.56976E-09
2 0.20174E-05-0.12075E+08 0.43301E-07
3 0.52574E-09
4 0.52574E-09 -0.21322E-07
5 0.52574E-09 -0.21409E-07
3.etable解
ELEM MYI MYJ
1 -0.48580E+07-0.48580E+07
2 -19965. -19965.
3 0.48181E+07 0.48181E+07
4 0.96561E+07 0.96561E+07
三种解中的符号这样规定:若弯矩轴(在这里是Y轴)的正向指向屏幕内,则弯矩以绕隔离体逆时针旋转为正;若若弯矩轴(在这里是Y轴)的正向指向屏幕外,则弯矩以绕隔离体顺时针旋转为正。
单元解和节点解在首尾节点1,2处的MY值是一样的,节点解在其他节点的MY值为零,这是因为节点解对中间节点做了叠加。
etable解中MY中的Y指的是单元坐标系的Y,这里单元坐标系的Y刚好与总体坐标系的Y同向。etable每个单元中的I、J端的值相等,且其绝对值等于单元解中每个单元I、J绝对值和的平均。这是由188单元的计算方法决定的。
一般来说list结果用单元解,而画弯矩图用etable。