ANSYS中彎矩、剪力圖的繪制

-.z.ansys中如何生成命令流方法：GUI是：UtilityMenu>File>WriteDBLogFile怎么用ansys繪制彎矩，剪力圖：GUI:GeneralPostproc->lotResult->ContourPlot->LineElementResult

彈出畫單元結(jié)果的對話框，分別在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,SMISC,6!單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12!單元J點(diǎn)彎矩!更新單元表ETABLE,REFL

!畫軸力分布圖/TITLE,A*ialforcediagram

PLLS,NI,NJ,1.0,0/image,save,'A*ial_force_%T%',jpg

!畫剪力分布圖/TITLE,Shearingforcediagram

PLLS,QI,QJ,1.0,0/image,save,'Shearing_force_%T%',jpg

!畫彎矩分布圖/TITLE,Bendingmomentdiagram

PLLS,MI,MJ,-0.8,0

/image,save,'Bending_moment_%T%',jpgANSYS中彎矩、剪力圖的繪制GUI:

GeneralPostproc-plotResult-ContourPlot-LineElementResult

彈出畫單元結(jié)果的對話框，分別在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,SMISC,6!單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12!單元J點(diǎn)彎矩!更新單元表ETABLE,REFL

!畫軸力分布圖/TITLE,A*ialforcediagram

PLLS,NI,NJ,1.0,0/image,save,'A*ial_force_%T%',jpg

!畫剪力分布圖/TITLE,Shearingforcediagram

PLLS,QI,QJ,1.0,0/image,save,'Shearing_force_%T%',jpg

