ANSYS中弯矩 剪力图的绘制(梁)

合集下载

ANSYS绘制弯矩

ANSYS绘制弯矩

ANSYS绘制弯矩、剪力、轴力图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,MJ2.标注弯矩图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!','NOUTER'!无需处理的节点包*ASK,UNASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NODE!(NONE IF THERE'S NO SUCH COMPONENT)','NONE'/POST1!输入信息:内力类型,欲处理单元的集合,欲处理节点的集合!ITFTYPE='MMOMY'!EASSEMBLY='EOUTER'!NASSEMBLY='NOUTER'!按内力类型确定ANSYS输出信息SMISC的编号*IF,ITFTYPE,EQ,'MFORX',THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOMX',THENITFINUM=4ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6ITFJNUM=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,1UNNOD(I+1)=NDNEXT(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!选出所需的单元和节点包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!获得当前选中单元总数(存入变量SELELENUM)*GET,SELELENUM,ELEM,0,COUNT!定义长度为SELELENUM的数组(ELENUM),以存放选中单元的单元编号*DIM,ELENUM,ARRAY,SELELENUM!将选中单元的编号按顺序存入数组ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNEXT(I)*ENDDO!获得当前选中节点总数(存入变量SELNODNUM)*GET,SELNODNUM,NODE,0,COUNT!定义长度为SELNODNUM的数组(NODNUM),以存放选中单元的单元编号*DIM,NODNUM,ARRAY,SELNODNUM!将选中单元的编号按顺序存入数组NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNEXT(I)*ENDDO!定义所需的线性单元内力ETABLE,节点I的内力存入数组ITNFI,!节点J的内力存入数组ITNFJETABLE,ITNFI,SMISC,ITFINUMETABLE,ITNFJ,SMISC,ITFJNUM!定义所需的结果数组,并将其置零ETABLE,ITNFINEO,SMISC,5SADD,ITNFINEO,ITNFI,,1ETABLE,ITNFJNEO,SMISC,11SADD,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),THENINDEX=1*ELSE*ENDIF*ENDDO*ELSE*ENDIF*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,I节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,J节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THENELEJ=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,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!END OF ITFAVG.MAC(2)对体和面来说,ANSYS默认的结果输出格式是云图格式,而这种彩色云图打印为黑白图像时对比很不明显,无法表达清楚,对于发表文章非常不便。

ansys教程——梁解析

ansys教程——梁解析
第五章

5. 梁
• 梁单元 是线单元,用来创建3-D结构的一维理想体。 • 梁单元比实体和壳单元更有效,经常用于工业领域中:
– 建筑结构 – 桥梁和道路 – 公共交通 (有轨电车, 火车, 公共汽车) – 等
Training Manual
January 30, 2001 Inventory #001443 5-2

...梁网格划分
步骤1: 线属性 • 梁网格划分的线属性包括:
– 材料号 – 横截面号 – 定位关键点
Training Manual
INTRODUCTION TO ANSYS 5.7 - Part 2
• 相对于梁轴线,横截面是怎样定位。
• 必须指定所有横截面类型. • 单个关键点可以分配给多条线 ( 即,不需要为每条线指定单个关键 点 )。
横截面 • 对 BEAM188 和 189单元的完整定义包括对横 截面属性的定义。 • BeamTool提供了方便的操作.
– Preprocessor > Sections > Common Sectns...
Training Manual
INTRODUCTION TO ANSYS 5.7 - Part 2

...梁网格划分
步骤2 :线分隔 • 对 BEAM188 和 189 单元,不主张把整个梁当作一个单元 • 使用 Mesh Tool的 “Size Controls”指定想要的线分隔数 (或用 LESIZE 命令).
Training Manual
January 30, 2001 Inventory #001443 5-13

...梁网格划分
Training Manual

Ansys绘制梁的弯矩图和剪力图(图文教程+APDL命令流)

Ansys绘制梁的弯矩图和剪力图(图文教程+APDL命令流)

Y ZX
-.625E-09 3482.91
6965.81
10448.7
13931.6
17414.5
20897.4
24380.3
27863.2
31346.2
Step17 弯矩图
9
1 LINE STRESS STEP=1 SUB =1 TIME=1 SZ_I SZ_J MIN =-15673.1 ELEM=18 MAX =15076.9 ELEM=61
6
Step11 施加集中载荷 Step12 施加梯度载荷
7
Step13 力学模型
Step14 计算
Step15 计算完成
8
1 NODAL SOLUTION
STEP=1
SUB =1
TIME=1
UY
(AVG)
RSYS=0
DMX =.073299
SMN =-.073299
Y MZX X
MN
-.073299
*do,i,31,45,1 sfbeam,i,2,pres,1e3/15*(i-31),1e3/15*(i-31+1) *enddo allsel,all !!!!计算 /solu
11
solve finish /post1 ETABLE,Mz_I,SMISC, 3 ETABLE,Mz_J,SMISC, 16 PLLS, Mz_I, Mz_J ETABLE,Sz_I,SMISC, 6 ETABLE,Sz_J,SMISC, 19 PLLS, Sz_I, Sz_J
27863.2
31346.2
图 1 结果图
1
图文教程
Step 1 选择单元类型 Step 2 设定材料属性
2

