ansys中内力图正负号与结构力学中的内力图正负号之探讨
本帖最后由 西风独自凉 于 2011-3-28 10:29 编辑下面是结构力学中的例题
利用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
楼主知道如何加倾斜梁单元吗？麻烦上传一下 强悍 我喜欢试学了一下 谢谢 以后楼主多多发一些啊{:4_126:} 支持你lz 非常喜欢 以后一定要多发点关于此类的好东西 好厉害啊膜拜
页:
[1]