内力和弯矩——实体模型和梁模型得到的结果不一致，求助

anew2008

!建立长10m，宽1m，高1m的矩形梁，简支结构，受向下的单位均布荷载1kN/m2.
!求指定断面的内力,如支座反力，断面的剪力、弯矩。
/prep7
et,1,solid65  !使用实体单元
mp,ex,1,3.60e10　!弹性模量
mp,prxy,1,0.2　　!泊松比
block,0,1,0,1,0,10 !建立模型
/VIEW,1,1,1,1 　
/pnum,line,1　!线条显示
lplot
!网格划分——长度方向划份100个，宽度和高度划分各10个（每0.1m分为1段，共10000个单元）
LESIZE, 5,,,10
LESIZE, 6,,,10
LESIZE, 9,,,100
MSHKEY,1
vmesh,1
save
/solu
!添加约束，简支结构
/pnum,node,1
nsel,s,loc,y,0
nsel,r,loc,z,0
d,all,ux,0
d,all,uy,0
nsel,s,loc,y,0
nsel,r,loc,z,10
d,all,uy,0
allsel,all
!添加均布荷载。
/pnum,area,1
APLOT
SFA,4,,pres,-1000
solve
finish
save
!后处理，查看内力和弯矩
/post1
nsel,s,loc,z,0
nsum
!应该得到的支座反力为FY＝5000N,弯矩Mx＝0
!ansys计算结果正确：支座反力为FY＝5000N,弯矩Mx约为0

nsel,s,loc,z,5　!跨中断面的内力和弯矩
nsum
!应该得到的断面剪力为FY＝0,弯矩Mx＝12500N/m
!ansys计算结果不对。但是又不知道什么地方不对。

lned2003

是不是单元的实常数没有设置好？

[ 本帖最后由 lned2003 于 2008-9-24 14:15 编辑 ]

ywh362568

实体单元没有实常数设定的，差多少呢，有差距正常的

geni

边界条件不一致所致

幽幽深谷

/post1
nsel,s,loc,z,0
nsum
!应该得到的支座反力为FY＝5000N,弯矩Mx＝0
!ansys计算结果正确：支座反力为FY＝5000N,弯矩Mx约为0
nsel,s,loc,z,5　!跨中断面的内力和弯矩
nsum
!应该得到的断面剪力为FY＝0,弯矩Mx＝12500N/m
!ansys计算结果不对。但是又不知道什么地方不对。

首先，我没有见过nsum的命令，倒是有fsum，用以计算结点的内力和，而不包括外力
nsel,s,loc,z,0选择支座节点，内力和当然等于支座反力
nsel,s,loc,z,5选择中间节点，内力当然为0了，没有弯矩

beyondboy

力学概念不清楚，直接读取节点力其实是所有跟该节点相连单元在该节点处力的总和，而当要得到截面弯矩时，是需要将面分为两边（姑且称为左截面和右截面，如果要计算左侧位置处的弯矩意思是该位置处左侧截面的弯矩），两边截面位置处的弯矩和轴力都是相反数，你得到的结果恰恰证明ansys给你算对了！

xiaoniu520

我好象也明白一点了。直接在后处理中nsel,s,loc,z,5
选出来的节点是一个单独存在的东西，他当然是受力平衡的。所有都为0

支座处的节点由于有支座的附加反力，所以不为0

jeffcs

力学概念不清楚，直接读取节点力其实是所有跟该节点相连单元在该节点处力的总和，而当要得到截面弯矩时，是需要将面分为两边（姑且称为左截面和右截面，如果要计算左侧位置处的弯矩意思是该位置处左侧截面的弯矩）， ...
beyondboy 发表于 2008-9-27 18:42

非常赞同，受教了，呵呵

wszdm0503

不错不错，受教了