ansys查轴力弯矩新版

ansys查轴力弯矩新版

查轴力:首先定义单元表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,MJ2.标注弯矩图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',THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOMX',THENITFINUM=4ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6ITFJNUM=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,1UNNOD(I+1)=NDNEXT(I)*ENDDOUNNODNUM=0*ENDIF!选出所需的单元和节点包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!获得当前选中单元总数(存入变量SELELENUM)*GET,SELELENUM,ELEM,0,COUNT!定义长度为SELELENUM的数组(ELENUM),以存放选中单元的单元编号*DIM,ELENUM,ARRAY,SELELENUM!将选中单元的编号按顺序存入数组ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNEXT(I)*ENDDO!获得当前选中节点总数(存入变量SELNODNUM)*GET,SELNODNUM,NODE,0,COUNT!定义长度为SELNODNUM的数组(NODNUM),以存放选中单元的单元编号*DIM,NODNUM,ARRAY,SELNODNUM!将选中单元的编号按顺序存入数组NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNEXT(I)*ENDDO!定义所需的线性单元内力ETABLE,节点I的内力存入数组ITNFI,!节点J的内力存入数组ITNFJETABLE,ITNFI,SMISC,ITFINUMETABLE,ITNFJ,SMISC,ITFJNUM!定义所需的结果数组,并将其置零ETABLE,ITNFINEO,SMISC,5SADD,ITNFINEO,ITNFI,,1ETABLE,ITNFJNEO,SMISC,11SADD,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),THENINDEX=1*ELSE*ENDIF*ENDDO*ELSE*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,I节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,J节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THENELEJ=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,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!END OF ITFAVG.MAC(2)对体和面来说,ANSYS默认的结果输出格式是云图格式,而这种彩色云图打印为黑白图像时对比很不明显,无法表达清楚,对于发表文章非常不便。

ansys APDL操作计算弯矩

ansys APDL操作计算弯矩

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,MJ2.标注弯矩图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!','NOUTER'!无需处理的节点包*ASK,UNASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NODE!(NONE IF THERE'S NO SUCH COMPONENT)','NONE'/POST1!输入信息:内力类型,欲处理单元的集合,欲处理节点的集合!ITFTYPE='MMOMY'!EASSEMBLY='EOUTER'!NASSEMBLY='NOUTER'!按内力类型确定ANSYS输出信息SMISC的编号*IF,ITFTYPE,EQ,'MFORX',THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOMX',THENITFINUM=4ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6ITFJNUM=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,1UNNOD(I+1)=NDNEXT(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!选出所需的单元和节点包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!获得当前选中单元总数(存入变量SELELENUM)*GET,SELELENUM,ELEM,0,COUNT!定义长度为SELELENUM的数组(ELENUM),以存放选中单元的单元编号*DIM,ELENUM,ARRAY,SELELENUM!将选中单元的编号按顺序存入数组ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNEXT(I)*ENDDO!获得当前选中节点总数(存入变量SELNODNUM)*GET,SELNODNUM,NODE,0,COUNT!定义长度为SELNODNUM的数组(NODNUM),以存放选中单元的单元编号*DIM,NODNUM,ARRAY,SELNODNUM!将选中单元的编号按顺序存入数组NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNEXT(I)*ENDDO!定义所需的线性单元内力ETABLE,节点I的内力存入数组ITNFI,!节点J的内力存入数组ITNFJETABLE,ITNFI,SMISC,ITFINUMETABLE,ITNFJ,SMISC,ITFJNUM!定义所需的结果数组,并将其置零ETABLE,ITNFINEO,SMISC,5SADD,ITNFINEO,ITNFI,,1ETABLE,ITNFJNEO,SMISC,11SADD,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),THENINDEX=1*ELSE*ENDIF*ENDDO*ELSE*ENDIF*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,I节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,J节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THENELEJ=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,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!END OF ITFAVG.MAC(2)对体和面来说,ANSYS默认的结果输出格式是云图格式,而这种彩色云图打印为黑白图像时对比很不明显,无法表达清楚,对于发表文章非常不便。

快速绘制梁的剪力图和弯矩图

快速绘制梁的剪力图和弯矩图

快速绘制梁的剪力图和弯矩图哎呀,今天咱们来聊聊一个很有意思的话题:快速绘制梁的剪力图和弯矩图。

