ANSYS中弯矩剪力图的绘制.docx
《ANSYS中弯矩剪力图的绘制.docx》由会员分享,可在线阅读,更多相关《ANSYS中弯矩剪力图的绘制.docx(25页珍藏版)》请在冰豆网上搜索。
![ANSYS中弯矩剪力图的绘制.docx](https://file1.bdocx.com/fileroot1/2023-7/7/77e0d6fc-d35c-4b90-8dd3-807cc9f1f20a/77e0d6fc-d35c-4b90-8dd3-807cc9f1f20a1.gif)
ANSYS中弯矩剪力图的绘制
ansys中如何生成命令流方法:
GUI就是:
UtilityMenu〉DBLogFile
怎么用ansys绘制弯矩,剪力图:
GUI:
GeneralPostproc—>
lotResult->ContourPlot—>LineElementResultﻫ弹出画单元结果得对话框,分别在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 momentdiagram
PLLS,MI,MJ,—0、8,0
/image,save,’Bending_moment_%T%',jpg
ANSYS中弯矩、剪力图得绘制
GUI:
ﻫGeneralPostproc—plotResult—ContourPlot—LineElementResult
弹出画单元结果得对话框,分别在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
另:
自定义截面梁剪力弯矩显示
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
ansys如何绘制弯矩图
Ansy中弯矩图,云图绘制总结
在回答别人问题时,利用前人得回复与总结,自己进行了总结改正,发表在这里,供各位参考
(1)ANSYS弯矩等可以直接标注在图上吗?
如何实现?
如果三维问题,在剖面上标出某一结构得轴心力、弯矩等,如何实现
(2)后处理图形,其等值线得数值能否直接标注在图上,而不就是采用图例得形式 ﻫ后处理结果往往用云图表示,下跟一图例表示数值大小,能够实现等值线直接标注在图上
回答 ﻫ(1) ﻫ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,imom,jmom ﻫ2、标注弯矩图
PLOTCTRLS>>NUMBERING〉〉SVALON即可在画出弯矩图得同时在图上标出弯矩值得大小 ﻫ3、调整弯矩图
如果弯矩图方向错误,则绘制弯矩图命令为
plls,imom,jmom,-1 ﻫ同一个节点处两边得单元内力有细微差别, ﻫ导致内力数字标注出现重影。
观察上面整体轴力图也可以发现, ﻫ一段一段得,好像马赛克,其实上面整体弯矩图也就是,不过不就是 ﻫ很明显罢了.这就是EULER-BEONOULI梁理论以及ANSYS输出定义造成 ﻫ得(详细原因就不展开了,瞧瞧梁理论得书与ANSYS得说明吧)。
为了修正重影与节点两边内力值不一样得问题,遍制了宏文件ITFAVG、MAC
命令文件内容如下:
!
—-———-—--——-————---—----—----——-——--—--——---———--—--——-—-—--——--—-—-—
!
宏:
ITFAVG、MAC(INTERNALFORCEAVERAGE MACRO)
!
获取线性单元内力,并对单元边界处得内力进行平衡
!
输入信息
!
内力类型:
MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ ﻫ*ASK,ITFTYPE,’PLEASEINPUTTHETYPEOFINTERNALFORCE','MMOMY’
!
需处理得单元包 ﻫ*ASK,EASSEMBLY,’PLEASE INPUTTHEPONENTNAMEOFELEMENTSTOBEPROCESSED!
’, ﻫ'EOUTER’
!
需处理得节点包
*ASK,NASSEMBLY,'PLEASEINPUTTHEPONENT NAME OF NODE TOBEPROCESSED!
',’NOU
TER’
!
无需处理得节点包 ﻫ*ASK,UNASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFTHEUNCHANGED NODE!
(NONEI
FTHERE’SNO SUCHPONENT)','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 ﻫ*ELSE
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 ﻫ*ENDIF
*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
!
ENDOFITFAVG、MAC
(2) ﻫ对体与面来说,ANSYS默认得结果输出格式就是云图格式,而这种彩色云图打印为黑白图像时对比很不明显,无法表达清楚,对于发表文章非常不便。
发文章所用得结果图最好就是等值线图,并且最好就是黑白得等值线图。
一般借用photoshop等第三方软件,很麻烦,效果也不好。
ansys本身也能实现这项功能。
1、将要输出得结果调出,这时为彩色云图; ﻫ2、将云图转换为等值线图得形式 ﻫGUI:
plotCtrls-〉DeviceOptions—>[/DEVI]中得vectormode选为on
命令:
/DEVICE,VECTOR,1
这时结果为彩色等值线,若直接输出,打印为黑白图像时仍然不清晰,为此需进行以下几步将图像转换为黑白形式; ﻫ3、将背景变为白色
命令:
jpgprf,500,100,1
/rep
4、对等值线中得等值线符号(图中为A,B,C等)得疏密进行调整 ﻫGUI:
plotCtrls—〉Style—〉Contours-〉ContoursLabeling在KeyVectormodecontour label中选中on every Nthelem,然后在N= 输入框中输入合适得数值,例如5,多试几次,直到疏密合适 ﻫ命令:
/clabel,1,5
5、将彩色等值线变为黑色 ﻫGUI:
plotCtrls—>Style—〉Colors—〉ContoursColors 将ItemsNumbered1,Items Numbered2等复选框中得颜色均选为黑色,图像即可变为黑白等值线图像 ﻫ命令:
/color,cntr,whit,1等等 ﻫ6、修改
在显示等值线时,图形中太多得标识字母(A、B等)使得等值线图完全瞧不清楚.请问有没有什么方法可以去掉图形得这些标识字母呢
在plotCtrls下面得style下面得contours下面得contourLabeling,在弹出得contourlabeling options里面得vectormodecontourlabels选off ﻫ后处理中画出了结构得等值线图,可否将等值线边得字母改为相应得数值标注
/PNUM,SVAL,1 ﻫ/replot
如果还想将等值线得字母去掉 ﻫ/PNUM,SVAL,1 ﻫ/CLABEL,1,-1
/replot
7、出图 ﻫGUI:
plotCtrls—〉CaptureImage
ANSYS中如何加弯矩或扭矩
A转矩一般有三种施加得方法:
ﻫ第一种,将矩转换成一对一对得力偶,直接施加在对应得节点上面、ﻫ第二种,在构件中心部位建立一个节点,定义为MASS21单元,然后跟其她受力节点藕荷,形成刚性区域,就就是用CERIG命令、然后直接加转矩到主节点,即中心节点上面
第三种,使用MPC184单元、就是在构件中心部位建立一个节点,跟其她受力节点分别形成多根刚性梁,,从而形成刚性面、最后也就是直接加载荷到中心节点上面,通过刚性梁来传递载荷、
上面三种方法计算得结果基本一致,我做过实验得、
只不过就是后两种情况都就是形成刚性区域,但就是CERIG命令就是要在小变形或者小旋转才能用,只支持静力,线形分析、
而第三种方法适用多种情况,不仅支持大应变,还支持非线形情况、
如果您需要例子,我下次在发给您瞧
A
思路1:
矩或扭矩说白了就就是矩,所谓矩就就是力与力臂得乘积。
施加矩可以等效为施加力;
思路2:
直接施加弯矩或扭矩,此时需要引入一个具有旋转自由度得节点;可以选择单元21,或者184
1。
将矩转换成一对得力偶,直接施加在对应得节点上面。
2。
在构件中心部位建立一个节点,定义为mass21单元,然后跟其她受力节点耦合,形成刚性区域,就就是用cerig命令。
然后直接加转矩到主节点,即中心节点上面.
3。
使用mpc184单元。
就是在构件中心部位建立一个节点,跟其她受力节点分别形成多根刚性梁,从而形成刚性面。
最后也就是直接加载荷到中心节点上面,通过刚性梁来传递载荷。
4.通过rbe3命令.该方法与方法2很接近。
5。
基于表面边界法:
主要通过定义一个接触表面与一个目标节点接触来实现,弯矩荷载可以通过在目标节点上用“F"命令施加。
对于方法1,通过转换为集中力或均布力,比如施加扭矩,把端面节点改成柱坐标,然后等效为施加环向得节点力;而施加弯矩,可以将力矩转化为端面得剪切均布力;但这种方法比较容易出现应力集中现象;
方法2,定义局部刚性区域,施加过程venture讲得很详细,这里就不在赘述.根据她得例子,我在下面给出了一段命令流。
该方法有个不足,它在端面额外得增加了一定得刚度,只能适用于小变形分析。
方法3,相对方法2来说,采用刚性梁单元,适用范围更广一些,对于大应变分析也能很好得适用。
但在小应变分析下,方法2与方法3没有什么区别。
方法4,定义一个主节点,施加了分布力面,应该说跟实际比较接近一点,但端面得结果好像不就是很理想,结果有点偏大,在远离端面处得位置跟实际很符合。
方法5,它具体得受力形式有如下两种:
刚性表面边界(Rigidsurfaceconstraint)-认为接触面就是刚性得,没有变形,与通过节点耦合命令CERIG比较相似;
分布力边界(Force—distributedconstraint)-允许接触面得变形,与边界定义命令RBE3相似。
使用这种方法,需要用KEYOPT
(2)=2打开接触单元得MPC(多点接触边界)算法
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,MJ ﻫ2、标注弯矩图 ﻫPLOTCTRLS>〉NUMBERING〉>SVAL ON即可在画出弯矩图得同时在图上标出弯矩值得大小
3、调整弯矩图
如果弯矩图方向错误,则绘制弯矩图命令为
plls,imom,jmom,—1 ﻫ同一个节点处两边得单元内力有细微差别, ﻫ导致内力数字标注出现重影.观察上面整体轴力图也可以发现, ﻫ一段一段得,好像马赛克,其实上面整体弯矩图也就是,不过不就是 ﻫ很明显罢了。
这就是EULER-BEONOULI梁理论以及ANSYS输出定义造成 ﻫ得(详细原因就不展开了,瞧瞧梁理论得书与ANSYS得说明吧)。
为了修正重影与节点两边内力值不一样得问题,遍制了宏文件ITFAVG、MAC
命令文件内容如下:
!
—--—-——--------——--———-———--——-——---—-———-——-——--—-———-—-—-————-——---
!
宏:
ITFAVG、MAC(INTERNALFORCE AVERAGEMACRO)
!
获取线性单元内力,并对单元边界处得内力进行平衡
!
输入信息 ﻫ!
内力类型:
MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ ﻫ*ASK,ITFTYPE,'PLEASEINPUT THETYPEOFINTERNALFORCE’,’MMOMY' ﻫ!
需处理得单元包 ﻫ*ASK,EASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFELEMENTSTOBEPROCESSED!
’, ﻫ’EOUTER' ﻫ!
需处理得节点包 ﻫ*ASK,NASSEMBLY,’PLEASE INPUTTHEPONENTNAMEOFNODE TOBEPROCESSED!
’,’NOU
TER'
!
无需处理得节点包 ﻫ*ASK,UNASSEMBLY,'PLEASE INPUTTHE PONENTNAME OFTHE UNCHANGEDNODE!
(NONEI
FTHERE’SNOSUCH PONENT)’,'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',T