(完整版)ANSYS中彎矩、剪力圖的繪制

1、ansys中如何生成命令流方法：GUI是：UtilityMenuFileWriteDBLogFile怎么用ansys繪制彎矩，剪力圖：GUI:GeneralPostproc->lotResult->ContourPlot->LineElementResult彈出畫(huà)單元結(jié)果的對(duì)話框，分別在Labi和Labj依次選取SMIS6和SMIS12（彎矩圖）、SMIS1和SMIS7（軸力圖）、SMIS2和SMIS8（剪力圖）!建立單元表ETABLE,NI,SMISC,1!單元I點(diǎn)軸力ETABLE,NJ,SMISC,7!單元J點(diǎn)軸力ETABLE,QI,SMISC,2!單元I點(diǎn)剪力ETABLE

2、,QJ,SMISC,8!單元J點(diǎn)剪力ETABLE,MI,SMISC,6!單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12!單元J點(diǎn)彎矩!更新單元表ETABLE,REFL!畫(huà)軸力分布圖/TITLE,AxialforcediagramPLLS,NI,NJ,1.0,0/image,save,Axial_force_%T%,jpg!畫(huà)剪力分布圖/TITLE,ShearingforcediagramPLLS,QI,QJ,1.0,0/image,save,Shearing_force_%T%,jpg!畫(huà)彎矩分布圖/TITLE,BendingmomentdiagramPLLS,MI,MJ,-0.8,0/im

3、age,save,Bending_moment_%T%,jpgANSYS中彎矩、剪力圖的繪制GUI:GeneralPostproc-plotResult-ContourPlot-LineElementResult彈出畫(huà)單元結(jié)果的對(duì)話框，分別在Labi和Labj依次選取SMIS6和SMIS12（彎矩圖）、SMIS1和SMIS7（軸力圖）、SMIS2和SMIS8（剪力圖）!建立單元表ETABLE,NI,SMISC,1!單元I點(diǎn)軸力ETABLE,NJ,SMISC,7!單元J點(diǎn)軸力ETABLE,QI,SMISC,2!單元I點(diǎn)剪力ETABLE,QJ,SMISC,8!單元J點(diǎn)剪力ETABLE,MI,SMI

4、SC,6!單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12!單元J點(diǎn)彎矩!更新單元表ETABLE,REFL!畫(huà)軸力分布圖/TITLE,AxialforcediagramPLLS,NI,NJ,1.0,0/image,save,Axial_force_%T%,jpg!畫(huà)剪力分布圖/TITLE,ShearingforcediagramPLLS,QI,QJ,1.0,0/image,save,Shearing_force_%T%,jpg!畫(huà)彎矩分布圖/TITLE,BendingmomentdiagramPLLS,MI,MJ,-0.8,0/image,save,Bending_moment_%T%,jpg

5、另：自定義截面梁剪力彎矩顯示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!讀入截面文

6、件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,

7、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中彎矩圖，云圖繪制總結(jié)在回答別人問(wèn)題時(shí)，利用前人的回復(fù)和總結(jié)，自己進(jìn)行了總結(jié)改正，發(fā)表在這里，供各位參考(1) ANSYS彎矩等可以直接標(biāo)注在圖上嗎

8、？如何實(shí)現(xiàn)？如果三維問(wèn)題，在剖面上標(biāo)出某一結(jié)構(gòu)的軸心力、彎矩等，如何實(shí)現(xiàn)(2) 后處理圖形，其等值線的數(shù)值能否直接標(biāo)注在圖上，而不是采用圖例的形式后處理結(jié)果往往用云圖表示，下跟一圖例表示數(shù)值大小，能夠?qū)崿F(xiàn)等值線直接標(biāo)注在圖上回答(1)1. 繪制彎矩圖建立彎矩單元表。例如梁?jiǎn)卧猧節(jié)點(diǎn)單元表名稱為imom,j節(jié)點(diǎn)單元表名稱為jmom,ETABLE,NI,SMISC,1!單元I點(diǎn)軸力ETABLE,NJ,SMISC,7!單元J點(diǎn)軸力ETABLE,QI,SMISC,2!單元I點(diǎn)剪力ETABLE,QJ,SMISC,8!單元J點(diǎn)剪力ETABLE,MI,SMISC,6!單元I點(diǎn)彎矩ETABLE,MJ,SMIS