这可不是啥高深莫测的技术活儿,就是让大家了解一下梁的结构特点,然后用简单的方法画出它的受力情况。

别看这个过程好像挺复杂的,其实只要掌握了一定的技巧,就能轻松搞定!咱们得了解一下梁的基本结构。

梁是建筑物中常见的一种构件,它的主要作用是承受垂直于其长度方向的荷载,将这些荷载传递到支撑结构上。

梁通常由两个平行的钢板或混凝土板组成,这两个板之间的距离叫做梁的宽度。

梁的长度可以是任意值,但通常会根据建筑物的实际需求进行设计。

接下来,我们来说说如何绘制梁的剪力图。

剪力图是一种表示梁在不同位置受到剪力作用时的力的分布情况的图形。

要画出剪力图,我们需要知道梁的几何尺寸、材料的抗剪强度以及荷载的大小。

有了这些信息,我们就可以用一些简单的公式来计算梁在不同位置受到的剪力大小了。

我们要把梁划分成若干个小的单元,这样才能更好地观察和分析梁的受力情况。

一般来说,我们可以把梁划分成若干个矩形或梯形的小单元。

然后,我们可以用一些公式来计算每个小单元受到的剪力大小。

这些公式包括:1. 矩形单元的剪力大小 = 矩形面积× 材料抗剪强度× 剪切模量2. 梯形单元的剪力大小 = (梯形面积× 材料抗剪强度× 剪切模量) ÷ 2有了这些公式,我们就可以开始画剪力图了。

我们在纸上画出一个坐标系,然后用坐标轴表示梁的长度和位置。

接着,我们可以用不同的颜色或线条来表示不同位置的剪力大小。

例如,我们可以用红色表示剪力较大的区域,用蓝色表示剪力较小的区域。

我们可以在图上加上一些注释,说明各个区域的具体含义。

现在,我们再来说说如何绘制梁的弯矩图。

弯矩图是一种表示梁在不同位置受到弯曲作用时的力的分布情况的图形。

要画出弯矩图,我们同样需要知道梁的几何尺寸、材料的抗弯强度以及荷载的大小。

有了这些信息,我们就可以用一些简单的公式来计算梁在不同位置受到的弯矩大小了。

ANSYS中弯矩、剪力图的绘制

ANSYS中弯矩、剪力图的绘制

.ansys中如何生成命令流方法:GUI是:Utility Menu>File>Write DB Log File怎么用ansys绘制弯矩,剪力图:GUI: General Postproc-&gt 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 diagramPLLS,NI,NJ,1.0,0/image,save,'Axial_force_%T%',jpg! 画剪力分布图/TITLE,Shearing force diagramPLLS,QI,QJ,1.0,0/image,save,'Shearing_force_%T%',jpg! 画弯矩分布图/TITLE,Bending moment diagramPLLS,MI,MJ,-0.8,0/image,save,'Bending_moment_%T%',jpgANSYS中弯矩、剪力图的绘制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点剪力ETABLE,MI,SMISC,6 !单元I点弯矩ETABLE,MJ,SMISC,12 !单元J点弯矩! 更新单元表ETABLE,REFL! 画轴力分布图/TITLE,Axial force diagramPLLS,NI,NJ,1.0,0/image,save,'Axial_force_%T%',jpg! 画剪力分布图/TITLE,Shearing force diagramPLLS,QI,QJ,1.0,0/image,save,'Shearing_force_%T%',jpg! 画弯矩分布图/TITLE,Bending moment diagramPLLS,MI,MJ,-0.8,0/image,save,'Bending_moment_%T%',jpg另:自定义截面梁剪力弯矩显示finish/clear/verify/replot!自定义截面/prep7et,1,plane82rectng,0,1.0,0,0.6,cyl4,0.28,0.25,0.18,-180,cyl4,0.28,0.35,0.18,180,cyl4,0.72,0.25,0.18,-180,cyl4,0.72,0.35,0.18,180,rectng,0.1,0.46,0.25,0.35,rectng,0.54,0.9,0.25,0.35,asel,u,,,1.cm,area0,areaallsel,allasba,1,area0esize,0.1amesh,all!.读入截面文件secwrite,jiemian,sect,,1aclear,alladele,all,,1ldele,all,,,1finish/clear/prep7et,1,beam44keyopt,1,6,1mp,dens,1,2600mp,ex,1,3.06e10mp,prxy,1,0.2sectype,1,beam,mesh,sect1secoffset,cent,,,secread,'jiemian','sect','',meshk,1k,2,10k,3,0,3lstr,1,2latt,1,,1,,3,,1lesize,all,0.5lmesh,all/eshape,1eplotdk,1,ux,0,,,uy,uzdk,2,uy,0,,,uzf,12,fy,-1/soluantype,staticsolvefinish/post1pldisp,2plnsol,u,y,2!显示剪力etable,sheari,smisc,3etable,shearj,smisc,9plls,sheari,shearj,-1!显示弯矩etable,mforcei,smisc,5etable,mforcej,smisc,11plls,mforcei,mforcej,-1ansys如何绘制弯矩图Ansy中弯矩图,云图绘制总结在回答别人问题时,利用前人的回复和总结,自己进行了总结改正,发表在这里,供各位参考(1)ANSYS弯矩等可以直接标注在图上吗?如何实现?如果三维问题,在剖面上标出某一结构的轴心力、弯矩等,如何实现(2)后处理图形,其等值线的数值能否直接标注在图上,而不是采用图例的形式后处理结果往往用云图表示,下跟一图例表示数值大小,能够实现等值线直接标注在图上回答(1)1.绘制弯矩图建立弯矩单元表。

