ANSYS中彎矩、剪力圖的繪制

1、Ansys應(yīng)用技巧ansys中如何生成命令流方法：GUI是：Utility Menu>File>Write DB Log FileANSYS中彎矩、剪力圖的繪制GUI:General Postproc-plot Result-Contour Plot-Line Element Result彈出畫(huà)單元結(jié)果的對(duì)話框，分別在Labi和Labj依次選取SMIS6和SMIS12（彎矩圖）、SMIS1和SMIS7（軸力圖）、SMIS2和SMIS8（剪力圖）! 建立單元表ETABLE,NI,SMISC,1     !單元I點(diǎn)軸力ETABLE,NJ,SMISC,7 &

2、#160;   !單元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     ! 畫(huà)軸力分布圖/TITLE,Axial   force  

3、; diagramPLLS,NI,NJ,1.0,0     /image,save,'Axial_force_%T%',jpg     ! 畫(huà)剪力分布圖/TITLE,Shearing   force   diagramPLLS,QI,QJ,1.0,0 /image,save,'Shearing_force_%T%',jpg     ! 畫(huà)彎矩分布圖/TITLE,Bending 

4、0; 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

5、,0.25,0.35,asel,u,1cm,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,2600 mp,ex,1,3.06e10           mp,prxy,1,0.2  

6、;    sectype,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,sh

7、eari,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)注在圖上嗎？如何實(shí)現(xiàn)？ 如果三維問(wèn)題，在剖面上標(biāo)出某一結(jié)構(gòu)的軸心力、彎矩等，如何實(shí)現(xiàn) （2）后處理圖形，其等值線的數(shù)值能否直接標(biāo)注在圖上，而不是采

8、用圖例的形式 后處理結(jié)果往往用云圖表示，下跟一圖例表示數(shù)值大小，能夠?qū)崿F(xiàn)等值線直接標(biāo)注在圖上 回答 （1） 1.繪制彎矩圖 建立彎矩單元表。例如梁?jiǎn)卧?#160;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,SMI

9、SC,8     !單元J點(diǎn)剪力ETABLE,MI,SMISC,6     !單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12 !單元J點(diǎn)彎矩plls,imom,jmom 2.標(biāo)注彎矩圖 PLOTCTRLS>>NUMBERING>>SVAL ON即可在畫(huà)出彎矩圖的同時(shí)在圖上標(biāo)出彎矩值的大小 3.調(diào)整彎矩圖 如果彎矩圖方向錯(cuò)誤，則繪制彎矩圖命令為 plls,imom,jmom，-1 同一個(gè)節(jié)點(diǎn)處兩邊的單元內(nèi)力有細(xì)微差別， 導(dǎo)

10、致內(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.MAC 命令檔內(nèi)容如下：!-!宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO） !獲取線性單元內(nèi)力，并對(duì)單元邊界處的內(nèi)力進(jìn)行平衡!輸入信息 !內(nèi)力類型：MFORX,MFORY,MFO

11、RZ,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'!需處理的節(jié)點(diǎn)包 *ASK,NASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF NODE

12、TO BE PROCESSED!','NOU TER'!無(wú)需處理的節(jié)點(diǎn)包 *ASK,UNASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NODE!(NONE I F THERE'S NO SUCH COMPONENT)','NONE'/POST1!輸入信息：內(nèi)力類型，欲處理單元的集合，欲處理節(jié)點(diǎn)的集合 !ITFTYPE='MMOMY' !EASSEMBLY='EOUTER'

13、0;!NASSEMBLY='NOUTER'!按內(nèi)力類型確定ANSYS輸出信息SMISC的編號(hào) *IF,ITFTYPE,EQ,'MFORX',THENITFINUM=1 ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2 ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3 ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOMX',THENITFINUM=4 I

14、TFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5 ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=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的

15、數(shù)組（UNNOD），以存放選中單元的單元編號(hào) *DIM,UNNOD,ARRAY,UNNODNUM!將選中單元的編號(hào)按順序存入數(shù)組UNNOD *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ù)（存入變數(shù)SELELENUM） *GET,SELELENUM,ELEM,0,COUNT!定義長(zhǎng)度為

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

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

18、0;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),EQ,NO

19、DNUM(K),THEN ELEI=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),THEN ELEJ=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的單元

20、的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,ETAVE DETAB,ELEJ,ITNFJNEO,ETAVE*ELSE *ENDIF*ELSE *EN

21、DIF*ENDDO/UDOC,1,LOGO,OFF PLLS,ITNFINEO,ITNFJNEO!END OF ITFAVG.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：plotCtrls>D

22、evice Options>/DEVI中的vector mode 選為on 命令：/DEVICE,VECTOR,1 這時(shí)結(jié)果為彩色等值線，若直接輸出，打印為黑白圖像時(shí)仍然不清晰，為此需進(jìn)行以下幾步將圖像轉(zhuǎn)換為黑白形式； 3. 將背景變?yōu)榘咨?#160;命令：jpgprf,500,100,1 /rep 4.對(duì)等值線中的等值線符號(hào)（圖中為A,B,C等）的疏密進(jìn)行調(diào)整 GUI：plotCtrls>Style>Contours> Contours Labeling 在Key Vector mode contour la