9、C,12!單元J點(diǎn)彎矩plls,imom,jmom2. 標(biāo)注彎矩圖PLOTCTRLSNUMBERINGSVALON即可在畫(huà)出彎矩圖的同時(shí)在圖上標(biāo)出彎矩值的大小3. 調(diào)整彎矩圖如果彎矩圖方向錯(cuò)誤，則繪制彎矩圖命令為plls,imom,jmom,-1同一個(gè)節(jié)點(diǎn)處兩邊的單元內(nèi)力有細(xì)微差別，導(dǎo)致內(nèi)力數(shù)字標(biāo)注出現(xiàn)重影。觀察上面整體軸力圖也可以發(fā)現(xiàn)，一段一段的，好像馬賽克，其實(shí)上面整體彎矩圖也是，不過(guò)不是很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成的(詳細(xì)原因就不展開(kāi)了，看看梁理論的書(shū)和ANSYS的說(shuō)明吧)。為了修正重影和節(jié)點(diǎn)兩邊內(nèi)力值不一樣的問(wèn)題，遍制了宏文件ITFAVG

10、.MAC命令文件內(nèi)容如下：!宏：ITFAVG.MAC(INTERNALFORCEAVERAGEMACRO)!獲取線性單元內(nèi)力，并對(duì)單元邊界處的內(nèi)力進(jìn)行平衡!輸入信息!內(nèi)力類型：MFORX,MFORY,MFOR乙MMOMX,MMOMY,MMOMZ*ASK,ITFTYPE,PLEASEINPUTTHETYPEOFINTERNALFORCE,MMOMY!需處理的單元包*ASK,EASSEMBLY,PLEASEINPUTTHECOMPONENTNAMEOFELEMENTSTOBEPROCESSED!,EOUTER!需處理的節(jié)點(diǎn)包*ASK,NASSEMBLY,PLEASEINPUTTHECOMPONEN

11、TNAMEOFNODETOBEPROCESSED!,NOUTER!無(wú)需處理的節(jié)點(diǎn)包*ASK,UNASSEMBLY,PLEASEINPUTTHECOMPONENTNAMEOFTHEUNCHANGEDNODE!(NONEIFTHERESNOSUCHCOMPONENT),NONE/POST1!輸入信息：內(nèi)力類型，欲處理單元的集合，欲處理節(jié)點(diǎn)的集合!ITFTYPE=MMOMY!EASSEMBLY=EOUTER!NASSEMBLY=NOUTER!按內(nèi)力類型確定ANSYS輸出信息SMISC的編號(hào)*IF,ITFTYPE,EQ,MFORX,THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTY

12、PE,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!對(duì)不需平均的節(jié)點(diǎn)進(jìn)行處理*IF,UNASSEMBLY,NE,NONE,THEN!選出不進(jìn)行處理的節(jié)點(diǎn)包并獲取不進(jìn)行處理節(jié)點(diǎn)的數(shù)

13、目CMSEL,S,UNASSEMBLY*GET,UNNODNUM,NODE,0,COUNT!定義長(zhǎng)度為UNNODNUM的數(shù)組（UNNOD）,以存放選中單元的單元編號(hào)*DIM,UNNOD,ARRAY,UNNODNUM!將選中單元的編號(hào)按順序存入數(shù)組UNNOD*DO,I,0,UNNODNUM-1,1UNNOD(I+1)=NDNEXT(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!選出所需的單元和節(jié)點(diǎn)包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!獲得當(dāng)前選中單元總數(shù)(存入變量SELELENUM)*GET,SELELENUM,ELEM,0,COUNT!定義長(zhǎng)度為

14、SELELENUM的數(shù)組(ELENUM),以存放選中單元的單元編號(hào)*DIM,ELENUM,ARRAY,SELELENUM!將選中單元的編號(hào)按順序存入數(shù)組ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNEXT(I)*ENDDO!獲得當(dāng)前選中節(jié)點(diǎn)總數(shù)(存入變量SELNODNUM)*GET,SELNODNUM,NODE,0,COUNT!定義長(zhǎng)度為SELNODNUM的數(shù)組(NODNUM),以存放選中單元的單元編號(hào)*DIM,NODNUM,ARRAY,SELNODNUM!將選中單元的編號(hào)按順序存入數(shù)組NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I

15、+1)=NDNEXT(I)*ENDDO!定義所需的線性單元內(nèi)力ETABLE，節(jié)點(diǎn)I的內(nèi)力存入數(shù)組ITNFI,!節(jié)點(diǎn)J的內(nèi)力存入數(shù)組ITNFJETABLE,ITNFI,SMISC,ITFINUMETABLE,ITNFJ,SMISC,ITFJNUM!定義所需的結(jié)果數(shù)組，并將其置零ETABLE,ITNFINEO,SMISC,5SADD,ITNFINEO,ITNFI,1ETABLE,ITNFJNEO,SMISC,11SADD,ITNFJNEO,ITNFJ,1*DO,K,1,SELNODNUM,1!處理不需平均的節(jié)點(diǎn)INDEX=0*IF,UNNODNUM,GE,1,THEN*DO,J,1,UNNODNU