ansys查轴力弯矩新版[精品文档]

ansys查轴力弯矩新版[精品文档]

查轴力:首先定义单元表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,MJ2.标注弯矩图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',THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOMX',THENITFINUM=4ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6ITFJNUM=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,1UNNOD(I+1)=NDNEXT(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!选出所需的单元和节点包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!获得当前选中单元总数(存入变量SELELENUM)*GET,SELELENUM,ELEM,0,COUNT!定义长度为SELELENUM的数组(ELENUM),以存放选中单元的单元编号*DIM,ELENUM,ARRAY,SELELENUM!将选中单元的编号按顺序存入数组ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNEXT(I)*ENDDO!获得当前选中节点总数(存入变量SELNODNUM)*GET,SELNODNUM,NODE,0,COUNT!定义长度为SELNODNUM的数组(NODNUM),以存放选中单元的单元编号*DIM,NODNUM,ARRAY,SELNODNUM!将选中单元的编号按顺序存入数组NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNEXT(I)*ENDDO!定义所需的线性单元内力ETABLE,节点I的内力存入数组ITNFI,!节点J的内力存入数组ITNFJETABLE,ITNFI,SMISC,ITFINUMETABLE,ITNFJ,SMISC,ITFJNUM!定义所需的结果数组,并将其置零ETABLE,ITNFINEO,SMISC,5SADD,ITNFINEO,ITNFI,,1ETABLE,ITNFJNEO,SMISC,11SADD,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),THENINDEX=1*ELSE*ENDIF*ENDDO*ELSE*ENDIF*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,I节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,J节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THENELEJ=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,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!END OF ITFAVG.MAC(2)对体和面来说,ANSYS默认的结果输出格式是云图格式,而这种彩色云图打印为黑白图像时对比很不明显,无法表达清楚,对于发表文章非常不便。

ansys查轴力弯矩新版

ansys查轴力弯矩新版

ansys查轴力弯矩新版查轴力:首先定义单元表grneralpostproc>elementtable>definetableadd左侧选bysequencenum,右侧选择smisc,在下面输入smisc,1然后在plotresults>contourplot》lineelemres查看弯矩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,mj2.标注弯矩图Plotctrls>>编号>>svalon可以绘制弯矩图,同时在图上标记弯矩值。

3调整弯矩图如果弯矩图方向错误,则绘制弯矩图命令为plls,imom,jmom,-1同一节点两侧构件的内力略有差异,导致内力数字标注出现重影。

观察上面整体轴力图也可以发现,一段一段的,好像马赛克,其实上面整体弯矩图也是,不过不是很明显罢了。

这是euler-beonouli梁理论以及ansys输出定义造成的(详细原因就不展开了,看看梁理论的书和ansys的说明吧)。