23、bel 中選中on every Nth elem，然后在N= 輸入框中輸入合適的數(shù)值，例如5，多試幾次，直到疏密合適 命令：/clabel,1,5 5.將彩色等值線變?yōu)楹谏?#160;GUI：plotCtrls>Style>Colors>Contours Colors 將Items Numbered 1，Items Numbered 2等復(fù)選框中的顏色均選為黑色，圖像即可變?yōu)楹诎椎戎稻€圖像 命令：/color,cntr,whit,1 等等 6.修改 在顯示等值線時(shí)，圖形中太多的標(biāo)識(shí)字母（A、B等）使得等值線圖完全看不清楚。請(qǐng)問(wèn)

24、有沒(méi)有什么方法可以去掉圖形的這些標(biāo)識(shí)字母呢 在plotCtrls下面的style下面的contours下面的contour Labeling，在彈出的contour labeling options里面的vector mode contour labels選off 后處理中畫(huà)出了結(jié)構(gòu)的等值線圖，可否將等值線邊的字母改為相應(yīng)的數(shù)值標(biāo)注 /PNUM,SVAL,1 /replot 如果還想將等值線的字母去掉 /PNUM,SVAL,1 /CLABEL,1,-1 /replot 7.出圖 GUI：plot

25、Ctrls>Capture ImageANSYS中如何加彎矩或扭矩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),跟其它受力節(jié)點(diǎn)分別形成多根剛性梁,從而形成剛性面.最后也是直接加載荷到中心節(jié)點(diǎn)上面,通過(guò)剛性梁來(lái)傳遞載荷.上面三種方法計(jì)算的結(jié)果基本一致,我做過(guò)實(shí)驗(yàn)的.只不過(guò)是后兩種情況都是形成剛性區(qū)域,但是CERIG命令是要

26、在小變形或者小旋轉(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)件中心部位建立一個(gè)節(jié)點(diǎn)，定義為mass21單元，然后跟其它受力節(jié)點(diǎn)耦合，形成剛性區(qū)域，就是用cerig命令。然后直接加轉(zhuǎn)矩到主節(jié)點(diǎn)，即中心節(jié)點(diǎn)上面。3

27、使用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)化為端面的剪切均布力；但這種方法比較容易出現(xiàn)應(yīng)力集中現(xiàn)象；方法2，定義局部剛性區(qū)域，施加過(guò)程venture講的很詳細(xì)，這里就不在贅述。根據(jù)他的例子，我在下面給出

28、了一段命令流。該方法有個(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，它具體的受力形式有如下兩種：剛性表面邊界（Rigid surface constraint）認(rèn)為接觸面是剛性的，沒(méi)有變形，和通過(guò)節(jié)點(diǎn)耦合命令CERIG比較相似； 分布力邊界（Force-distributed constraint）允

29、許接觸面的變形，和邊界定義命令RBE3相似。使用這種方法，需要用KEYOPT(2) = 2打開(kāi)接觸單元的MPC（多點(diǎn)接觸邊界）算法ANSYS繪制彎矩、剪力、軸力圖1.繪制彎矩圖 建立彎矩單元表。例如梁?jiǎn)卧?#160;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

30、60;   !單元J點(diǎn)剪力ETABLE,MI,SMISC,6     !單元I點(diǎn)彎矩ETABLE,MJ,SMISC,12 !單元J點(diǎn)彎矩plls,MI,MJ 2.標(biāo)注彎矩圖 PLOTCTRLS>>NUMBERING>>SVAL ON即可在畫(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)重影。觀察上面整體軸力圖也

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

32、FORZ,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' !需處理的節(jié)點(diǎn)包 *ASK,NASSEMBLY,'PLEASE INPUT THE COMPONENT

33、 NAME OF NODE TO BE PROCESSED!','NOU TER' !無(wú)需處理的節(jié)點(diǎn)包 *ASK,UNASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NODE!(NONE I F THERE'S NO SUCH COMPONENT)','NONE' /POST1 !輸入信息：內(nèi)力類型，欲處理單元的集合，欲處理節(jié)點(diǎn)的集合 !ITFTYPE='MMOMY' 

34、;!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&#

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

36、BLY,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,1 UNNOD(I+1)=NDNEXT(I) *ENDDO *ELSE UNNOD

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

38、EXT(I) *ENDDO !獲得當(dāng)前選中節(jié)點(diǎn)總數(shù)（存入變數(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,

39、0;!節(jié)點(diǎn)J的內(nèi)力存入數(shù)組ITNFJ ETABLE,ITNFI,SMISC,ITFINUM 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&

40、#160;*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),EQ,NODNUM(K),THEN ELEI=ELENUM(J) *EXIT *ELSE *ENDIF

41、0;*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),THEN ELEJ=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)力存入結(jié)果數(shù)組中 DETAB,ELEI,ITNFINEO,ETAVE DETAB,ELEJ,ITNFJNEO,ETAVE *EL