16、M*IF,NODNUM(K),EQ,UNNOD(J),THENINDEX=1*ELSE*ENDIF*ENDDO*ELSE*ENDIF*DO,J,1,SELELENUM,1!選出和節(jié)點(diǎn)K相連的線性單元中，I節(jié)點(diǎn)(對(duì)線性單元而言)為節(jié)點(diǎn)K的單元編號(hào)*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!選出和節(jié)點(diǎn)K相連的線性單元中，J節(jié)點(diǎn)(對(duì)線性單元而言)為節(jié)點(diǎn)K的單元編號(hào)*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THENELEJ=E

17、LENUM(J)*EXIT*ELSE*ENDIF*ENDDO*IF,INDEX,EQ,0,THEN*IF,ELEJ,NE,0,THEN!有可能出現(xiàn)ELEJ為0的情況!取出I節(jié)點(diǎn)為節(jié)點(diǎn)K的單元的I節(jié)點(diǎn)端的內(nèi)力放入?yún)?shù)ETELEI*GET,ETELEI,ELEM,ELEI,ETAB,ITNFI!取出J節(jié)點(diǎn)為節(jié)點(diǎn)K的單元的J節(jié)點(diǎn)端的內(nèi)力放入?yún)?shù)ETELEJ*GET,ETELEJ,ELEM,ELEJ,ETAB,ITNFJ!平均節(jié)點(diǎn)K的單元的I節(jié)點(diǎn)端的內(nèi)力和節(jié)點(diǎn)K的單元的J節(jié)點(diǎn)端的內(nèi)力ETAVE=(ETELEI+ETELEJ)/2!將平均后的內(nèi)力存入結(jié)果數(shù)組中DETAB,ELEI,ITNFINEO,ET

18、AVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!ENDOFITFAVG.MAC(2)對(duì)體和面來(lái)說(shuō)，ANSYS默認(rèn)的結(jié)果輸出格式是云圖格式，而這種彩色云圖打印為黑白圖像時(shí)對(duì)比很不明顯，無(wú)法表達(dá)活楚，對(duì)于發(fā)表文章非常不便。發(fā)文章所用的結(jié)果圖最好是等值線圖，并且最好是黑白的等值線圖。一般借用photoshop等第三方軟件，很麻煩，效果也不好。ansys本身也能實(shí)現(xiàn)這項(xiàng)功能。1. 將要輸出的結(jié)果調(diào)出，這時(shí)為彩色云圖；2. 將云圖轉(zhuǎn)換為等值線圖的形式GUI:plo

19、tCtrlsDeviceOptions/DEVI中的vectormode選為on命令：/DEVICE,VECTOR,1這時(shí)結(jié)果為彩色等值線，若直接輸出，打印為黑白圖像時(shí)仍然不活晰，為此需進(jìn)行以下幾步將圖像轉(zhuǎn)換為黑白形式；3. 將背景變?yōu)榘咨睿簀pgprf,500,100,1/rep4. 對(duì)等值線中的等值線符號(hào)（圖中為A,B,C等）的疏密進(jìn)行調(diào)整GUI:plotCtrlsStyleContourLContoursLabeling在KeyVectormodecontourlabel中選中oneveryNthelemi,然后在N=輸入框中輸入合適的數(shù)值，例如5,多試幾次，直到疏密合適命令：/cl

