查轴力：首先定义单元表grneralpostproc>element table >define table add 左侧选by sequence num,右侧选择smisc, 在下面输入smisc,1 然后在plot results>contour plot》line elem res 查看
弯矩
1.绘制弯矩图
建立弯矩单元表。例如梁单元
i节点单元表名称为imom，j节点单元表名称为jmom，
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点弯矩
plls,MI,MJ
2.标注弯矩图
PLOTCTRLS>>NUMBERING>>SVAL ON即可在画出弯矩图的同时在图上标出弯矩值的大小
3.调整弯矩图
如果弯矩图方向错误，则绘制弯矩图命令为
plls,imom,jmom，-1
同一个节点处两边的单元内力有细微差别，
导致内力数字标注出现重影。观察上面整体轴力图也可以发现，
一段一段的，好像马赛克，其实上面整体弯矩图也是，不过不是
很明显罢了。这是EULER-BEONOULI梁理论以及ANSYS输出定义造成
的（详细原因就不展开了，看看梁理论的书和ANSYS的说明吧）。
为了修正重影和节点两边内力值不一样的问题，遍制了宏文件ITFAVG.MAC
命令文件内容如下：
!---------------------------------------------------------------------
!宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO）
!获取线性单元内力，并对单元边界处的内力进行平衡
!输入信息
!内力类型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ
*ASK,ITFTYPE,'PLEASE INPUT THE TYPE OF INTERNAL FORCE','MMOMY'
!需处理的单元包
*ASK,EASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF ELEMENTS TO BE PROCESSED!',
'EOUTER'
!需处理的节点包
*ASK,NASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF NODE TO BE PROCESSED!','NOU TER'
!无需处理的节点包
*ASK,UNASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NODE!(NONE I F THERE'S NO SUCH COMPONENT)','NONE'
/POST1
!输入信息：内力类型，欲处理单元的集合，欲处理节点的集合
!ITFTYPE='MMOMY'
!EASSEMBLY='EOUTER'
!NASSEMBLY='NOUTER'
!按内力类型确定ANSYS输出信息SMISC的编号
*IF,ITFTYPE,EQ,'MFORX',THEN
ITFINUM=1
ITFJNUM=7
*ELSEIF,ITFTYPE,EQ,'MFORY',THEN
ITFINUM=2
ITFJNUM=8
*ELSEIF,ITFTYPE,EQ,'MFORZ',THEN
ITFINUM=3
ITFJNUM=9
*ELSEIF,ITFTYPE,EQ,'MMOMX',THEN
ITFINUM=4
ITFJNUM=10
*ELSEIF,ITFTYPE,EQ,'MMOMY',THEN
ITFINUM=5
ITFJNUM=11
*ELSEIF,ITFTYPE,EQ,'MMOMZ',THEN
ITFINUM=6
ITFJNUM=12
*ELSE
*ENDIF
!对不需平均的节点进行处理
*IF,UNASSEMBLY,NE,'NONE',THEN
!选出不进行处理的节点包并获取不进行处理节点的数目
CMSEL,S,UNASSEMBLY
*GET,UNNODNUM,NODE,0,COUNT
!定义长度为UNNODNUM的数组（UNNOD），以存放选中单元的单元编号
*DIM,UNNOD,ARRAY,UNNODNUM
!将选中单元的编号按顺序存入数组UNNOD
*DO,I,0,UNNODNUM-1,1
UNNOD(I+1)=NDNEXT(I)
*ENDDO
UNNODNUM=0
*ENDIF
!选出所需的单元和节点包
CMSEL,S,EASSEMBLY
CMSEL,S,NASSEMBLY
!获得当前选中单元总数（存入变量SELELENUM）
*GET,SELELENUM,ELEM,0,COUNT
!定义长度为SELELENUM的数组（ELENUM），以存放选中单元的单元编号
*DIM,ELENUM,ARRAY,SELELENUM
!将选中单元的编号按顺序存入数组ELENUM
*DO,I,0,SELELENUM-1,1
ELENUM(I+1)=ELNEXT(I)
*ENDDO
!获得当前选中节点总数（存入变量SELNODNUM）
*GET,SELNODNUM,NODE,0,COUNT
!定义长度为SELNODNUM的数组（NODNUM），以存放选中单元的单元编号*DIM,NODNUM,ARRAY,SELNODNUM
!将选中单元的编号按顺序存入数组NODNUM
*DO,I,0,SELNODNUM-1,1
NODNUM(I+1)=NDNEXT(I)
*ENDDO
!定义所需的线性单元内力ETABLE，节点I的内力存入数组ITNFI,
!节点J的内力存入数组ITNFJ
ETABLE,ITNFI,SMISC,ITFINUM
ETABLE,ITNFJ,SMISC,ITFJNUM
!定义所需的结果数组，并将其置零
ETABLE,ITNFINEO,SMISC,5
SADD,ITNFINEO,ITNFI,,1
ETABLE,ITNFJNEO,SMISC,11
SADD,ITNFJNEO,ITNFJ,,1
*DO,K,1,SELNODNUM,1
!处理不需平均的节点
INDEX=0
*IF,UNNODNUM,GE,1,THEN
*DO,J,1,UNNODNUM
*IF,NODNUM(K),EQ,UNNOD(J),THEN
INDEX=1
*ELSE
*ENDIF
*ENDDO
*ELSE
*DO,J,1,SELELENUM,1
!选出和节点K相连的线性单元中，I节点（对线性单元而言）为节点K的单元编号
*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THEN
ELEI=ELENUM(J)
*EXIT
*ELSE
*ENDIF
*ENDDO
*DO,J,1,SELELENUM,1
!选出和节点K相连的线性单元中，J节点（对线性单元而言）为节点K的单元编号
*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THEN
ELEJ=ELENUM(J)
*EXIT
*ELSE
*ENDIF
*ENDDO
*IF,INDEX,EQ,0,THEN
*IF,ELEJ,NE,0,THEN !有可能出现ELEJ为0的情况
!取出I节点为节点K的单元的I节点端的内力放入参数ETELEI
*GET,ETELEI,ELEM,ELEI,ETAB,ITNFI
!取出J节点为节点K的单元的J节点端的内力放入参数ETELEJ
*GET,ETELEJ,ELEM,ELEJ,ETAB,ITNFJ
!平均节点K的单元的I节点端的内力和节点K的单元的J节点端的内力
ETAVE=(ETELEI+ETELEJ)/2
!将平均后的内力存入结果数组中
DETAB,ELEI,ITNFINEO,ETAVE
DETAB,ELEJ,ITNFJNEO,ETAVE
*ELSE
*ENDIF
*ELSE
*ENDIF
*ENDDO
/UDOC,1,LOGO,OFF
PLLS,ITNFINEO,ITNFJNEO
!END OF ITFAVG.MAC
（2）
对体和面来说，ANSYS默认的结果输出格式是云图格式，而这种彩色云图打印为黑白图像时对比很不明显，无法表达清楚，对于发表文章非常不便。发文章所用的结果图最好是等值线图，并且最好是黑白的等值线图。一般借用photoshop等第三方软件，很麻烦，效果也不好。ansys本身也能实现这项功能。
1.将要输出的结果调出，这时为彩色云图；
2.将云图转换为等值线图的形式
GUI：plotCtrls—>Device Options—>[/DEVI]中的vector mode 选为on
命令：/DEVICE,VECTOR,1
这时结果为彩色等值线，若直接输出，打印为黑白图像时仍然不清晰，为此需进行以下几步将图像转换为黑白形式；
3. 将背景变为白色
命令：jpgprf,500,100,1
/rep
4.对等值线中的等值线符号（图中为A,B,C等）的疏密进行调整
GUI：plotCtrls—>Style—>Contours—> Contours Labeling 在Key Vector mode contour label 中选中on every Nth elem，然后在N= 输入框中输入合适的数值，例如5，多试几次，直到疏密合适
命令：/clabel,1,5
5.将彩色等值线变为黑色
GUI：plotCtrls—>Style—>Colors—>Contours Colors 将Items Numbered 1，Items Numbered 2等复选框中的颜色均选为黑色，图像即可变为黑白等值线图像
命令：/color,cntr,whit,1 等等
6.修改
在显示等值线时，图形中太多的标识字母（A、B等）使得等值线图完全看不清楚。请问有没有什么方法可以去掉图形的这些标识字母呢
在plotCtrls下面的style下面的contours下面的contour Labeling，在弹出的contour labeling options里面的vector mode contour labels选off
后处理中画出了结构的等值线图，可否将等值线边的字母改为相应的数值标注
/PNUM,SVAL,1
/replot
如果还想将等值线的字母去掉
/PNUM,SVAL,1
/replot
如果还想将等值线的字母去掉
/PNUM,SVAL,1
/CLABEL,1,-1
/replot
7.出图
GUI：plotCtrls—>Capture Image
! 建立单元表
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,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,-0.8,0
/image,save,'Bending_moment_%T%',jp
查轴力：首先定义单元表grneralpostproc>element table >define table add 左侧选by sequence num,右侧选择smisc, 在下面输入smisc,1 然后在plot results>contour plot》line elem res 查看 弯矩 1.绘制弯矩图 建立弯矩单元表。例如梁单元 i节点单元表名称为imom，j节点单元表名称为jmom， 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点弯矩 plls,MI,MJ 2.标注弯矩图 PLOTCTRLS>>NUMBERING>>SVAL ON即可在画出弯矩图的同时在图上标出弯矩值的大小 3.调整弯矩图 如果弯矩图方向错误，则绘制弯矩图命令为 plls,imom,jmom，-1 同一个节点处两边的单元内力有细微差别， 导致内力数字标注出现重影。观察上面整体轴力图也可以发现， 一段一段的，好像马赛克，其实上面整体弯矩图也是，不过不是 很明显罢了。这是EULER-BEONOULI梁理论以及ANSYS输出定义造成 的（详细原因就不展开了，看看梁理论的书和ANSYS的说明吧）。 为了修正重影和节点两边内力值不一样的问题，遍制了宏文件ITFAVG.MAC 命令文件内容如下： !--------------------------------------------------------------------- !宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO） !获取线性单元内力，并对单元边界处的内力进行平衡 !输入信息 !内力类型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ *ASK,ITFTYPE,'PLEASE INPUT THE TYPE OF INTERNAL FORCE','MMOMY' !需处理的单元包 *ASK,EASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF ELEMENTS TO BE PROCESSED!', 'EOUTER' !需处理的节点包
实体单元弯矩和轴力的提取 对于一个实体梁 长10m，宽0.6m，高1.0m 程序代码 finish$/clear finish$/prep7 blc4,0,0,10,0.6,1 ET,1,SOLID45$MP,EX,1,200e6$MP,PRXY,1,0.3!定义单元类型、材料特性 ESIZE,0.5$VMESH,ALL FINISH$/SOLU asel,s,loc,x,0 nsla,s,1$d,all,all nsel,s,loc,x,10 *get,Nnode,node,0,count f,all,fx,100/Nnode
f,all,fy,30/Nnode f,all,fz,50/Nnode SOLVE finish$/post1 ! define surface wpcsys,-1 wpoffs,5,0.3,0.5$wprota,0,0,90 sucr,xsurf1,cplane,3$!supl,xsurf1 !define varible needed sumap,sx,s,x$!supl,xsurf1,sx sumap,sy,s,y$!supl,xsurf1,sy sumap,sz,s,z$!supl,xsurf1,sz sumap,sxy,s,xy$!supl,xsurf1,sxy sumap,syz,s,yz$!supl,xsurf1,syz sumap,sxz,s,xz$!supl,xsurf1,sxz !--------The following code varies in the normal direction of section definded------ !------------------------This is for the direction of x aixs------------------ ! get section force: Ax SFy SFz sueval,Ax,sx,intg sueval,SFy,sxy,intg sueval,SFz,sxz,intg ! get section center sueval,A,DA,sum! GCX, GCY, GCZ---global Cartesian coordinates at each point on the surface. sueval,Sx,GCx,intg$dx=Sx/A sueval,Sy,GCy,intg$dy=Sy/A sueval,Sz,GCz,intg$dz=Sz/A ! get the section moment: Mz,My,TOR
ansys中如何生成命令流方法： GUI是：Utility Menu>File>Write DB Log File 怎么用ansys绘制弯矩，剪力图：GUI: General Postproc-> lot Result->Contour Plot- >Line Element Result 弹出画单元结果的对话框，分别在Labi和Labj依次选取SMIS6和SMIS12（弯矩图）、SMIS1和SMIS7（轴力图）、SMIS2和SMIS8（剪力图） ! 建立单元表 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,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,-0.8,0 /image,save,'Bending_moment_%T%',jpg ANSYS中弯矩、剪力图的绘制 GUI: General Postproc-plot Result-Contour Plot-Line Element Result 弹出画单元结果的对话框，分别在Labi和Labj依次选取SMIS6和SMIS12（弯矩图）、SMIS1和SMIS7（轴力图）、SMIS2和SMIS8（剪力图）
查轴力：首先定义单元表 grneral postproc>element table >define table add 左侧选by sequence num,右侧选择smisc, 在下面输入smisc,1 然后在plot results>contour plot》line elem res 查看 弯矩 1.绘制弯矩图 建立弯矩单元表。例如梁单元 i节点单元表名称为imom，j节点单元表名称为jmom， 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点弯矩 plls,MI,MJ 2.标注弯矩图 PLOTCTRLS>>NUMBERING>>SVAL ON即可在画出弯矩图的同时在图上标出弯矩值的大小 3.调整弯矩图 如果弯矩图方向错误，则绘制弯矩图命令为 plls,imom,jmom，-1 同一个节点处两边的单元力有细微差别， 导致力数字标注出现重影。观察上面整体轴力图也可以发现， 一段一段的，好像马赛克，其实上面整体弯矩图也是，不过不是 很明显罢了。这是EULER-BEONOULI梁理论以及ANSYS输出定义造成 的（详细原因就不展开了，看看梁理论的书和ANSYS的说明吧）。 为了修正重影和节点两边力值不一样的问题，遍制了宏文件ITFAVG.MAC 命令文件容如下： !--------------------------------------------------------------------- !宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO） !获取线性单元力，并对单元边界处的力进行平衡 !输入信息 !力类型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ *ASK,ITFTYPE,'PLEASE INPUT THE TYPE OF INTERNAL FORCE','MMOMY' !需处理的单元包 *ASK,EASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF ELEMENTS TO BE PROCESSED!', 'EOUTER' !需处理的节点包
ABAQUS中如何通过cutting surface和section print输出桩的轴力 经过两个星期的摸索与学习，今天终于学会了桩轴力的输出。现总结如下： 1.主要步骤是先定义截面cutting surface，然后用section print输出轴力sof。 2.所有操作均是在inp文件中进行修改的，而不是ABAQUS/CAE中的编辑关键词（edit keywords）。 原因：在CAE中编辑关键词是可以修改inp文件，但CAE并不能识别所有的inp文件关键词，下面将举例说明。 3.最后提交的inp文件也不是在CAE中导入模型文件（import model），然后提交job进行运算的，而是在ABAQUS命令窗口（小黑屏）进行的。 原因同2中的一样，CAE并不能识别关键词*section print。 好了，下面开始详细的步骤讲解吧！ 第一步：定义截面（cutting surface），具体的关键语句为： *surface,type=cutting surface,name=cutsurf-1 0.6,25,0,0,1,0 Set-pile 解读： 第一行，定义surface、surface类型以及名称。 第二行，定义截面上的一点（0.6,25,0）以及截面的法向量（0,1,0）。法向量不一定是单位向量。 第三行，截面所在的单元或集合。这个集合可以是事先在CAE里定义好的。 需要注意的是，这个cutting surface是垂直于桩径方向的一个桩截面，而不是桩的侧表面。我一开始理解错了！ 此关键句在inp文件中的位置是在*Assembly, name=Assembly这一行之后，我试了下在放在材料定义之后，运算不成功。估计是因为我的Set-pile是在assembly里定义的，而不是在part里定义的单元集。若是在part 里定义了桩的集合，是不是可以将此关键句移动到材料定义之后，这个我倒没试。 如果是在CAE中通过编辑关键词来添加上述语句，将会有下面的错误提示：
ansys中如何生成命令流方法： 令狐采学 GUI是：Utility Menu>File>Write DB Log File 怎么用ansys绘制弯矩，剪力图：GUI:General Postproc- > lot Result->Contour Plot->Line Element Result弹出画单元结果的对话框，分别在Labi和Labj依次选取SMIS6和SMIS12（弯矩图）、SMIS1和SMIS7（轴力图）、SMIS2和SMIS8（剪力图） ! 建立单元表 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,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,-0.8,0 /image,save,'Bending_moment_%T%',jpg ANSYS中弯矩、剪力图的绘制 GUI: General Postproc-plot Result-Contour Plot-Line Element Result弹出画单元结果的对话框，分别在Labi和Labj依次选取SMIS6和SMIS12（弯矩图）、SMIS1和SMIS7（轴力图）、SMIS2和SMIS8（剪力图） ! 建立单元表 ETABLE,NI,SMISC,1 !单元I点轴力 ETABLE,NJ,SMISC,7 !单元J点轴力 ETABLE,QI,SMISC,2 !单元I点剪力 ETABLE,QJ,SMISC,8 !单元J点剪力
第三章静定梁与静定刚架 目的要求：熟练掌握静定梁和静定刚架的内力计算和内力图的绘制方法，熟练掌握绘制弯矩图的叠加法及内力图的形状特征，掌握绘制弯矩图的技巧。掌握多跨静定梁的几何组成特点和受力特点。能恰当选取隔离体和平衡方程计算静定结构的内力。 重点：截面法、微分关系的应用、简支梁叠加法。 难点：简支梁叠加法，绘制弯矩图的技巧 §3-1 单跨静定梁 1．反力 常见的单跨静定梁有简支梁、伸臂梁和悬臂梁三种，如图3-1(a)、(b)、(c)所示，其支座反力都只有三个，可取全梁为隔离体，由三个平衡条件求出。 图3-１ 2．内力 截面法是将结构沿所求内力的截面截开，取截面任一侧的部分为隔离体，由平衡条件计算截面内力的一种基本方法。 （1）内力正负号规定 轴力以拉力为正；剪力以绕隔离体有顺时 针转动趋势者为正；弯矩以使梁的下侧纤维受 拉者为正，如图3-2(b)所示。 （2）梁的内力与截面一侧外力的关系图3-2 1) 轴力的数值等于截面一侧的所有外力(包括荷载和反力)沿截面法线方向的投影代数和。 2) 剪力的数值等于截面一侧所有外力沿截面方向的投影代数和。 3) 弯矩的数值等于截面一侧所有外力对截面形心的力矩代数和。 3．利用微分关系作内力图 表示结构上各截面内力数值的图形称为内力图。内力图常用平行于杆轴线的坐标表示截面位置(此坐标轴常称为基线)，而用垂直于杆轴线的坐标(亦称竖标)表示内力的数值而绘出的。弯矩图要画在杆件的受拉侧，不标注正负号；剪力图和轴力图将正值的竖标绘在基线的上方，同时要标注正负号。绘内力图的基本方法是先写出内力方程，即以变量x表示任意截面的位置并由截面法写出所求内力与x之间的函数关系式，然后由方程作图。但通常采用的是利用微分关系来作内力图的方法。 （1）荷载与内力之间的微分关系