!畫彎矩分布圖/TITLE,Bendingmomentdiagram

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,2600mp,e*,1,3.06e10mp,pr*y,1,0.2sectype,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,u*,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,-1ansys如何繪制彎矩圖Ansy中彎矩圖，云圖繪制總結(jié)在回答別人問題時(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.繪制彎矩圖建立彎矩單元表。例如梁單元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,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,imom,jmom2.標(biāo)注彎矩圖PLOTCTRLS>>NUMBERING>>SVALON即可在畫出彎矩圖的同時(shí)在圖上標(biāo)出彎矩值的大小3.調(diào)整彎矩圖如果彎矩圖方向錯(cuò)誤，則繪制彎矩圖命令為plls,imom,jmom，-1同一個(gè)節(jié)點(diǎn)處兩邊的單元力有細(xì)微差別，導(dǎo)致力數(shù)字標(biāo)注出現(xiàn)重影。觀察上面整體軸力圖也可以發(fā)現(xiàn)，一段一段的，好像馬賽克，其實(shí)上面整體彎矩圖也是，不過不是很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成的（詳細(xì)原因就不展開了，看看梁理論的書和ANSYS的說明吧）。為了修正重影和節(jié)點(diǎn)兩邊力值不一樣的問題，遍制了宏文件ITFAVG.MAC命令文件容如下：!---------------------------------------------------------------------!宏：ITFAVG.MAC（INTERNALFORCEAVERAGEMACRO）!獲取線性單元力，并對單元邊界處的力進(jìn)行平衡!輸入信息!力類型：MFOR*,MFORY,MFORZ,MMOM*,MMOMY,MMOMZ*ASK,ITFTYPE,'PLEASEINPUTTHETYPEOFINTERNALFORCE','MMOMY'!需處理的單元包*ASK,EASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFELEMENTSTOBEPROCESSED!','EOUTER'!需處理的節(jié)點(diǎn)包*ASK,NASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFNODETOBEPROCESSED!','NOUTER'!無需處理的節(jié)點(diǎn)包*ASK,UNASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFTHEUNCHANGEDNODE!(NONEIFTHERE'SNOSUCHPONENT)','NONE'/POST1!輸入信息：力類型，欲處理單元的集合，欲處理節(jié)點(diǎn)的集合!ITFTYPE='MMOMY'!EASSEMBLY='EOUTER'!NASSEMBLY='NOUTER'!按力類型確定ANSYS輸出信息SMISC的編號*IF,ITFTYPE,EQ,'MFOR*',THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOM*',THENITFINUM=4ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6ITFJNUM=12*ELSE*ENDIF!對不需平均的節(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!定義長度為UNNODNUM的數(shù)組（UNNOD），以存放選中單元的單元編號*DIM,UNNOD,ARRAY,UNNODNUM!將選中單元的編號按順序存入數(shù)組UNNOD*DO,I,0,UNNODNUM-1,1UNNOD(I+1)=NDNE*T(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!選出所需的單元和節(jié)點(diǎn)包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!獲得當(dāng)前選中單元總數(shù)（存入變量SELELENUM）*GET,SELELENUM,ELEM,0,COUNT!定義長度為SELELENUM的數(shù)組（ELENUM），以存放選中單元的單元編號*DIM,ELENUM,ARRAY,SELELENUM!將選中單元的編號按順序存入數(shù)組ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNE*T(I)*ENDDO!獲得當(dāng)前選中節(jié)點(diǎn)總數(shù)（存入變量SELNODNUM）*GET,SELNODNUM,NODE,0,COUNT!定義長度為SELNODNUM的數(shù)組（NODNUM），以存放選中單元的單元編號*DIM,NODNUM,ARRAY,SELNODNUM!將選中單元的編號按順序存入數(shù)組NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNE*T(I)*ENDDO!定義所需的線性單元力ETABLE，節(jié)點(diǎn)I的力存入數(shù)組ITNFI,!節(jié)點(diǎn)J的力存入數(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)INDE*=0*IF,UNNODNUM,GE,1,THEN*DO,J,1,UNNODNUM*IF,NODNUM(K),EQ,UNNOD(J),THENINDE*=1*ELSE*ENDIF*ENDDO*ELSE*ENDIF*DO,J,1,SELELENUM,1!選出和節(jié)點(diǎn)K相連的線性單元中，I節(jié)點(diǎn)（對線性單元而言）為節(jié)點(diǎn)K的單元編號*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*E*IT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!選出和節(jié)點(diǎn)K相連的線性單元中，J節(jié)點(diǎn)（對線性單元而言）為節(jié)點(diǎn)K的單元編號*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THENELEJ=ELENUM(J)*E*IT*ELSE*ENDIF*ENDDO*IF,INDE*,EQ,0,THEN*IF,ELEJ,NE,0,THEN!有可能出現(xiàn)ELEJ為0的情況!取出I節(jié)點(diǎn)為節(jié)點(diǎn)K的單元的I節(jié)點(diǎn)端的力放入?yún)?shù)ETELEI*GET,ETELEI,ELEM,ELEI,ETAB,ITNFI!取出J節(jié)點(diǎn)為節(jié)點(diǎn)K的單元的J節(jié)點(diǎn)端的力放入?yún)?shù)ETELEJ*GET,ETELEJ,ELEM,ELEJ,ETAB,ITNFJ!平均節(jié)點(diǎn)K的單元的I節(jié)點(diǎn)端的力和節(jié)點(diǎn)K的單元的J節(jié)點(diǎn)端的力ETAVE=(ETELEI+ETELEJ)/2!將平均后的力存入結(jié)果數(shù)組中DETAB,ELEI,ITNFINEO,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!ENDOFITFAVG.MAC（2）對體和面來說，ANSYS默認(rèn)的結(jié)果輸出格式是云圖格式，而這種彩色云圖打印為黑白圖像時(shí)對比很不明顯，無法表達(dá)清楚，對于發(fā)表文章非常不便。發(fā)文章所用的結(jié)果圖最好是等值線圖，并且最好是黑白的等值線圖。一般借用photoshop等第三方軟件，很麻煩，效果也不好。ansys本身也能實(shí)現(xiàn)這項(xiàng)功能。1.將要輸出的結(jié)果調(diào)出，這時(shí)為彩色云圖；2.將云圖轉(zhuǎn)換為等值線圖的形式GUI：plotCtrls—>DeviceOptions—>[/DEVI]中的vectormode選為on命令：/DEVICE,VECTOR,1這時(shí)結(jié)果為彩色等值線，若直接輸出，打印為黑白圖像時(shí)仍然不清晰，為此需進(jìn)行以下幾步將圖像轉(zhuǎn)換為黑白形式；3.將背景變?yōu)榘咨睿簀pgprf,500,100,1/rep4.對等值線中的等值線符號（圖中為A,B,C等）的疏密進(jìn)行調(diào)整GUI：plotCtrls—>Style—>Contours—>ContoursLabeling在KeyVectormodecontourlabel中選中oneveryNthelem，然后在N=輸入框中輸入合適的數(shù)值，例如5，多試幾次，直到疏密合適命令：/clabel,1,55.將彩色等值線變?yōu)楹谏獹UI：plotCtrls—>Style—>Colors—>ContoursColors將ItemsNumbered1，ItemsNumbered2等復(fù)選框中的顏色均選為黑色，圖像即可變?yōu)楹诎椎戎稻€圖像命令：/color,tr,whit,1等等6.修改在顯示等值線時(shí)，圖形中太多的標(biāo)識字母（A、B等）使得等值線圖完全看不清楚。請問有沒有什么方法可以去掉圖形的這些標(biāo)識字母呢在plotCtrls下面的style下面的contours下面的contourLabeling，在彈出的contourlabelingoptions里面的vectormodecontourlabels選off后處理中畫出了結(jié)構(gòu)的等值線圖，可否將等值線邊的字母改為相應(yīng)的數(shù)值標(biāo)注/PNUM,SVAL,1/replot如果還想將等值線的字母去掉/PNUM,SVAL,1/CLABEL,1,-1/replot7.出圖GUI：plotCtrls—>CaptureImageANSYS中如何加彎矩或扭矩A轉(zhuǎn)矩一般有三種施加的方法:第一種,將矩轉(zhuǎn)換成一對一對的力偶,直接施加在對應(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)變,還支持非線形情況.

如果你需要例子,我下次在發(fā)給你看A思路1：矩或扭矩說白了就是矩，所謂矩就是力和力臂的乘積。施加矩可以等效為施加力；思路2：直接施加彎矩或扭矩，此時(shí)需要引入一個(gè)具有旋轉(zhuǎn)自由度的節(jié)點(diǎn)；可以選擇單元21，或者1841．將矩轉(zhuǎn)換成一對的力偶，直接施加在對應(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．基于表面邊界法：主要通過定義一個(gè)接觸表面和一個(gè)目標(biāo)節(jié)點(diǎn)接觸來實(shí)現(xiàn)，彎矩荷載可以通過在目標(biāo)節(jié)點(diǎn)上用“F”命令施加。對于方法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，相對方法2來說，采用剛性梁單元，適用圍更廣一些，對于大應(yīng)變分析也能很好的適用。但在小應(yīng)變分析下，方法2和方法3沒有什么區(qū)別。方法4，定義一個(gè)主節(jié)點(diǎn)，施加了分布力面，應(yīng)該說跟實(shí)際比較接近一點(diǎn)，但端面的結(jié)果好像不是很理想，結(jié)果有點(diǎn)偏大，在遠(yuǎn)離端面處的位置跟實(shí)際很符合。方法5，它具體的受力形式有如下兩種：剛性表面邊界（Rigidsurfaceconstraint）－認(rèn)為接觸面是剛性的，沒有變形，和通過節(jié)點(diǎn)耦合命令CERIG比較相似；分布力邊界（Force-distributedconstraint）－允許接觸面的變形，和邊界定義命令RBE3相似。使用這種方法，需要用KEYOPT(2)=2打開接觸單元的MPC（多點(diǎn)接觸邊界）算法ANSYS繪制彎矩、剪力、軸力圖1.繪制彎矩圖建立彎矩單元表。例如梁單元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,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)注彎矩圖PLOTCTRLS>>NUMBERING>>SVALON即可在畫出彎矩圖的同時(shí)在圖上標(biāo)出彎矩值的大小3.調(diào)整彎矩圖如果彎矩圖方向錯(cuò)誤，則繪制彎矩圖命令為plls,imom,jmom，-1同一個(gè)節(jié)點(diǎn)處兩邊的單元力有細(xì)微差別，導(dǎo)致力數(shù)字標(biāo)注出現(xiàn)重影。觀察上面整體軸力圖也可以發(fā)現(xiàn)，一段一段的，好像馬賽克，其實(shí)上面整體彎矩圖也是，不過不是很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成的（詳細(xì)原因就不展開了，看看梁理論的書和ANSYS的說明吧）。為了修正重影和節(jié)點(diǎn)兩邊力值不一樣的問題，遍制了宏文件ITFAVG.MAC命令文件容如下：!---------------------------------------------------------------------!宏：ITFAVG.MAC（INTERNALFORCEAVERAGEMACRO）!獲取線性單元力，并對單元邊界處的力進(jìn)行平衡!輸入信息!力類型：MFOR*,MFORY,MFORZ,MMOM*,MMOMY,MMOMZ*ASK,ITFTYPE,'PLEASEINPUTTHETYPEOFINTERNALFORCE','MMOMY'!需處理的單元包*ASK,EASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFELEMENTSTOBEPROCESSED!','EOUTER'!需處理的節(jié)點(diǎn)包*ASK,NASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFNODETOBEPROCESSED!','NOUTER'!無需處理的節(jié)點(diǎn)包*ASK,UNASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFTHEUNCHANGEDNODE!(NONEIFTHERE'SNOSUCHPONENT)','NONE'/POST1!輸入信息：力類型，欲處理單元的集合，欲處理節(jié)點(diǎn)的集合!ITFTYPE='MMOMY'!EASSEMBLY='EOUTER'!NASSEMBLY='NOUTER'!按力類型確定ANSYS輸出信息SMISC的編號*IF,ITFTYPE,EQ,'MFOR*',THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOM*',THENITFINUM=4ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6ITFJNUM=12*ELSE*ENDIF!對不需平均的節(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!定義長度為UNNODNUM的數(shù)組（UNNOD），以存放選中單元的單元編號*DIM,UNNOD,ARRAY,UNNODNUM!將選中單元的編號按順序存入數(shù)組UNNOD*DO,I,0,UNNODNUM-1,1UNNOD(I+1)=NDNE*T(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!選出所需的單元和節(jié)點(diǎn)包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!獲得當(dāng)前選中單元總數(shù)（存入變量SELELENUM）*GET,SELELENUM,ELEM,0,COUNT!定義長度為SELELENUM的數(shù)組（ELENUM），以存放選中單元的單元編號*DIM,ELENUM,ARRAY,SELELENUM!將選中單元的編號按順序存入數(shù)組ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNE*T(I)*ENDDO!獲得當(dāng)前選中節(jié)點(diǎn)總數(shù)（存入變量SELNODNUM）*GET,SELNODNUM,NODE,0,COUNT!定義長度為SELNODNUM的數(shù)組（NODNUM），以存放選中單元的單元編號*DIM,NODNUM,ARRAY,SELNODNUM!將選中單元的編號按順序存入數(shù)組NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNE*T(I)*ENDDO!定義所需的線性單元力ETABLE，節(jié)點(diǎn)I的力存入數(shù)組ITNFI,!節(jié)點(diǎn)J的力存入數(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)INDE*=0*IF,UNNODNUM,GE,1,THEN*DO,J,1,UNNODNUM*IF,NODNUM(K),EQ,UNNOD(J),THENINDE*=1*ELSE*ENDIF*ENDDO*ELSE*ENDIF*DO,J,1,SELELENUM,1!選出和節(jié)點(diǎn)K相連的線性單元中，I節(jié)點(diǎn)（對線性單元而言）為節(jié)點(diǎn)K的單元編號*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*E*IT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!選出和節(jié)點(diǎn)K相連的線性單元中，J節(jié)點(diǎn)（對線性單元而言）為節(jié)點(diǎn)K的單元編號*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THENELEJ=ELENUM(J)*E*IT*ELSE*ENDIF*ENDDO*IF,INDE*,EQ,0,THEN*IF,ELEJ,NE,0,THEN!有可能出現(xiàn)ELEJ為0的情況!取出I節(jié)點(diǎn)為節(jié)點(diǎn)K的單元的I節(jié)點(diǎn)端的力放入?yún)?shù)ETELEI*GET,ETELEI,ELEM,ELEI,ETAB,ITNFI!取出J節(jié)點(diǎn)為節(jié)點(diǎn)K的單元的J節(jié)點(diǎn)端的力放入?yún)?shù)ETELEJ*GET,ETELEJ,ELEM,ELEJ,ETAB,ITNFJ!平均節(jié)點(diǎn)K的單元的I節(jié)點(diǎn)端的力和節(jié)點(diǎn)K的單元的J節(jié)點(diǎn)端的力ETAVE=(ETELEI+ETELEJ)/2!將平均后的力存入結(jié)果數(shù)組中DETAB,ELEI,ITNFINEO,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!ENDOFITFAVG.MAC（2）對體和面來說，ANSYS默認(rèn)的結(jié)果輸出格式是云圖格式，而這種彩色云圖打印為黑白圖像時(shí)對比很不明顯，無法表達(dá)清楚，對于發(fā)表文章非常不便。發(fā)文章所用的結(jié)果圖最好是等值線圖，并且最好是黑白的等值線圖。一般借用photoshop等第三方軟件，很麻煩，效果也不好。ansys本身也能實(shí)現(xiàn)這項(xiàng)功能。1.將要輸出的結(jié)果調(diào)出，這時(shí)為彩色云圖；2.將云圖轉(zhuǎn)換為等值線圖的形式GUI：plotCtrls—>DeviceOptions—>[/DEVI]中的vectormode選為on命令：/DEVICE,VECTOR,1這時(shí)結(jié)果為彩色等值線，若直接輸出，打印為黑白圖像時(shí)仍然不清晰，為此需進(jìn)行以下幾步將圖像轉(zhuǎn)換為黑白形式；3.將背景變?yōu)榘咨睿簀pgprf,500,100,1/rep4.對等值線中的等值線符號（圖中為A,B,C等）的疏密進(jìn)行調(diào)整GUI：plotCtrls—>Style—>Contours—>ContoursLabeling在KeyVectormodecontourlabel中選中oneveryNthelem，然后在N=輸入框中輸入合適的數(shù)值，例如5，多試幾次，直到疏密合適命令：/clabel,1,55.將彩色等值線變?yōu)楹谏獹UI：plotCtrls—>Style—>Colors—>ContoursColors將ItemsNumbered1，ItemsNumbered2等復(fù)選框中的顏色均選為黑色，圖像即可變?yōu)楹诎椎戎稻€圖像命令：/color,tr,whit,1等等6.修改在顯示等值線時(shí)，圖形中太多的標(biāo)識字母（A、B等）使得等值線圖完全看不清楚。請問有沒有什么方法可以去掉圖形的這些標(biāo)識字母呢在plotCtrls下面的style下面的contours下面的contourLabeling，在彈出的contourlabelingoptions里面的vectormodecontourlabels選off后處理中畫出了結(jié)構(gòu)的等值線圖，可否將等值線邊的字母改為相應(yīng)的數(shù)值標(biāo)注/PNUM,SVAL,1/replot如果還想將等值線的字母去掉/PNUM,SVAL,1/CLABEL,1,-1/replot7.出圖GUI：plotCtrls—>CaptureImage============================================================================【原創(chuàng)】下面部分為作者自己總結(jié)對于beam3單元其他的單元類型，可以參考ansys幫助文檔中給出的相關(guān)單元的參數(shù)。附上一個(gè)小例子，下圖所示為一懸臂梁，長100，左端完全約束，上面加載均勻載荷q=100，繪制這個(gè)梁的剪力和彎矩圖。彎矩圖剪力圖1.繪制彎矩圖建立彎矩單元表。例如梁單元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,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)注彎矩圖PLOTCTRLS>>NUMBERING>>SVALON即可在畫出彎矩圖的同時(shí)在圖上標(biāo)出彎矩值的大小3.調(diào)整彎矩圖如果彎矩圖方向錯(cuò)誤，則繪制彎矩圖命令為plls,imom,jmom，-1同一個(gè)節(jié)點(diǎn)處兩邊的單元力有細(xì)微差別，導(dǎo)致力數(shù)字標(biāo)注出現(xiàn)重影。觀察上面整體軸力圖也可以發(fā)現(xiàn)，一段一段的，好像馬賽克，其實(shí)上面整體彎矩圖也是，不過不是很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成的（詳細(xì)原因就不展開了，看看梁理論的書和ANSYS的說明吧）。為了修正重影和節(jié)點(diǎn)兩邊力值不一樣的問題，遍制了宏文件ITFAVG.MAC命令文件容如下：!---------------------------------------------------------------------!宏：ITFAVG.MAC（INTERNALFORCEAVERAGEMACRO）!獲取線性單元力，并對單元邊界處的力進(jìn)行平衡!輸入信息!力類型：MFOR*,MFORY,MFORZ,MMOM*,MMOMY,MMOMZ*ASK,ITFTYPE,'PLEASEINPUTTHETYPEOFINTERNALFORCE','MMOMY'!需處理的單元包*ASK,EASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFELEMENTSTOBEPROCESSED!','EOUTER'!需處理的節(jié)點(diǎn)包*ASK,NASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFNODETOBEPROCESSED!','NOUTER'!無需處理的節(jié)點(diǎn)包*ASK,UNASSEMBLY,'PLEASEINPUTTHEPONENTNAMEOFTHEUNCHANGEDNODE!(NONEIFTHERE'SNOSUCHPONENT)','NONE'/POST1!輸入信息：力類型，欲處理單元的集合，欲處理節(jié)點(diǎn)的集合!ITFTYPE='MMOMY'!EASSEMBLY='EOUTER'!NASSEMBLY='NOUTER'!按力類型確定ANSYS輸出信息SMISC的編號*IF,ITFTYPE,EQ,'MFOR*',THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOM*',THENITFINUM=4ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6ITFJNUM=12*ELSE*ENDIF!對不需平均的節(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!定義長度為UNNODNUM的數(shù)組（UNNOD），以存放選中單元的單元編號*DIM,UNNOD,ARRAY,UNNODNUM!將選中單元的編號按順序存入數(shù)組UNNOD*DO,I,0,UNNODNUM-1,1UNNOD(I+1)=NDNE*T(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!選出所需的單元和節(jié)點(diǎn)包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!獲得當(dāng)前選中單元總數(shù)（存入變量SELELENUM）*GET,SELELENUM,ELEM,0,COUNT!定義長度為SELELENUM的數(shù)組（ELENUM），以存放選中單元的單元編號*DIM,ELENUM,ARRAY,SELELENUM!將選中單元的編號按順序存入數(shù)組ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNE*T(I)*ENDDO!獲得當(dāng)前選中節(jié)點(diǎn)總數(shù)（存入變量SELNODNUM）*GET,SELNODNUM,NODE,0,COUNT!定義長度為SELNODNUM的數(shù)組（NODNUM），以存放選中單元的單元編號*DIM,NODNUM,ARRAY,SELNODNUM!將選中單元的編號按順序存入數(shù)組NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNE*T(I)*ENDDO!定義所需的線性單元力ETABLE，節(jié)點(diǎn)I的力存入數(shù)組ITNFI,!節(jié)點(diǎn)J的力存入數(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)INDE*=0*IF,UNNODNUM,GE,1,THEN*DO,J,1,UNNODNUM*IF,NODNUM(K),EQ,UNNOD(J),THENINDE*=1*ELSE*ENDIF*ENDDO*ELSE*ENDIF*DO,J,1,SELELENUM,1!選出