为了纠正节点两侧的重影和不同内力值问题,在MAC命令文件的所有内容中都创建了宏文件itfavg,如下所示:!---------------------------------------------------------------------!宏:itfavg.mac(internalforceaveragemacro)!获取线性单元内力,并对单元边界处的内力进行平衡!输入信息! 内力类型:mforx、mfory、mforz、mmomx、mmomy、mmomz*ask,itftype,'pleaseinputthetypeofinternalforce','mmomy'!需处理的单元包! 要处理的节点包!无需处理的节点包! 输入信息:内力类型、要处理的单元集、要处理的节点集!itftype='mmomy'!易组装!nassembly='nouter'!按内力类型确定ansys输出信息smisc的编号*if,itftype,eq,'mforx',thenitfinum=1itfjnum=7*elseif,itftype,eq,'mfory',然后itfinum=2itfjnum=8*elseif,itftype,eq,'mforz',thenitfinum=3itfjnum=9*elseif,ITF类型,eq,'mmomx',然后ITFINUM=4itfjnum=10*elseif,itftype,eq,'mmomy',thenitfinum=5itfjnum=11*elseif,itftype,eq,'mmomz',然后itfinum=6itfjnum=12*else*endif!对不需平均的节点进行处理*if,unassembly,ne,'none',then! 选择将不被处理的节点包,并获取将不被处理的节点数cmsel,s,unassembly*get,unnodnum,node,0,count!定义长度为unnodnum的数组(unnod),以存放选中单元的单元编号*dim,unnod,array,unnodnum! 将所选单元格的编号存储在数组unod*do,I,0,unodnum-1,1 unod(I+1)=ndnext(I)*enddo中*elseunnodnum=0*endif!选出所需的单元和节点包cmsel,s,eassemblycmsel,s,nassembly! 获取当前选定单元格的总数(存储在变量selenium中)*Get、selenium、elem、0、count!定义长度为selelenum的数组(elenum),以存放选中单元的单元编号*dim,elenum,array,selelenum!将选中单元的编号按顺序存入数组elenum*do,i,0,selelenum-1,1elenum(i+1)=elnext(i)*enddo! 获取当前选定节点的总数(存储在变量selnodnum中)*Get,selnodnum,node,0,count!定义长度为selnodnum的数组(nodnum),以存放选中单元的单元编号*dim,nodnum,array,selnodnum! 将所选单元格的编号存储在数组nodnum*do,I,0,selnodnum-1,1nodnum(I+1)=ndnext(I)*enddo中!定义所需的线性单元内力etable,节点i的内力存入数组itnfi,!节点j的内力存入数组itnfjetable,itnfi,smisc,itfinumetable,itnfj,smisc,itfjnum! 定义所需的结果数组,并将其设置为0 etable、itnfineo、smisc、5sadd、itnfineo、itnfi、、1etable、itnfjneo、smisc、11sadd、itnfjneo、itnfj、、1*do、K、1、selnodnum、1!正在处理不需要平均索引=0的节点*if,unnodnum,ge,1,then*do,j,1,unnodnum*如果,nodnum(k),eq,unnod(j),那么index=1*else*endif*enddo*else*endif*do,j,1,selelenum,1!选出和节点k相连的线性单元中,i节点(对线性单元而言)为节点k的单元编号*if,nelem(elenum(j),1),eq,nodnum(k),thenelei=elenum(j)*exit*else*endif*enddo*do,j,1,selelenum,1!选出和节点k相连的线性单元中,j节点(对线性单元而言)为节点k的单元编号*if,nelem(elenum(j),2),eq,nodnum(k),thenelej=elenum(j)*exit*else*endif*enddo*如果,索引,等式,0,那么*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,etavedetab,elej,itnfjneo,etave*else*endif*else*endif*enddo/udoc,1,标志,关闭plls,itnfineo,itnfjneo!endofitfavg.mac(2)对于体积和表面,ANSYS的默认结果输出格式为云图格式,当该彩色云图打印为黑白图像时,对比度不明显,无法清晰表达,这对发表文章非常不方便。

ansys提出实体单元建立弯矩

ansys提出实体单元建立弯矩

对于3D实体单元的内力计算,可采用面操作(应该是V8以上吧,没有细考!)。

例如用悬臂梁采用3D实体单元,其某个截面的的内力计算如下,并具有详细解释。

这种方法较原来的更方便,大家不放一用。

!******************************************************!EX4.31 面操作及悬臂梁的内力计算FINISHCLEAR/PREP7ET,1,SOLID95MP,EX,1,2E11MP,PRXY,1,0.3 !定义单元类型、材料特性BLC4,2,3,0.2,0.3,4DA,2,ALLFK,1,FY,-2E4FK,3,FY,-2E4 !创建几何模型、加约束和荷载FK,3,FX,0.8E4FK,4,FX,0.8E4SFA,1,1,PRES,1E6 !施加荷载ESIZE,0.05VMESH,ALLFINISH/SOLUSOLVE !生成有限元模型并求解FINISH/POST1 !进入后处理层WPOFF,,,2SUCR,SUZ2,CPLANE,3 !移动工作平面、创建面SUZ2SUMAP,MYSX,S,XSUMAP,MYSY,S,Y !映射X和Y方向应力SUMAP,MYSZ,S,ZSUMAP,MYSXY,S,XY !映射Z和XY方向应力SUMAP,MYSYZ,S,YZSUMAP,MYSXZ,S,XZ !映射YZ和XZ方向应力SUPL,SUZ2SUPL,SUZ2,MYSZ !显示面本身、面项MYSZSUPL,SUZ2,MYSYZSUPL,SUZ2,MYS !显示面项MYSYZ,矢量显示应力SUPR,ALL,MYSZ !列表显示MYSZ面项SUEVAL,XFORCE,MYSXZ,INTG !求截面上FX,理论结果为-16000,误差1%SUEVAL,YFORCE,MYSYZ,INTG !求截面上FY,理论结果为40000,误差0.5%SUEVAL,ZFORCE,MYSZ,INTG !求截面上FZ,理论结果为-6000,没有误差SUEVAL,MYA,DA,SUM !求截面面积并赋给MYA变量SUEVAL,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在面上积分得到MX1SUCALC,SZGCX,MYSZ,MULT,GCX !计算MYSZ×GCZ,并赋给面项SZGCXSUEVAL,MY1,SZGCX,INTG !对面项SZGCX在面上积分得到MY1!上述弯矩基于总体直角坐标系原点而言的,应对面积重心取矩,将内力简化到面积重心上MX1=MX1-ZFORCE*MYYA !理论结果为80000,误差为0.08%MY1=MY1-ZFORCE*MYXA !理论结果-32000,误差为0.2%。