20、abel,1,55. 將彩色等值線變?yōu)楹谏獹UI:plotCtrlsStyleColorsContoursColors將ItemsNumbered1,ItemsNumbered2等復(fù)選框中的顏色均選為黑色，圖像即可變?yōu)楹诎椎戎稻€圖像命令：/color,cntr,whit,1等等6修改在顯示等值線時(shí)，圖形中太多的標(biāo)識(shí)字母（A、B等）使得等值線圖完全看不活楚。請(qǐng)問(wèn)有沒(méi)有什么方法可以去掉圖形的這些標(biāo)識(shí)字母呢在plotCtrls下面的style下面的contours下面的contourLabeling,在彈出的contourlabelingoptions里面的vectormodecontourlabe

21、lsoff后處理中畫(huà)出了結(jié)構(gòu)的等值線圖，可否將等值線邊的字母改為相應(yīng)的數(shù)值標(biāo)注/PNUM,SVAL,1/replot如果還想將等值線的字母去掉/PNUM,SVAL,1/CLABEL,1,-1/replot7.出圖GUI:plotCtrlsCaptureImageANSYS中如何加彎矩或扭矩A轉(zhuǎn)矩一般有三種施加的方法：第一種，將矩轉(zhuǎn)換成一對(duì)一對(duì)的力偶，直接施加在對(duì)應(yīng)的節(jié)點(diǎn)上面.第二種，在構(gòu)件中心部位建立一個(gè)節(jié)點(diǎn)，定義為MASS21單元，然后跟其他受力節(jié)點(diǎn)藕荷,形成剛性區(qū)域，就是用CERIG命令.然后直接加轉(zhuǎn)矩到主節(jié)點(diǎn),即中心節(jié)點(diǎn)上面第三種，使用MPC184單元.是在構(gòu)件中心部位建立一個(gè)節(jié)點(diǎn)，跟其

22、他受力節(jié)點(diǎn)分別形成多根剛性梁，,從而形成剛性面.最后也是直接加載荷到中心節(jié)點(diǎn)上面，通過(guò)剛性梁來(lái)傳遞載荷.上面三種方法計(jì)算的結(jié)果基本一致,我做過(guò)實(shí)驗(yàn)的.只不過(guò)是后兩種情況都是形成剛性區(qū)域，但是CERIG命令是要在小變形或者小旋轉(zhuǎn)才能用，只支持靜力,線形分析.而第三種方法適用多種情況,不僅支持大應(yīng)變,還支持非線形情況.如果你需要例子，我下次在發(fā)給你看A思路1:矩或扭矩說(shuō)白了就是矩，所謂矩就是力和力臂的乘積。施加矩可以等效為施加力；思路2:直接施加彎矩或扭矩，此時(shí)需要引入一個(gè)具有旋轉(zhuǎn)自由度的節(jié)點(diǎn)；可以選擇單元21,或者1841 .將矩轉(zhuǎn)換成一對(duì)的力偶，直接施加在對(duì)應(yīng)的節(jié)點(diǎn)上面。2 .在構(gòu)件中心部位建

23、立一個(gè)節(jié)點(diǎn)，定義為mass21單元，然后跟其他受力節(jié)點(diǎn)耦合，形成剛性區(qū)域，就是用cerig命令。然后直接加轉(zhuǎn)矩到主節(jié)點(diǎn)，即中心節(jié)點(diǎn)上面。3 .使用mpc184單元。是在構(gòu)件中心部位建立一個(gè)節(jié)點(diǎn)，跟其他受力節(jié)點(diǎn)分別形成多根剛性梁，從而形成剛性面。最后也是直接加載荷到中心節(jié)點(diǎn)上面，通過(guò)剛性梁來(lái)傳遞載荷。4 .通過(guò)rbe3命令。該方法與方法2很接近。5.基丁表面邊界法：主要通過(guò)定義一個(gè)接觸表面和一個(gè)目標(biāo)節(jié)點(diǎn)接觸來(lái)實(shí)現(xiàn)，彎矩荷載可以通過(guò)在目標(biāo)節(jié)點(diǎn)上用F”命令施加。對(duì)丁方法1,通過(guò)轉(zhuǎn)換為集中力或均布力，比如施加扭矩，把端面節(jié)點(diǎn)改成柱坐標(biāo)，然后等效為施加環(huán)向的節(jié)點(diǎn)力；而施加彎矩，可以將力矩轉(zhuǎn)化為端面的剪切

