ANSYS中彎矩剪力圖的繪制92290

1、ansys中如何生成命令流方法：GUI是：Utility MenuFileWrite DB Log File怎么用ansys繪制彎矩，剪力圖：GUI: General Postproc->lot Result->Contour Plot->Line Element Result彈出畫單元結(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

2、 !單元I點(diǎn)剪力ETABLE,QJ,SMISC,8 !單元J點(diǎn)剪力ETABLE,MI,SMISC,6 !單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12 !單元J點(diǎn)彎矩 ! 更新單元表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 mo

3、ment diagramPLLS,MI,MJ,-0.8,0/image,save,Bending_moment_%T%,jpgANSYS中彎矩、剪力圖的繪制GUI:General Postproc-plot Result-Contour Plot-Line Element Result彈出畫單元結(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 !單元

4、I點(diǎn)剪力ETABLE,QJ,SMISC,8 !單元J點(diǎn)剪力ETABLE,MI,SMISC,6 !單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12 !單元J點(diǎn)彎矩 ! 更新單元表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

5、 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

6、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,1l

7、esize,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中彎矩圖，云圖繪制總結(jié)在回答別

8、人問題時(shí)，利用前人的回復(fù)和總結(jié)，自己進(jìn)行了總結(jié)改正，發(fā)表在這里，供各位參考（1）ANSYS彎矩等可以直接標(biāo)注在圖上嗎？如何實(shí)現(xià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)剪力ETA

9、BLE,QJ,SMISC,8 !單元J點(diǎn)剪力ETABLE,MI,SMISC,6 !單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12 !單元J點(diǎn)彎矩plls,imom,jmom2.標(biāo)注彎矩圖PLOTCTRLSNUMBERINGSVAL ON即可在畫出彎矩圖的同時(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í)上面整體彎矩圖也是，不過不是很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成的（詳

10、細(xì)原因就不展開了，看看梁理論的書和ANSYS的說明吧）。為了修正重影和節(jié)點(diǎn)兩邊內(nèi)力值不一樣的問題，遍制了宏文件ITFAVG.MAC命令文件內(nèi)容如下：!-!宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO）!獲取線性單元內(nèi)力，并對(duì)單元邊界處的內(nèi)力進(jìn)行平衡!輸入信息!內(nèi)力類型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ*ASK,ITFTYPE,PLEASE INPUT THE TYPE OF INTERNAL FORCE,MMOMY!需處理的單元包*ASK,EASSEMBLY,PLEASE INPUT THE COMPONENT NAM

11、E OF ELEMENTS TO BE PROCESSED!,EOUTER!需處理的節(jié)點(diǎn)包*ASK,NASSEMBLY,PLEASE INPUT THE COMPONENT NAME OF NODE TO BE PROCESSED!,NOUTER!無需處理的節(jié)點(diǎn)包*ASK,UNASSEMBLY,PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NODE!(NONE IF THERES NO SUCH COMPONENT),NONE/POST1!輸入信息：內(nèi)力類型，欲處理單元的集合，欲處理節(jié)點(diǎn)的集合!ITFTYPE=MMOMY!EASSEMBLY=

12、EOUTER!NASSEMBLY=NOUTER!按內(nèi)力類型確定ANSYS輸出信息SMISC的編號(hào)*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,

13、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ù)目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!

14、選出所需的單元和節(jié)點(diǎn)包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!獲得當(dāng)前選中單元總數(shù)（存入變量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)度為SELNODN

15、UM的數(shù)組（NODNUM），以存放選中單元的單元編號(hào)*DIM,NODNUM,ARRAY,SELNODNUM!將選中單元的編號(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,ITNFJNE

16、O,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!選出和節(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*END

17、DO*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=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,EL

18、EJ,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,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!END OF ITFAVG.MAC（2）對(duì)體和面來說，ANSYS默認(rèn)的結(jié)果輸出格式是云圖格式，而這種彩色云圖打印為黑白圖像時(shí)對(duì)比很不明顯，無法表達(dá)清楚，對(duì)于發(fā)表文章非常不便。發(fā)文章所用的結(jié)果圖最好是等

19、值線圖，并且最好是黑白的等值線圖。一般借用photoshop等第三方軟件，很麻煩，效果也不好。ansys本身也能實(shí)現(xiàn)這項(xiàng)功能。1.將要輸出的結(jié)果調(diào)出，這時(shí)為彩色云圖；2.將云圖轉(zhuǎn)換為等值線圖的形式GUI：plotCtrlsDevice Options/DEVI中的vector mode 選為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：plotCtrlsStyle

20、Contours Contours Labeling 在Key Vector mode contour label 中選中on every Nth elem，然后在N= 輸入框中輸入合適的數(shù)值，例如5，多試幾次，直到疏密合適命令：/clabel,1,55.將彩色等值線變?yōu)楹谏獹UI：plotCtrlsStyleColorsContours Colors 將Items Numbered 1，Items Numbered 2等復(fù)選框中的顏色均選為黑色，圖像即可變?yōu)楹诎椎戎稻€圖像命令：/color,cntr,whit,1 等等6.修改在顯示等值線時(shí)，圖形中太多的標(biāo)識(shí)字母（A、B等）使得等值線圖完全看

