ansys中内力图正负号与结构力学中的内力图正负号之探讨

下面是结构力学中的例题

利用ansys求解该例题的命令流如下：
FINISH
/CLEAR
! 1 在前处理模型中建模
/PREP7          ! 进入前处理模块Prep7
N,1,0,0,0       ! 定义节点1位置(0,1,0)
N,2,0,5,0       ! 定义节点2位置(1,1,0)
N,3,2.5,5,0       ! 定义节点3位置(1,0,0)
N,4,5,5,0     ! 定义节点4位置(0.6,1,0)
ET,1,beam3              ! 定义单元为二维梁单元beam3
R,1,0.5,1/24,   ! 定义单元实常数--截面积，抗弯惯性矩和截面高度
MP,EX,1,3e10            ! 定义材料弹性模量EX=2e11 N/m^2
E,1,2       ! 定义1,4节点组成的单元1
E,2,3       ! 定义4,2节点组成的单元2
E,3,4      ! 定义2,5节点组成的单元3

! 2 定义位移约束条件和荷载
D,1,ALL     ! 固定1号节点的所有自由度
D,4,ALL     ! 固定3号节点的所有自由度
F,2,FX,-10000     ! 2号节点沿着X轴方向的集中力100N（向右）
F,2,FY,-6000    ! 2号节点沿着-Y轴方向的集中力300N（向下）
F,2,MZ,-10000     ! 2号节点绕着Z轴顺时针转向的集中力偶50N.m
F,3,Fy,-8000     ! 2号节点沿着X轴方向的集中力100N（向右）
SFBEAM,1,1,PRESS,12000     ! 1号单元上沿着局部坐标系的-Y轴方向的均布荷载800N/m（向下）
FINISH    ! 退出前处理模块
! 3 在求解模块求解
/SOLU      ! 进入求解模块
SOLVE      ! 开始求解
FINISH     ! 退出求解模块
! 4 在后处理模块显示计算结果
/POST1     ! 进入后处理模块
PLDISP,2   ! 带原有结构轮廓线显示变形图
PRDISP      ! 列表显示节点位移计算结果
! 建立单元表
ETABLE,NI,SMISC,1     !单元I点轴力
ETABLE,NJ,SMISC,7     !单元J点轴力
ETABLE,QI,SMISC,2     !单元I点剪力
ETABLE,QJ,SMISC,8     !单元J点剪力
ETABLE,MI,SMISC,6     !单元I点弯矩
ETABLE,MJ,SMISC,12 !单元J点弯矩
etable,smax1,nmisc,1
etable,smax2,nmisc,3
etable,smin1,nmisc,2
etable,smin2,nmisc,4
PRETAB              ! 显示单元表结果
! 更新单元表
ETABLE,REFL
! 画轴力分布图
/TITLE,Axial   force   diagram  ! 定义图形窗口标题
PLLS,NI,NJ,1.0,0    !ansys默认设置
/image,save,'Axial_force_%T%',jpg
! 画剪力分布图
/TITLE,Shearing   force   diagram
PLLS,QI,QJ,1.0,0!ansys默认设置
/image,save,'Shearing_force_%T%',jpg
! 画弯矩分布图
/TITLE,Bending   moment diagram
PLLS,MI,MJ,1,0!ansys默认设置
/image,save,'Bending_moment_%T%',jpg
plls,smax1,smax2
plls,smin1,smin2


                                                    弯矩图
    说明：由于命令“PLLS”只使用了单元两端的内力，所以只能绘制出内力线性变化情况〔这对于只有节点载荷作用情况是准确的)。如果单元中间有分布载荷，则单元中的内力不再是线性变化的时候，该内力图在单元中间是不真实的，这时可以在分布载荷所在的结构上多划分些单元来实现。

                                                     剪力图

                                                     轴力图
由此可以看出:
在ansys默认情况下，轴力图的方向与结构力学中轴力图方向一致，弯矩和剪力的方向与结构力学中规定的方向相反。