24、均布力；但這種方法比較容易出現(xiàn)應(yīng)力集中現(xiàn)象；方法2,定義局部剛性區(qū)域，施加過(guò)程venture講的很詳細(xì)，這里就不在贅述。根據(jù)他的例子，我在下面給出了一段命令流。該方法有個(gè)不足，它在端面額外的增加了一定的剛度，只能適用丁小變形分析。方法3,相對(duì)方法2來(lái)說(shuō)，采用剛性梁?jiǎn)卧m用范圍更廣一些，對(duì)丁大應(yīng)變分析也能很好的適用。但在小應(yīng)變分析下，方法2和方法3沒(méi)有什么區(qū)別。方法4,定義一個(gè)主節(jié)點(diǎn)，施加了分布力面，應(yīng)該說(shuō)跟實(shí)際比較接近一點(diǎn)，但端面的結(jié)果好像不是很理想，結(jié)果有點(diǎn)偏大，在遠(yuǎn)離端面處的位置跟實(shí)際很符合。方法5,它具體的受力形式有如下兩種：剛性表面邊界(Rigidsurfaceconstraint)

25、認(rèn)為接觸面是剛性的，沒(méi)有變形，和通過(guò)節(jié)點(diǎn)耦合命令CERIG比較相似；分布力邊界(Force-distributedconstraint)允許接觸面的變形，和邊界定義命令RBE3相似。使用這種方法，需要用KEYOPT(2)=2打開(kāi)接觸單元的MPC(多點(diǎn)接觸邊界)算法ANSYS繪制彎矩、剪力、軸力圖1. 繪制彎矩圖建立彎矩單元表。例如梁?jiǎn)卧猧節(jié)點(diǎn)單元表名稱為imom,j節(jié)點(diǎn)單元表名稱為jmom,ETABLE,NI,SMISC,1!單元I點(diǎn)軸力ETABLE,NJ,SMISC,7!單元J點(diǎn)軸力ETABLE,QI,SMISC,2!單元I點(diǎn)剪力ETABLE,QJ,SMISC,8!單元J點(diǎn)剪力ETABLE,

26、MI,SMISC,6!單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12!單元J點(diǎn)彎矩plls,MI,MJ2. 標(biāo)注彎矩圖PLOTCTRLSNUMBERINGSVALON即可在畫(huà)出彎矩圖的同時(shí)在圖上標(biāo)出彎矩值的大小3. 調(diào)整彎矩圖如果彎矩圖方向錯(cuò)誤，則繪制彎矩圖命令為plls,imom,jmom,-1同一個(gè)節(jié)點(diǎn)處兩邊的單元內(nèi)力有細(xì)微差別，導(dǎo)致內(nèi)力數(shù)字標(biāo)注出現(xiàn)重影。觀察上面整體軸力圖也可以發(fā)現(xiàn)，一段一段的，好像馬賽克，其實(shí)上面整體彎矩圖也是，不過(guò)不是很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成的(詳細(xì)原因就不展開(kāi)了，看看梁理論的書(shū)和ANSYS的說(shuō)明吧)。為了修正重影

27、和節(jié)點(diǎn)兩邊內(nèi)力值不一樣的問(wèn)題，遍制了宏文件ITFAVG.MAC命令文件內(nèi)容如下：!宏：ITFAVG.MAC(INTERNALFORCEAVERAGEMACRO)!獲取線性單元內(nèi)力，并對(duì)單元邊界處的內(nèi)力進(jìn)行平衡!輸入信息!內(nèi)力類型：MFORX,MFORY,MFOR乙MMOMX,MMOMY,MMOMZ*ASK,ITFTYPE,PLEASEINPUTTHETYPEOFINTERNALFORCE,MMOMY!需處理的單元包*ASK,EASSEMBLY,PLEASEINPUTTHECOMPONENTNAMEOFELEMENTSTOBEPROCESSED!,EOUTER!需處理的節(jié)點(diǎn)包*ASK,NASSE

28、MBLY,PLEASEINPUTTHECOMPONENTNAMEOFNODETOBEPROCESSED!,NOUTER!無(wú)需處理的節(jié)點(diǎn)包*ASK,UNASSEMBLY,PLEASEINPUTTHECOMPONENTNAMEOFTHEUNCHANGEDNODE!(NONEIFTHERESNOSUCHCOMPONENT),NONE/POST1!輸入信息：內(nèi)力類型，欲處理單元的集合，欲處理節(jié)點(diǎn)的集合!ITFTYPE=MMOMY!EASSEMBLY=EOUTER!NASSEMBLY=NOUTER!按內(nèi)力類型確定ANSYS輸出信息SMISC的編號(hào)*IF,ITFTYPE,EQ,MFORX,THENITFI