ik4ANSYS中弯矩、剪力图的绘制2

ik4ANSYS中弯矩、剪力图的绘制2

ansys中如何生成命令流方法:GUI是:Utility Menu>File>Write DB Log File怎么用ansys绘制弯矩,剪力图:GUI: General Postproc-&gt 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 diagramPLLS,NI,NJ,1.0,0/image,save,'Axial_force_%T%',jpg! 画剪力分布图/TITLE,Shearing force diagramPLLS,QI,QJ,1.0,0/image,save,'Shearing_force_%T%',jpg! 画弯矩分布图/TITLE,Bending moment diagramPLLS,MI,MJ,-0.8,0/image,save,'Bending_moment_%T%',jpgANSYS中弯矩、剪力图的绘制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点剪力ETABLE,MI,SMISC,6 !单元I点弯矩ETABLE,MJ,SMISC,12 !单元J点弯矩 ! 更新单元表ETABLE,REFL! 画轴力分布图/TITLE,Axial force diagramPLLS,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另:自定义截面梁剪力弯矩显示finish/clear/verify/replot!自定义截面/prep7et,1,plane82rectng,0,1.0,0,0.6,cyl4,0.28,0.25,0.18,-180,cyl4,0.28,0.35,0.18,180,cyl4,0.72,0.25,0.18,-180,cyl4,0.72,0.35,0.18,180,rectng,0.1,0.46,0.25,0.35,rectng,0.54,0.9,0.25,0.35,asel,u,,,1cm,area0,areaallsel,allasba,1,area0esize,0.1amesh,all!读入截面文件secwrite,jiemian,sect,,1 aclear,alladele,all,,1ldele,all,,,1finish/clear/prep7et,1,beam44keyopt,1,6,1mp,dens,1,2600mp,ex,1,3.06e10 mp,prxy,1,0.2 sectype,1,beam,mesh,sect1 secoffset,cent,,,secread,'jiemian','sect','',mesh k,1k,2,10k,3,0,3lstr,1,2latt,1,,1,,3,,1lesize,all,0.5lmesh,all/eshape,1eplotdk,1,ux,0,,,uy,uzdk,2,uy,0,,,uzf,12,fy,-1/soluantype,staticsolvefinish/post1pldisp,2plnsol,u,y,2!显示剪力etable,sheari,smisc,3etable,shearj,smisc,9plls,sheari,shearj,-1!显示弯矩etable,mforcei,smisc,5etable,mforcej,smisc,11plls,mforcei,mforcej,-1ansys如何绘制弯矩图Ansy中弯矩图,云图绘制总结在回答别人问题时,利用前人的回复和总结,自己进行了总结改正,发表在这里,供各位参考(1)ANSYS弯矩等可以直接标注在图上吗?如何实现?如果三维问题,在剖面上标出某一结构的轴心力、弯矩等,如何实现(2)后处理图形,其等值线的数值能否直接标注在图上,而不是采用图例的形式后处理结果往往用云图表示,下跟一图例表示数值大小,能够实现等值线直接标注在图上回答(1)1.绘制弯矩图建立弯矩单元表。

ANSYS中弯矩 剪力图的绘制(梁)