21、不清楚。請(qǐng)問有沒有什么方法可以去掉圖形的這些標(biāo)識(shí)字母呢在plotCtrls下面的style下面的contours下面的contour Labeling，在彈出的contour labeling options里面的vector mode contour labels選off后處理中畫出了結(jié)構(gòu)的等值線圖，可否將等值線邊的字母改為相應(yīng)的數(shù)值標(biāo)注/PNUM,SVAL,1/replot如果還想將等值線的字母去掉/PNUM,SVAL,1/CLABEL,1,-1/replot7.出圖GUI：plotCtrlsCapture ImageANSYS中如何加彎矩或扭矩A轉(zhuǎn)矩一般有三種施加的方法:第一種,將矩轉(zhuǎn)換

22、成一對(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),跟其他受力節(jié)點(diǎn)分別形成多根剛性梁,從而形成剛性面.最后也是直接加載荷到中心節(jié)點(diǎn)上面,通過剛性梁來傳遞載荷.上面三種方法計(jì)算的結(jié)果基本一致,我做過實(shí)驗(yàn)的.只不過是后兩種情況都是形成剛性區(qū)域,但是CERIG命令是要在小變形或者小旋轉(zhuǎn)才能用,只支持靜力,線形分析.而第三種方法適用多種情況,不僅支持大應(yīng)變,還支持非線形情況.如果你需要例子,我下

23、次在發(fā)給你看A思路1：矩或扭矩說白了就是矩，所謂矩就是力和力臂的乘積。 施加矩可以等效為施加力；思路2：直接施加彎矩或扭矩，此時(shí)需要引入一個(gè)具有旋轉(zhuǎn)自由度的節(jié)點(diǎn)；可以選擇單元21，或者1841將矩轉(zhuǎn)換成一對(duì)的力偶，直接施加在對(duì)應(yīng)的節(jié)點(diǎn)上面。2在構(gòu)件中心部位建立一個(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)上面，通過剛性梁來傳遞載荷。4通過rbe3命令。該方法與方法2很接近。5基

24、于表面邊界法：主要通過定義一個(gè)接觸表面和一個(gè)目標(biāo)節(jié)點(diǎn)接觸來實(shí)現(xiàn)，彎矩荷載可以通過在目標(biāo)節(jié)點(diǎn)上用“F”命令施加。對(duì)于方法1，通過轉(zhuǎn)換為集中力或均布力，比如施加扭矩，把端面節(jié)點(diǎn)改成柱坐標(biāo)，然后等效為施加環(huán)向的節(jié)點(diǎn)力；而施加彎矩，可以將力矩轉(zhuǎn)化為端面的剪切均布力；但這種方法比較容易出現(xiàn)應(yīng)力集中現(xiàn)象；方法2，定義局部剛性區(qū)域，施加過程venture講的很詳細(xì)，這里就不在贅述。根據(jù)他的例子，我在下面給出了一段命令流。該方法有個(gè)不足，它在端面額外的增加了一定的剛度，只能適用于小變形分析。方法3，相對(duì)方法2來說，采用剛性梁?jiǎn)卧?，適用范圍更廣一些，對(duì)于大應(yīng)變分析也能很好的適用。但在小應(yīng)變分析下，方法2和方法3

25、沒有什么區(qū)別。方法4，定義一個(gè)主節(jié)點(diǎn)，施加了分布力面，應(yīng)該說跟實(shí)際比較接近一點(diǎn)，但端面的結(jié)果好像不是很理想，結(jié)果有點(diǎn)偏大，在遠(yuǎn)離端面處的位置跟實(shí)際很符合。方法5，它具體的受力形式有如下兩種：剛性表面邊界（Rigid surface constraint）認(rèn)為接觸面是剛性的，沒有變形，和通過節(jié)點(diǎn)耦合命令CERIG比較相似； 分布力邊界（Force-distributed constraint）允許接觸面的變形，和邊界定義命令RBE3相似。使用這種方法，需要用KEYOPT(2) = 2打開接觸單元的MPC（多點(diǎn)接觸邊界）算法ANSYS繪制彎矩、剪力、軸力圖1.繪制彎矩圖建立彎矩單元表。例如梁?jiǎn)卧猧