29、NUM=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!對(duì)不需平均的節(jié)點(diǎn)進(jìn)行處理*IF,UNASSEMBLY,NE,NONE,T

30、HEN!選出不進(jìn)行處理的節(jié)點(diǎn)包并獲取不進(jìn)行處理節(jié)點(diǎn)的數(shù)目CMSEL,S,UNASSEMBLY*GET,UNNODNUM,NODE,0,COUNT!定義長(zhǎng)度為UNNODNUM的數(shù)組(UNNOD),以存放選中單元的單元編號(hào)*DIM,UNNOD,ARRAY,UNNODNUM!將選中單元的編號(hào)按順序存入數(shù)組UNNOD*DO,I,0,UNNODNUM-1,1UNNOD(I+1)=NDNEXT(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!選出所需的單元和節(jié)點(diǎn)包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!獲得當(dāng)前選中單元總數(shù)(存入變量SELELENUM)*GET,S

31、ELELENUM,ELEM,0,COUNT!定義長(zhǎng)度為SELELENUM的數(shù)組(ELENUM),以存放選中單元的單元編號(hào)*DIM,ELENUM,ARRAY,SELELENUM!將選中單元的編號(hào)按順序存入數(shù)組ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNEXT(I)*ENDDO!獲得當(dāng)前選中節(jié)點(diǎn)總數(shù)(存入變量SELNODNUM)*GET,SELNODNUM,NODE,0,COUNT!定義長(zhǎng)度為SELNODNUM的數(shù)組(NODNUM),以存放選中單元的單元編號(hào)*DIM,NODNUM,ARRAY,SELNODNUM!將選中單元的編號(hào)按順序存入數(shù)組NODNUM*D

32、O,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNEXT(I)*ENDDO!定義所需的線性單元內(nèi)力ETABLE，節(jié)點(diǎn)I的內(nèi)力存入數(shù)組ITNFI,!節(jié)點(diǎn)J的內(nèi)力存入數(shù)組ITNFJETABLE,ITNFI,SMISC,ITFINUMETABLE,ITNFJ,SMISC,ITFJNUM!定義所需的結(jié)果數(shù)組，并將其置零ETABLE,ITNFINEO,SMISC,5SADD,ITNFINEO,ITNFI,1ETABLE,ITNFJNEO,SMISC,11SADD,ITNFJNEO,ITNFJ,1*DO,K,1,SELNODNUM,1!處理不需平均的節(jié)點(diǎn)INDEX=0*IF,UNNODN

33、UM,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!選出和節(jié)點(diǎn)K相連的線性單元中，I節(jié)點(diǎn)(對(duì)線性單元而言)為節(jié)點(diǎn)K的單元編號(hào)*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!選出和節(jié)點(diǎn)K相連的線性單元中，J節(jié)點(diǎn)(對(duì)線性單元而言)為節(jié)點(diǎn)K的單元編號(hào)*IF,NELEM(ELENUM(J)

34、,2),EQ,NODNUM(K),THENELEJ=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*IF,INDEX,EQ,0,THEN*IF,ELEJ,NE,0,THEN!有可能出現(xiàn)ELEJ為0的情況!取出I節(jié)點(diǎn)為節(jié)點(diǎn)K的單元的I節(jié)點(diǎn)端的內(nèi)力放入?yún)?shù)ETELEI*GET,ETELEI,ELEM,ELEI,ETAB,ITNFI!取出J節(jié)點(diǎn)為節(jié)點(diǎn)K的單元的J節(jié)點(diǎn)端的內(nèi)力放入?yún)?shù)ETELEJ*GET,ETELEJ,ELEM,ELEJ,ETAB,ITNFJ!平均節(jié)點(diǎn)K的單元的I節(jié)點(diǎn)端的內(nèi)力和節(jié)點(diǎn)K的單元的J節(jié)點(diǎn)端的內(nèi)力ETAVE=(ETELEI+ETELEJ)/2!將平均后的內(nèi)力存入