节点力和节点位移在ANSYS中的转向是按Z坐标轴的指向来规定的，顺时针转向正好是Z轴的负方向，结果为负值。
ANSYS中的正剪力是使单元体产生绕Z轴产生正方向转动(逆时针)为正。


若想让ansys中内力正负号与结构力学中相统一，只需将后处理做如下修改：

! 画轴力分布图
/TITLE,Axial   force   diagram  ! 定义图形窗口标题
PLLS,NI,NJ,1.0,0    !ansys默认设置
/image,save,'Axial_force_%T%',jpg
! 画剪力分布图
/TITLE,Shearing   force   diagram
PLLS,QI,QJ,-1.0,0
/image,save,'Shearing_force_%T%',jpg
! 画弯矩分布图
/TITLE,Bending   moment diagram
PLLS,MI,MJ,-1,0







分布载荷所在的结构上多划分些单元的命令流如下：
FINISH
/CLEAR
! 1 在前处理模型中建模
/PREP7          ! 进入前处理模块Prep7
*do,i,1,11   
N,i,0,0.5*(i-1),
*enddo  
N,12,2.5,5,0       ! 定义节点3位置(1,0,0)
N,13,5,5,0     ! 定义节点4位置(0.6,1,0)
ET,1,beam3              ! 定义单元为二维梁单元beam3
R,1,0.5,1/24,   ! 定义单元实常数--截面积，抗弯惯性矩和截面高度
MP,EX,1,3e10            ! 定义材料弹性模量EX=2e11 N/m^2
*do,i,1,12
e,i,i+1,
*enddo
! 2 定义位移约束条件和荷载
D,1,ALL     ! 固定1号节点的所有自由度
D,13,ALL     ! 固定3号节点的所有自由度
F,11,FX,-10000     ! 2号节点沿着X轴方向的集中力100N（向右）
F,11,FY,-6000    ! 2号节点沿着-Y轴方向的集中力300N（向下）
F,11,MZ,-10000     ! 2号节点绕着Z轴顺时针转向的集中力偶50N.m
F,12,Fy,-8000     ! 2号节点沿着X轴方向的集中力100N（向右）
esel,s,,,1,10
SFBEAM,all,1,PRESS,12000     ! 1号单元上沿着局部坐标系的-Y轴方向的均布荷载800N/m（向下）
allsel
eplot
FINISH    ! 退出前处理模块
! 3 在求解模块求解
/SOLU      ! 进入求解模块
SOLVE      ! 开始求解
FINISH     ! 退出求解模块
! 4 在后处理模块显示计算结果
/POST1     ! 进入后处理模块
PLDISP,2   ! 带原有结构轮廓线显示变形图
PRDISP      ! 列表显示节点位移计算结果
! 建立单元表
ETABLE,NI,SMISC,1     !单元I点轴力
ETABLE,NJ,SMISC,7     !单元J点轴力
ETABLE,QI,SMISC,2     !单元I点剪力
ETABLE,QJ,SMISC,8     !单元J点剪力
ETABLE,MI,SMISC,6     !单元I点弯矩
ETABLE,MJ,SMISC,12 !单元J点弯矩
etable,smax1,nmisc,1
etable,smax2,nmisc,3
etable,smin1,nmisc,2
etable,smin2,nmisc,4
PRETAB              ! 显示单元表结果
! 更新单元表
ETABLE,REFL
! 画轴力分布图
/TITLE,Axial   force   diagram  ! 定义图形窗口标题
PLLS,NI,NJ,1.0,0   
/image,save,'Axial_force_%T%',jpg
! 画剪力分布图
/TITLE,Shearing   force   diagram
PLLS,QI,QJ,1.0,0
/image,save,'Shearing_force_%T%',jpg
! 画弯矩分布图
/TITLE,Bending   moment diagram
PLLS,MI,MJ,1,0
/image,save,'Bending_moment_%T%',jpg
plls,smax1,smax2
plls,smin1,smin2

楼主知道如何加倾斜梁单元吗？麻烦上传一下

强悍 我喜欢  试学了一下 谢谢 以后楼主多多发一些啊

支持你lz 非常喜欢 以后一定要多发点关于此类的好东西

好厉害啊  膜拜