26、節(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,SMISC,12 !單元J點(diǎn)彎矩plls,MI,MJ2.標(biāo)注彎矩圖PLOTCTRLSNUMBERINGSVAL ON即可在畫出彎矩圖的同時(shí)在圖上標(biāo)出彎矩值的大小3.調(diào)整彎矩圖如果彎矩圖方向錯(cuò)誤，則繪制彎矩圖命令為plls,imom,jmom，-1同一個(gè)節(jié)點(diǎn)處兩邊的單元內(nèi)力有細(xì)

27、微差別，導(dǎo)致內(nèi)力數(shù)字標(biāo)注出現(xiàn)重影。觀察上面整體軸力圖也可以發(fā)現(xiàn)，一段一段的，好像馬賽克，其實(shí)上面整體彎矩圖也是，不過不是很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成的（詳細(xì)原因就不展開了，看看梁理論的書和ANSYS的說明吧）。為了修正重影和節(jié)點(diǎn)兩邊內(nèi)力值不一樣的問題，遍制了宏文件ITFAVG.MAC命令文件內(nèi)容如下：!-!宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO）!獲取線性單元內(nèi)力，并對(duì)單元邊界處的內(nèi)力進(jìn)行平衡!輸入信息!內(nèi)力類型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ*ASK,ITFTY

28、PE,PLEASE INPUT THE TYPE OF INTERNAL FORCE,MMOMY!需處理的單元包*ASK,EASSEMBLY,PLEASE INPUT THE COMPONENT NAME OF ELEMENTS TO BE PROCESSED!,EOUTER!需處理的節(jié)點(diǎn)包*ASK,NASSEMBLY,PLEASE INPUT THE COMPONENT NAME OF NODE TO BE PROCESSED!,NOUTER!無需處理的節(jié)點(diǎn)包*ASK,UNASSEMBLY,PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NO

29、DE!(NONE IF THERES NO SUCH COMPONENT),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,ITFTYPE,EQ,MFORY,THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,MFORZ,THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYP

30、E,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ù)目CMSEL,S,UNASSEMBLY*GET,UNNODNUM,NODE,0,COUNT!定義長(zhǎng)度為UNNODNUM的數(shù)組（UNNOD），以存放選中單元的單元編號(hào)*DIM,UNNOD,AR

31、RAY,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)度為SELELENUM的數(shù)組（ELENUM），以存放選中單元的單元編號(hào)*DIM,ELENUM,ARRAY,SELELENUM!將選中單元的編號(hào)按順序存入數(shù)組ELENUM*DO,I,0,SELEL

32、ENUM-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+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,ITFINU

33、METABLE,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!選出和節(jié)點(diǎn)K相連的

34、線性單元中，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=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的單

35、元的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,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO

36、!END OF ITFAVG.MAC（2）對(duì)體和面來說，ANSYS默認(rèn)的結(jié)果輸出格式是云圖格式，而這種彩色云圖打印為黑白圖像時(shí)對(duì)比很不明顯，無法表達(dá)清楚，對(duì)于發(fā)表文章非常不便。發(fā)文章所用的結(jié)果圖最好是等值線圖，并且最好是黑白的等值線圖。一般借用photoshop等第三方軟件，很麻煩，效果也不好。ansys本身也能實(shí)現(xiàn)這項(xiàng)功能。1.將要輸出的結(jié)果調(diào)出，這時(shí)為彩色云圖；2.將云圖轉(zhuǎn)換為等值線圖的形式GUI：plotCtrlsDevice Options/DEVI中的vector mode 選為on命令：/DEVICE,VECTOR,1這時(shí)結(jié)果為彩色等值線，若直接輸出，打印為黑白圖像時(shí)仍然不清晰，為

37、此需進(jìn)行以下幾步將圖像轉(zhuǎn)換為黑白形式；3. 將背景變?yōu)榘咨睿簀pgprf,500,100,1/rep4.對(duì)等值線中的等值線符號(hào)（圖中為A,B,C等）的疏密進(jìn)行調(diào)整GUI：plotCtrlsStyleContours Contours Labeling 在Key Vector mode contour label 中選中on every Nth elem，然后在N= 輸入框中輸入合適的數(shù)值，例如5，多試幾次，直到疏密合適命令：/clabel,1,55.將彩色等值線變?yōu)楹谏獹UI：plotCtrlsStyleColorsContours Colors 將Items Numbered 1，Ite

38、ms Numbered 2等復(fù)選框中的顏色均選為黑色，圖像即可變?yōu)楹诎椎戎稻€圖像命令：/color,cntr,whit,1 等等6.修改在顯示等值線時(shí)，圖形中太多的標(biāo)識(shí)字母（A、B等）使得等值線圖完全看不清楚。請(qǐng)問有沒有什么方法可以去掉圖形的這些標(biāo)識(shí)字母呢在plotCtrls下面的style下面的contours下面的contour Labeling，在彈出的contour labeling options里面的vector mode contour labels選off后處理中畫出了結(jié)構(gòu)的等值線圖，可否將等值線邊的字母改為相應(yīng)的數(shù)值標(biāo)注/PNUM,SVAL,1/replot如果還想將等值線的