35、結(jié)果數(shù)組中DETAB,ELEI,ITNFINEO,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!ENDOFITFAVG.MAC（2）對(duì)體和面來(lái)說(shuō)，ANSYS默認(rèn)的結(jié)果輸出格式是云圖格式，而這種彩色云圖打印為黑白圖像時(shí)對(duì)比很不明顯，無(wú)法表達(dá)活楚，對(duì)于發(fā)表文章非常不便。發(fā)文章所用的結(jié)果圖最好是等值線圖，并且最好是黑白的等值線圖。一般借用photoshop等第三方軟件，很麻煩，效果也不好。ansys本身也能實(shí)現(xiàn)這項(xiàng)功能。1. 將要輸出的結(jié)果調(diào)出，這時(shí)為彩

36、色云圖；2. 將云圖轉(zhuǎn)換為等值線圖的形式GUI:plotCtrlsDeviceOptions/DEVI中的vectormode選為on命令：/DEVICE,VECTOR,1這時(shí)結(jié)果為彩色等值線，若直接輸出，打印為黑白圖像時(shí)仍然不活晰，為此需進(jìn)行以下幾步將圖像轉(zhuǎn)換為黑白形式；3. 將背景變?yōu)榘咨睿簀pgprf,500,100,1/rep4. 對(duì)等值線中的等值線符號(hào)（圖中為A,B,C等）的疏密進(jìn)行調(diào)整GUI:plotCtrlsStyleContourLContoursLabeling在KeyVectormodecontourlabel中選中oneveryNthelemi,然后在N=輸入框中輸入

37、合適的數(shù)值，例如5,多試幾次，直到疏密合適命令：/clabel,1,55. 將彩色等值線變?yōu)楹谏獹UI:plotCtrlsStyleColorsContoursColors將ItemsNumbered1,ItemsNumbered2等復(fù)選框中的顏色均選為黑色，圖像即可變?yōu)楹诎椎戎稻€圖像命令：/color,cntr,whit,1等等6修改在顯示等值線時(shí)，圖形中太多的標(biāo)識(shí)字母（A、B等）使得等值線圖完全看不活楚。請(qǐng)問(wèn)有沒(méi)有什么方法可以去掉圖形的這些標(biāo)識(shí)字母呢在plotCtrls下面的style下面的contours下面的contourLabeling,在彈出的contourlabelingopti

38、ons里面的vectormodecontourlabelsoff后處理中畫(huà)出了結(jié)構(gòu)的等值線圖，可否將等值線邊的字母改為相應(yīng)的數(shù)值標(biāo)注/PNUM,SVAL,1/replot如果還想將等值線的字母去掉/PNUM,SVAL,1/CLABEL,1,-1/replot7.出圖GUI:plotCtrlsCaptureImage【原創(chuàng)】下面部分為作者自己總結(jié)對(duì)丁beam3單元睥-1.!5皿野ISMISC-3|育炬MM0X1Z-W其他的單元類型，可以參考ansys幫助文檔中給出的相關(guān)單元的參數(shù)。附上一個(gè)小例子，下圖所示為一懸臂梁，長(zhǎng)100,左端完全約束，上面加載均勻載荷q=100,繪制這個(gè)梁的剪力和彎矩圖。L

39、；L上&Ipb扁hI彎矩圖TIMESTRESSANSYSSUB-ITIME=1HIMJHIKT=-500000ELEH-1MAX-A5E-07ELEM-107UY1920101&43普9-SOOaQO-308639-277773-166657T555&-4444W-33333J-E22222-112111.155EQ7剪力圖1. 繪制彎矩圖建立彎矩單元表。例如梁?jiǎn)卧猧節(jié)點(diǎn)單元表名稱為imom,j節(jié)點(diǎn)單元表名稱為jmom,ETABLE,NI,SMISC,1ETABLE,NJ,SMISC,7ETABLE,QI,SMISC,2ETABLE,QJ,SMISC,8ETABLE,MI,SMISC,6!單元