ANSYS中弯矩 剪力图的绘制(梁)

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点剪力ETABLE,MI,SMISC,6 !单元I点弯矩ETABLE,MJ,SMISC,12 !单元J点弯矩! 更新单元表ETABLE,REFL! 画轴力分布图/TITLE,Ax ial force diagramPLLS,NI,NJ,1.0,0/image,save,'Axial_force_%T%',jpg! 画剪力分布图/TITLE,Shearing force diagramPLLS,QI,QJ,1.0,0/image,save,'Shearing_force_%T%',jpg! 画弯矩分布图/TITLE,Bending moment diagramPLLS,MI,MJ,-0.8,0/image,save,'Bending_moment_%T%',jpg另:自定义截面梁剪力弯矩显示finish/clear/verify/replot!自定义截面/prep7et,1,plane82rectng,0,1.0,0,0.6,cyl4,0.28,0.25,0.18,-180,cyl4,0.28,0.35,0.18,180,cyl4,0.72,0.25,0.18,-180,cyl4,0.72,0.35,0.18,180,rectng,0.1,0.46,0.25,0.35,rectng,0.54,0.9,0.25,0.35,asel,u,,,1cm,area0,areaallsel,allasba,1,area0esize,0.1amesh,all!读入截面文件secwrite,jiemian,sect,,1 aclear,alladele,all,,1ldele,all,,,1finish/clear/prep7et,1,beam44keyopt,1,6,1mp,dens,1,2600mp,ex,1,3.06e10 mp,prxy,1,0.2 sectype,1,beam,mesh,sect1 secoffset,cent,,,secread,'jiemian','sect','',mesh k,1k,2,10k,3,0,3lstr,1,2latt,1,,1,,3,,1lesize,all,0.5lmesh,all/eshape,1eplotdk,1,ux,0,,,uy,uzdk,2,uy,0,,,uzf,12,fy,-1/soluantype,staticsolvefinish/post1pldisp,2plnsol,u,y,2!显示剪力etable,sheari,smisc,3 etable,shearj,smisc,9plls,sheari,shearj,-1!显示弯矩etable,mforcei,smisc,5 etable,mforcej,smisc,11 plls,mforcei,mforcej,-1。

ansys查轴力弯矩

ansys查轴力弯矩
*DO,I,0,UNNODNUM-1,1
UNNOD(I+1)=NDNEXT(I)
*ENDDO
*ELSE
UNNODNUM=0
*ENDIF
!选出所需的单元和节点包
CMSEL,S,EASSEMBLY
CMSEL,S,NASSEMBLY
!获得当前选中单元总数(存入变量SELELENUM)
*GET,SELELENUM,ELEM,0,COUNT
*GET,SELNODNUM,NODE,0,COUNT
!定义长度为SELNODNUM的数组(NODNUM),以存放选中单元的单元编号
*DIM,NODNUM,ARRAY,SELNODNUM
!将选中单元的编号按顺序存入数组NODNUM
*DO,I,0,SELNODNUM-1,1
NODNUM(I+1)=NDNEXT(I)
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,多试几次,直到疏密合适
/image,save,'Shearing_force_%T%',jpg
! 画弯矩分布图
/TITLE,Bending moment diagram
PLLS,MI,MJ,-0.8,0
/image,save,'Bending_moment_%T%',jp

ansys建立弯矩方法

ansys建立弯矩方法

对于3D实体单元的内力计算,可采用面操作(应该是V8以上吧,没有细考!)。

例如用悬臂梁采用3D实体单元,其某个截面的的内力计算如下,并具有详细解释。

这种方法较原来的更方便,大家不放一用。

!******************************************************!EX4.31 面操作及悬臂梁的内力计算FINISHCLEAR/PREP7ET,1,SOLID95MP,EX,1,2E11MP,PRXY,1,0.3 !定义单元类型、材料特性BLC4,2,3,0.2,0.3,4!指定角点位置生成矩形或长方体DA,2,ALL!在面2上施加约束FK,1,FY,-2E4FK,3,FY,-2E4 !创建几何模型、加约束和荷载FK,3,FX,0.8E4FK,4,FX,0.8E4SFA,1,1,PRES,1E6 !在面上施加荷载 SF在节点上施加面载荷ESIZE,0.05VMESH,ALLFINISH/SOLUSOLVE !生成有限元模型并求解FINISH/POST1 !进入后处理层WPOFF,,,2SUCR,SUZ2,CPLANE,3 !移动工作平面、创建面SUZ2SUMAP,MYSX,S,XSUMAP,MYSY,S,Y !映射X和Y方向应力SUMAP,MYSZ,S,ZSUMAP,MYSXY,S,XY !映射Z和XY方向应力SUMAP,MYSYZ,S,YZSUMAP,MYSXZ,S,XZ !映射YZ和XZ方向应力SUPL,SUZ2SUPL,SUZ2,MYSZ !显示面本身、面项MYSZSUPL,SUZ2,MYSYZSUPL,SUZ2,MYS !显示面项MYSYZ,矢量显示应力SUPR,ALL,MYSZ !列表显示MYSZ面项SUEVAL,XFORCE,MYSXZ,INTG !求截面上FX,理论结果为-16000,误差1%SUEVAL,YFORCE,MYSYZ,INTG !求截面上FY,理论结果为40000,误差0.5%SUEVAL,ZFORCE,MYSZ,INTG !求截面上FZ,理论结果为-6000,没有误差SUEVAL,MYA,DA,SUM !求截面面积并赋给MYA变量SUEVAL,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在面上积分得到MX1SUCALC,SZGCX,MYSZ,MULT,GCX !计算MYSZ×GCZ,并赋给面项SZGCXSUEVAL,MY1,SZGCX,INTG !对面项SZGCX在面上积分得到MY1!上述弯矩基于总体直角坐标系原点而言的,应对面积重心取矩,将内力简化到面积重心上MX1=MX1-ZFORCE*MYYA !理论结果为80000,误差为0.08%MY1=MY1-ZFORCE*MYXA !理论结果-32000,误差为0.2%施加弯矩扭矩方法总结施加弯矩扭矩的方法其实不只三种,有很多种方法,在这里介绍其中的5种,并进行比较:1.将矩转换成一对的力偶,直接施加在对应的节点上面。