39、字母去掉/PNUM,SVAL,1/CLABEL,1,-1/replot7.出圖GUI：plotCtrlsCapture Image=【原創(chuàng)】下面部分為作者自己總結(jié)對(duì)于beam3單元其他的單元類型，可以參考ansys幫助文檔中給出的相關(guān)單元的參數(shù)。附上一個(gè)小例子，下圖所示為一懸臂梁，長(zhǎng)100，左端完全約束，上面加載均勻載荷q=100，繪制這個(gè)梁的剪力和彎矩圖。彎矩圖剪力圖1.繪制彎矩圖 建立彎矩單元表。例如梁?jiǎn)卧?i節(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,SMI

40、SC,2 !單元I點(diǎn)剪力ETABLE,QJ,SMISC,8 !單元J點(diǎn)剪力ETABLE,MI,SMISC,6 !單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12 !單元J點(diǎn)彎矩plls,MI,MJ 2.標(biāo)注彎矩圖 PLOTCTRLSNUMBERINGSVAL ON即可在畫出彎矩圖的同時(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í)上面整體彎矩圖也是，不過不是 很明顯罷了。這是EULER-BEONOUL

41、I梁理論以及ANSYS輸出定義造成 的（詳細(xì)原因就不展開了，看看梁理論的書和ANSYS的說明吧）。 為了修正重影和節(jié)點(diǎn)兩邊內(nèi)力值不一樣的問題，遍制了宏文件ITFAVG.MAC 命令文件內(nèi)容如下： !- !宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO） !獲取線性單元內(nèi)力，并對(duì)單元邊界處的內(nèi)力進(jìn)行平衡 !輸入信息 !內(nèi)力類型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ *ASK,ITFTYPE,PLEASE INPUT THE TYPE OF INTERNAL FORCE,MMOMY !需處理的單元包 *ASK,EASSEMBLY

42、,PLEASE INPUT THE COMPONENT NAME OF ELEMENTS TO BE PROCESSED!, EOUTER !需處理的節(jié)點(diǎn)包 *ASK,NASSEMBLY,PLEASE INPUT THE COMPONENT NAME OF NODE TO BE PROCESSED!,NOU TER !無需處理的節(jié)點(diǎn)包 *ASK,UNASSEMBLY,PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NODE!(NONE I F THERES NO SUCH COMPONENT),NONE /POST1 !輸入信息：內(nèi)力類型，欲處

43、理單元的集合，欲處理節(jié)點(diǎn)的集合 !ITFTYPE=MMOMY !EASSEMBLY=EOUTER !NASSEMBLY=NOUTER !按內(nèi)力類型確定ANSYS輸出信息SMISC的編號(hào) *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,IT

44、FTYPE,EQ,MMOMY,THEN ITFINUM=5 ITFJNUM=11 *ELSEIF,ITFTYPE,EQ,MMOMZ,THEN ITFINUM=6 ITFJNUM=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ù)目 CMSEL,S,UNASSEMBLY *GET,UNNODNUM,NODE,0,COUNT !定義長(zhǎng)度為UNNODNUM的數(shù)組（UNNOD），以存放選中單元的單元編號(hào) *DIM,UNNOD,ARRAY,UNNODNUM !將選中單元的編號(hào)按順序存入數(shù)組U

45、NNOD *DO,I,0,UNNODNUM-1,1 UNNOD(I+1)=NDNEXT(I) *ENDDO *ELSE UNNODNUM=0 *ENDIF !選出所需的單元和節(jié)點(diǎn)包 CMSEL,S,EASSEMBLY CMSEL,S,NASSEMBLY !獲得當(dāng)前選中單元總數(shù)（存入變量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,1 ELENU

46、M(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,1 NODNUM(I+1)=NDNEXT(I) *ENDDO !定義所需的線性單元內(nèi)力ETABLE，節(jié)點(diǎn)I的內(nèi)力存入數(shù)組ITNFI, !節(jié)點(diǎn)J的內(nèi)力存入數(shù)組ITNFJ ETABLE,ITNFI,SMISC,ITFINUM

47、 ETABLE,ITNFJ,SMISC,ITFJNUM !定義所需的結(jié)果數(shù)組，并將其置零 ETABLE,ITNFINEO,SMISC,5 SADD,ITNFINEO,ITNFI,1 ETABLE,ITNFJNEO,SMISC,11 SADD,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),THEN INDEX=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),E