正在加载...
建立了一个受到端部集中荷载的悬臂梁，用面操作方法求解截面内力。梁的中部各个截面所得结果与理论值非常接近，但在梁的两端（自由端和固定端）内力与理论值出现了非常大的不同。附上代码，请大家帮忙看看！
申明：内容来自用户上传，著作权归原作者所有，如涉及侵权问题，请点击此处联系，我们将及时处理！
全部回答（1 ）
上面的附件有点问题，在施加节点荷载之前没有选中节点。我重新贴一下吧。
/CLEAR,START/prep7
ET,1,SOLID185
MP,EX,1,3.0E10 !材料1弹模MP,DENS,1,1800 !材料1密度MP,PRXY,1,0.25 !材料1泊松比
BLOCK,0,2,0,2,0,40
VSEL,ALLVATT,1,,1,0
MSHAPE,0,3DMSHKEY,1 !选择映射划分方式
LSEL,ALLESIZE,0.1
VSEL,ALLVMESH,ALLFINISH
/solu !进入求解
!底部固定NSEL,ALLNSEL,s,loc,z,0D,all,all,0!顶部施加水平集中力NSEL,ALLF,463,FX,5000F,443,FX,5000
allsel !选择所有solvefinish
/post1*do,ii,1,21WPCSYS,-1WPOFFS,1,1,(ii-1)*2SUCR,suz%ii%,cplane !面的建立
SUMAP,mysx,s,xSUMAP,mysy,s,ySUMAP,mysz,s,zSUMAP,mysxy,s,xySUMAP,mysyz,s,yzSUMAP,mysxz,s,xz
SUEVAL,xforce,mysxz,intg !求截面上FXSUEVAL,yforce,mysyz,intg !求截面上FYSUEVAL,zforce,mysz,intg !求截面上FZ
SUEVAL,mya,da,sum !求截面面积并赋给变量myaSUEVAL,myya,gcy,intg !求对X轴的面积矩并赋给变量myyamyya=myya/mya !得到面积重心到X轴的距离=面积矩/面积SUEVAL,myxa,gcx,intg !求对Y轴的面积矩并赋给变量myxamyxa=myxa/mya !得到面积重心到Y轴的距离=面积矩/面积
SUCALC,szgcy,mysz,mult,gcy !计算mysz乘以gcy，并赋给面项szgcySUEVAL,mx1,szgcy,intg !对面项szgcy在面上积分，得到X方向弯矩mx1
SUCALC,szgcx,mysz,mult,gcx !计算mysz乘以gcx，并赋给面项szgcxSUEVAL,my1,szgcx,intg !对面项szgcx在面上积分，得到Y方向弯矩my1
*cfopen,result%ii%,txt,,APPEND !输出计算结果*VWRITE,xforce,yforce,zforce,mx1,my1(5E16.8)
SUDEL,suz%ii%
*enddo
主题
回复
粉丝
混凝土结构
140.21 万条内容 · 2162 人订阅
遇到了一点小问题
内容不能包含下列词语
回帖成功
经验值 +10
全部回答（1 ）
上面的附件有点问题，在施加节点荷载之前没有选中节点。我重新贴一下吧。
/CLEAR,START
/prep7
ET,1,SOLID185
MP,EX,1,3.0E10 !材料1弹模
MP,DENS,1,1800 !材料1密度
MP,PRXY,1,0.25 !材料1泊松比
BLOCK,0,2,0,2,0,40
VSEL,ALL
VATT,1,,1,0
MSHAPE,0,3D
MSHKEY,1 !选择映射划分方式
LSEL,ALL
ESIZE,0.1
VSEL,ALL
VMESH,ALL
FINISH
/solu !进入求解
!底部固定
NSEL,ALL
NSEL,s,loc,z,0
D,all,all,0
!顶部施加水平集中力
NSEL,ALL
F,463,FX,5000
F,443,FX,5000
allsel !选择所有
solve
finish
/post1
*do,ii,1,21
WPCSYS,-1
WPOFFS,1,1,(ii-1)*2
SUCR,suz%ii%,cplane !面的建立
SUMAP,mysx,s,x
SUMAP,mysy,s,y
SUMAP,mysz,s,z
SUMAP,mysxy,s,xy
SUMAP,mysyz,s,yz
SUMAP,mysxz,s,xz
SUEVAL,xforce,mysxz,intg !求截面上FX
SUEVAL,yforce,mysyz,intg !求截面上FY
SUEVAL,zforce,mysz,intg !求截面上FZ
SUEVAL,mya,da,sum !求截面面积并赋给变量mya
SUEVAL,myya,gcy,intg !求对X轴的面积矩并赋给变量myya
myya=myya/mya !得到面积重心到X轴的距离=面积矩/面积
SUEVAL,myxa,gcx,intg !求对Y轴的面积矩并赋给变量myxa
myxa=myxa/mya !得到面积重心到Y轴的距离=面积矩/面积
SUCALC,szgcy,mysz,mult,gcy !计算mysz乘以gcy，并赋给面项szgcy
SUEVAL,mx1,szgcy,intg !对面项szgcy在面上积分，得到X方向弯矩mx1
SUCALC,szgcx,mysz,mult,gcx !计算mysz乘以gcx，并赋给面项szgcx
SUEVAL,my1,szgcx,intg !对面项szgcx在面上积分，得到Y方向弯矩my1
*cfopen,result%ii%,txt,,APPEND !输出计算结果
*VWRITE,xforce,yforce,zforce,mx1,my1
(5E16.8)
SUDEL,suz%ii%
*enddo