40、I點(diǎn)軸力!單元J點(diǎn)軸力!單元I點(diǎn)剪力!單元J點(diǎn)剪力!單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12!單元J點(diǎn)彎矩plls,MI,MJ2. 標(biāo)注彎矩圖PLOTCTRLSNUMBERINGSVALON即可在畫(huà)出彎矩圖的同時(shí)在圖上標(biāo)出彎矩值的大小3. 調(diào)整彎矩圖如果彎矩圖方向錯(cuò)誤，則繪制彎矩圖命令為plls,imom,jmom,-1同一個(gè)節(jié)點(diǎn)處兩邊的單元內(nèi)力有細(xì)微差別，導(dǎo)致內(nèi)力數(shù)字標(biāo)注出現(xiàn)重影。觀察上面整體軸力圖也可以發(fā)現(xiàn)，一段一段的，好像馬賽克，其實(shí)上面整體彎矩圖也是，不過(guò)不是很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成的(詳細(xì)原因就不展開(kāi)了，看看梁理論的書(shū)和AN

41、SYS的說(shuō)明吧)。為了修正重影和節(jié)點(diǎn)兩邊內(nèi)力值不一樣的問(wèn)題，遍制了宏文件ITFAVG.MAC命令文件內(nèi)容如下：!宏：ITFAVG.MAC(INTERNALFORCEAVERAGEMACRO)!獲取線性單元內(nèi)力，并對(duì)單元邊界處的內(nèi)力進(jìn)行平衡!輸入信息!內(nèi)力類型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ*ASK,ITFTYPE,PLEASEINPUTTHETYPEOFINTERNALFORCE,MMOMY!需處理的單元包*ASK,EASSEMBLY,PLEASEINPUTTHECOMPONENTNAMEOFELEMENTSTOBEPROCESSED!,EOUTER!需

42、處理的節(jié)點(diǎn)包*ASK,NASSEMBLY,PLEASEINPUTTHECOMPONENTNAMEOFNODETOBEPROCESSED!,NOUTER!無(wú)需處理的節(jié)點(diǎn)包*ASK,UNASSEMBLY,PLEASEINPUTTHECOMPONENTNAMEOFTHEUNCHANGEDNODE!(NONEIFTHERESNOSUCHCOMPONENT),NONE/POST1!輸入信息：內(nèi)力類型，欲處理單元的集合，欲處理節(jié)點(diǎn)的集合!ITFTYPE=MMOMY!EASSEMBLY=EOUTER!NASSEMBLY=NOUTER!按內(nèi)力類型確定ANSYS輸出信息SMISC的編號(hào)*IF,ITFTYPE,E

43、Q,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!對(duì)不需平均的節(jié)點(diǎn)進(jìn)行處理*IF,UNAS

44、SEMBLY,NE,NONE,THEN!選出不進(jìn)行處理的節(jié)點(diǎn)包并獲取不進(jìn)行處理節(jié)點(diǎn)的數(shù)目CMSEL,S,UNASSEMBLY*GET,UNNODNUM,NODE,0,COUNT!定義長(zhǎng)度為UNNODNUM的數(shù)組(UNNOD),以存放選中單元的單元編號(hào)*DIM,UNNOD,ARRAY,UNNODNUM!將選中單元的編號(hào)按順序存入數(shù)組UNNOD*DO,I,0,UNNODNUM-1,1UNNOD(I+1)=NDNEXT(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!選出所需的單元和節(jié)點(diǎn)包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!獲得當(dāng)前選中單元總數(shù)(存入變量

45、SELELENUM)*GET,SELELENUM,ELEM,0,COUNT!定義長(zhǎng)度為SELELENUM的數(shù)組(ELENUM),以存放選中單元的單元編號(hào)*DIM,ELENUM,ARRAY,SELELENUM!將選中單元的編號(hào)按順序存入數(shù)組ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNEXT(I)*ENDDO!獲得當(dāng)前選中節(jié)點(diǎn)總數(shù)(存入變量SELNODNUM)*GET,SELNODNUM,NODE,0,COUNT!定義長(zhǎng)度為SELNODNUM的數(shù)組(NODNUM),以存放選中單元的單元編號(hào)*DIM,NODNUM,ARRAY,SELNODNUM!將選中單元的編

46、號(hào)按順序存入數(shù)組NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNEXT(I)*ENDDO!定義所需的線性單元內(nèi)力ETABLE，節(jié)點(diǎn)I的內(nèi)力存入數(shù)組ITNFI,!節(jié)點(diǎn)J的內(nèi)力存入數(shù)組ITNFJETABLE,ITNFI,SMISC,ITFINUMETABLE,ITNFJ,SMISC,ITFJNUM!定義所需的結(jié)果數(shù)組，并將其置零ETABLE,ITNFINEO,SMISC,5SADD,ITNFINEO,ITNFI,1ETABLE,ITNFJNEO,SMISC,11SADD,ITNFJNEO,ITNFJ,1*DO,K,1,SELNODNUM,1!處理不需平均的節(jié)點(diǎn)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!選