ansys查轴力弯矩

ansys查轴力弯矩
'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
! 更新单元表
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
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
6.修改
在显示等值线时,图形中太多的标识字母(A、B等)使得等值线图完全看不清楚。请问有没有什么方法可以去掉图形的这些标识字母呢
在plotCtrls下面的style下面的contours下面的contour Labeling,在弹出的contour labeling options里面的vector mode contour labels选off

ansys查轴力弯矩新版

ansys查轴力弯矩新版

查轴力:首先定义单元表 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,MJ2.标注弯矩图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',THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOMX',THENITFINUM=4ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6ITFJNUM=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,1UNNOD(I+1)=NDNEXT(I)*ENDDOUNNODNUM=0*ENDIF!选出所需的单元和节点包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!获得当前选中单元总数(存入变量SELELENUM)*GET,SELELENUM,ELEM,0,COUNT!定义长度为SELELENUM的数组(ELENUM),以存放选中单元的单元编号*DIM,ELENUM,ARRAY,SELELENUM!将选中单元的编号按顺序存入数组ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNEXT(I)*ENDDO!获得当前选中节点总数(存入变量SELNODNUM)*GET,SELNODNUM,NODE,0,COUNT!定义长度为SELNODNUM的数组(NODNUM),以存放选中单元的单元编号*DIM,NODNUM,ARRAY,SELNODNUM!将选中单元的编号按顺序存入数组NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNEXT(I)*ENDDO!定义所需的线性单元力ETABLE,节点I的力存入数组ITNFI,!节点J的力存入数组ITNFJETABLE,ITNFI,SMISC,ITFINUMETABLE,ITNFJ,SMISC,ITFJNUM!定义所需的结果数组,并将其置零ETABLE,ITNFINEO,SMISC,5SADD,ITNFINEO,ITNFI,,1ETABLE,ITNFJNEO,SMISC,11SADD,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),THENINDEX=1*ELSE*ENDIF*ENDDO*ELSE*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,I节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中,J节点(对线性单元而言)为节点K的单元编号*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THENELEJ=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,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!END OF ITFAVG.MAC(2)对体和面来说,ANSYS默认的结果输出格式是云图格式,而这种彩色云图打印为黑白图像时对比很不明显,无法表达清楚,对于发表文章非常不便。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

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点剪力
ETABLE,MI,SMISC,6 !单元I点弯矩
ETABLE,MJ,SMISC,12 !单元J点弯矩
! 更新单元表
ETABLE,REFL
! 画轴力分布图
/TITLE,Ax ial 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
另:自定义截面梁剪力弯矩显示
finish
/clear
/verify
/replot
!自定义截面
/prep7
et,1,plane82
rectng,0,1.0,0,0.6,
cyl4,0.28,0.25,0.18,-180,
cyl4,0.28,0.35,0.18,180,
cyl4,0.72,0.25,0.18,-180,
cyl4,0.72,0.35,0.18,180,
rectng,0.1,0.46,0.25,0.35,
rectng,0.54,0.9,0.25,0.35,
asel,u,,,1
cm,area0,area
allsel,all
asba,1,area0
esize,0.1
amesh,all
!读入截面文件
secwrite,jiemian,sect,,1 aclear,all
adele,all,,1
ldele,all,,,1
finish
/clear
/prep7
et,1,beam44
keyopt,1,6,1
mp,dens,1,2600
mp,ex,1,3.06e10 mp,prxy,1,0.2 sectype,1,beam,mesh,sect1 secoffset,cent,,,
secread,'jiemian','sect','',mesh k,1
k,2,10
k,3,0,3
lstr,1,2
latt,1,,1,,3,,1
lesize,all,0.5
lmesh,all
/eshape,1
eplot
dk,1,ux,0,,,uy,uz
dk,2,uy,0,,,uz
f,12,fy,-1
/solu
antype,static
solve
finish
/post1
pldisp,2
plnsol,u,y,2
!显示剪力
etable,sheari,smisc,3 etable,shearj,smisc,9
plls,sheari,shearj,-1
!显示弯矩
etable,mforcei,smisc,5 etable,mforcej,smisc,11 plls,mforcei,mforcej,-1。

相关文档
最新文档