1、Good is good, but better carries it.精益求精，善益求善。ik4ANSYS中彎矩、剪力圖的繪制2ik4ANSYS中彎矩、剪力圖的繪制2ansys中如何生成命令流方法：GUI是：Utility Menu>File>Write DB Log File怎么用ansys繪制彎矩，剪力圖：GUI: General Postproc-&gtlot Result->Contour Plot->Line Element Result彈出畫單元結(jié)果的對話框，分別在Labi和Labj依次選取SMIS6和SMIS12（彎矩圖）、SMI

2、S1和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點彎矩   

3、;  ! 更新單元表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 /

4、image,save,'Shearing_force_%T%',jpg     ! 畫彎矩分布圖/TITLE,Bending   moment diagramPLLS,MI,MJ,-0.8,0/image,save,'Bending_moment_%T%',jpgANSYS中彎矩、剪力圖的繪制GUI:General Postproc-plot Result-Contour Plot-Line Element Result彈出畫單元結(jié)果的對話框，分別在Labi和Labj依次選取SMIS6和SMIS12（彎

5、矩圖）、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點彎矩 

6、0;   ! 更新單元表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

7、 /image,save,'Shearing_force_%T%',jpg     ! 畫彎矩分布圖/TITLE,Bending   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.

8、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!讀入截面文件secwrite,jiemian,sect,1aclear,alladele,all,1ldele,all,1finish/clear/prep7et,1,beam44keyopt,1,6,1mp,dens,1,2600

9、0;mp,ex,1,3.06e10           mp,prxy,1,0.2      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,1eplo

10、tdk,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é)在回答別人問題時，利用前人的回復(fù)和總結(jié)，自己進行了總結(jié)改正，發(fā)表在這里，供各位

11、參考 （1）ANSYS彎矩等可以直接標(biāo)注在圖上嗎？如何實現(xiàn)？ 如果三維問題，在剖面上標(biāo)出某一結(jié)構(gòu)的軸心力、彎矩等，如何實現(xiàn) （2）后處理圖形，其等值線的數(shù)值能否直接標(biāo)注在圖上，而不是采用圖例的形式 后處理結(jié)果往往用云圖表示，下跟一圖例表示數(shù)值大小，能夠?qū)崿F(xiàn)等值線直接標(biāo)注在圖上 回答 （1） 1.繪制彎矩圖 建立彎矩單元表。例如梁單元 i節(jié)點單元表名稱為imom，j節(jié)點單元表名稱為jmom， ETABLE,NI,SMISC,1     !單元I點軸力ETAB

12、LE,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.標(biāo)注彎矩圖 PLOTCTRLS>>NUMBERING>>SVAL ON即可在畫出彎矩圖

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

14、：!-!宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO） !獲取線性單元內(nèi)力，并對單元邊界處的內(nèi)力進行平衡!輸入信息 !內(nèi)力類型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ *ASK,ITFTYPE,'PLEASE INPUT THE TYPE OF INTERNAL FORCE','MMOMY'!需處理的單元包 *ASK,EASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF ELEMENTS TO BE PROC

15、ESSED!', 'EOUTER'!需處理的節(jié)點包 *ASK,NASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF NODE TO BE PROCESSED!','NOU TER'!無需處理的節(jié)點包 *ASK,UNASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NODE!(NONE I F THERE'S NO SUCH COMPONENT)','NO

16、NE'/POST1!輸入信息：內(nèi)力類型，欲處理單元的集合，欲處理節(jié)點的集合 !ITFTYPE='MMOMY' !EASSEMBLY='EOUTER' !NASSEMBLY='NOUTER'!按內(nèi)力類型確定ANSYS輸出信息SMISC的編號 *IF,ITFTYPE,EQ,'MFORX',THENITFINUM=1 ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2 ITFJNUM=8*ELSEIF,IT

17、FTYPE,EQ,'MFORZ',THENITFINUM=3 ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOMX',THENITFINUM=4 ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5 ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6 ITFJNUM=12*ELSE*ENDIF!對不需平均的節(jié)點進行處理 *IF,UNASSEMBLY,NE,'NO

18、NE',THEN!選出不進行處理的節(jié)點包并獲取不進行處理節(jié)點的數(shù)目 CMSEL,S,UNASSEMBLY *GET,UNNODNUM,NODE,0,COUNT!定義長度為UNNODNUM的數(shù)組（UNNOD），以存放選中單元的單元編號 *DIM,UNNOD,ARRAY,UNNODNUM!將選中單元的編號按順序存入數(shù)組UNNOD *DO,I,0,UNNODNUM-1,1 UNNOD(I+1)=NDNEXT(I) *ENDDO *ELSE UNNODNUM=0 *ENDIF!選出所需的單元和節(jié)點包

19、60;CMSEL,S,EASSEMBLY CMSEL,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,1 ELENUM(I+1)=ELNEXT(I) *ENDDO!獲得當(dāng)前選中節(jié)點總數(shù)（存入變量SELNODNUM） *GET,S

20、ELNODNUM,NODE,0,COUNT!定義長度為SELNODNUM的數(shù)組（NODNUM），以存放選中單元的單元編號 *DIM,NODNUM,ARRAY,SELNODNUM!將選中單元的編號按順序存入數(shù)組NODNUM *DO,I,0,SELNODNUM-1,1 NODNUM(I+1)=NDNEXT(I) *ENDDO!定義所需的線性單元內(nèi)力ETABLE，節(jié)點I的內(nèi)力存入數(shù)組ITNFI, !節(jié)點J的內(nèi)力存入數(shù)組ITNFJ ETABLE,ITNFI,SMISC,ITFINUM ETABLE,ITNFJ,SMISC,ITFJN

21、UM!定義所需的結(jié)果數(shù)組，并將其置零 ETABLE,ITNFINEO,SMISC,5 SADD,ITNFINEO,ITNFI,1 ETABLE,ITNFJNEO,SMISC,11 SADD,ITNFJNEO,ITNFJ,1*DO,K,1,SELNODNUM,1!處理不需平均的節(jié)點 INDEX=0 *IF,UNNODNUM,GE,1,THEN *DO,J,1,UNNODNUM*IF,NODNUM(K),EQ,UNNOD(J),THEN INDEX=1 *ELSE *ENDIF*ENDDO 

22、;*ELSE *ENDIF*DO,J,1,SELELENUM,1!選出和節(jié)點K相連的線性單元中，I節(jié)點（對線性單元而言）為節(jié)點K的單元編號 *IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THEN ELEI=ELENUM(J) *EXIT *ELSE *ENDIF*ENDDO*DO,J,1,SELELENUM,1!選出和節(jié)點K相連的線性單元中，J節(jié)點（對線性單元而言）為節(jié)點K的單元編號 *IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THEN ELEJ=ELENUM(

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

24、!將平均后的內(nèi)力存入結(jié)果數(shù)組中 DETAB,ELEI,ITNFINEO,ETAVE DETAB,ELEJ,ITNFJNEO,ETAVE*ELSE *ENDIF*ELSE *ENDIF*ENDDO/UDOC,1,LOGO,OFF PLLS,ITNFINEO,ITNFJNEO!END OF ITFAVG.MAC （2） 對體和面來說，ANSYS默認(rèn)的結(jié)果輸出格式是云圖格式，而這種彩色云圖打印為黑白圖像時對比很不明顯，無法表達清楚，對于發(fā)表文章非常不便。發(fā)文章所用的結(jié)果圖最好是等值線圖，并且最好是黑白的等值線圖。一般借用photo

25、shop等第三方軟件，很麻煩，效果也不好。ansys本身也能實現(xiàn)這項功能。 1.將要輸出的結(jié)果調(diào)出，這時為彩色云圖； 2.將云圖轉(zhuǎn)換為等值線圖的形式 GUI：plotCtrls>Device Options>/DEVI中的vector mode 選為on 命令：/DEVICE,VECTOR,1 這時結(jié)果為彩色等值線，若直接輸出，打印為黑白圖像時仍然不清晰，為此需進行以下幾步將圖像轉(zhuǎn)換為黑白形式； 3. 將背景變?yōu)榘咨?#160;命令：jpgprf,500,100,1 /rep 4.對等值線中的等值線符號

26、（圖中為A,B,C等）的疏密進行調(diào)整 GUI：plotCtrls>Style>Contours> Contours Labeling 在Key Vector mode contour label 中選中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ù)選框

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

28、L,1 /replot 如果還想將等值線的字母去掉 /PNUM,SVAL,1 /CLABEL,1,-1 /replot 7.出圖 GUI：plotCtrls>Capture ImageANSYS中如何加彎矩或扭矩A轉(zhuǎn)矩一般有三種施加的方法: 第一種,將矩轉(zhuǎn)換成一對一對的力偶,直接施加在對應(yīng)的節(jié)點上面.第二種,在構(gòu)件中心部位建立一個節(jié)點,定義為MASS21單元,然后跟其他受力節(jié)點藕荷,形成剛性區(qū)域,就是用CERIG命令.然后直接加轉(zhuǎn)矩到主節(jié)點,即中心節(jié)點上面第三種,使用MPC184單元.是在構(gòu)件中心部位建立一個

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

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

31、標(biāo)，然后等效為施加環(huán)向的節(jié)點力；而施加彎矩，可以將力矩轉(zhuǎn)化為端面的剪切均布力；但這種方法比較容易出現(xiàn)應(yīng)力集中現(xiàn)象；方法2，定義局部剛性區(qū)域，施加過程venture講的很詳細(xì)，這里就不在贅述。根據(jù)他的例子，我在下面給出了一段命令流。該方法有個不足，它在端面額外的增加了一定的剛度，只能適用于小變形分析。方法3，相對方法2來說，采用剛性梁單元，適用范圍更廣一些，對于大應(yīng)變分析也能很好的適用。但在小應(yīng)變分析下，方法2和方法3沒有什么區(qū)別。方法4，定義一個主節(jié)點，施加了分布力面，應(yīng)該說跟實際比較接近一點，但端面的結(jié)果好像不是很理想，結(jié)果有點偏大，在遠(yuǎn)離端面處的位置跟實際很符合。方法5，它具體的受力形式有

32、如下兩種：剛性表面邊界（Rigid surface constraint）認(rèn)為接觸面是剛性的，沒有變形，和通過節(jié)點耦合命令CERIG比較相似； 分布力邊界（Force-distributed constraint）允許接觸面的變形，和邊界定義命令RBE3相似。使用這種方法，需要用KEYOPT(2) = 2打開接觸單元的MPC（多點接觸邊界）算法ANSYS繪制彎矩、剪力、軸力圖1.繪制彎矩圖 建立彎矩單元表。例如梁單元 i節(jié)點單元表名稱為imom，j節(jié)點單元表名稱為jmom， ETABLE,NI,SMISC,1     !單元I點軸力E

33、TABLE,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.標(biāo)注彎矩圖 PLOTCTRLS>>NUMBERING>>SVAL ON即可在畫出彎矩圖的同時在圖上標(biāo)出彎矩值的大小 3

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

35、!宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO） !獲取線性單元內(nèi)力，并對單元邊界處的內(nèi)力進行平衡 !輸入信息 !內(nèi)力類型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ *ASK,ITFTYPE,'PLEASE INPUT THE TYPE OF INTERNAL FORCE','MMOMY' !需處理的單元包 *ASK,EASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF ELEMENTS T

36、O BE PROCESSED!', 'EOUTER' !需處理的節(jié)點包 *ASK,NASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF NODE TO BE PROCESSED!','NOU TER' !無需處理的節(jié)點包 *ASK,UNASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF THE UNCHANGED NODE!(NONE I F THERE'S NO SUCH CO

37、MPONENT)','NONE' /POST1 !輸入信息：內(nèi)力類型，欲處理單元的集合，欲處理節(jié)點的集合 !ITFTYPE='MMOMY' !EASSEMBLY='EOUTER' !NASSEMBLY='NOUTER' !按內(nèi)力類型確定ANSYS輸出信息SMISC的編號 *IF,ITFTYPE,EQ,'MFORX',THEN ITFINUM=1 ITFJNUM=7 *ELSEIF,ITFTYPE,EQ,

38、9;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,&

39、#39;MMOMZ',THEN ITFINUM=6 ITFJNUM=12 *ELSE *ENDIF !對不需平均的節(jié)點進行處理 *IF,UNASSEMBLY,NE,'NONE',THEN !選出不進行處理的節(jié)點包并獲取不進行處理節(jié)點的數(shù)目 CMSEL,S,UNASSEMBLY *GET,UNNODNUM,NODE,0,COUNT !定義長度為UNNODNUM的數(shù)組（UNNOD），以存放選中單元的單元編號 *DIM,UNNOD,ARRAY,UNNODNUM

40、0;!將選中單元的編號按順序存入數(shù)組UNNOD *DO,I,0,UNNODNUM-1,1 UNNOD(I+1)=NDNEXT(I) *ENDDO *ELSE UNNODNUM=0 *ENDIF !選出所需的單元和節(jié)點包 CMSEL,S,EASSEMBLY CMSEL,S,NASSEMBLY !獲得當(dāng)前選中單元總數(shù)（存入變量SELELENUM） *GET,SELELENUM,ELEM,0,COUNT !定義長度為SELELENUM的數(shù)組（ELENUM），以存放選中單元的單元編號

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

42、160;*DO,I,0,SELNODNUM-1,1 NODNUM(I+1)=NDNEXT(I) *ENDDO !定義所需的線性單元內(nèi)力ETABLE，節(jié)點I的內(nèi)力存入數(shù)組ITNFI, !節(jié)點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&

43、#160;SADD,ITNFJNEO,ITNFJ,1 *DO,K,1,SELNODNUM,1 !處理不需平均的節(jié)點 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é)點K相連的線性單元中，I節(jié)點（對線性單元而言）為節(jié)點K

44、的單元編號 *IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THEN ELEI=ELENUM(J) *EXIT *ELSE *ENDIF *ENDDO *DO,J,1,SELELENUM,1 !選出和節(jié)點K相連的線性單元中，J節(jié)點（對線性單元而言）為節(jié)點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 !有可能出現(xiàn)ELEJ為0的情況 !取出I節(jié)點為節(jié)點K的單元的I節(jié)點端的內(nèi)力放入?yún)?shù)ETELEI *GET,ETELEI,ELEM,ELEI,ETAB,ITNFI !取出J節(jié)點為節(jié)點K的單元的J節(jié)點端的內(nèi)力放入?yún)?shù)ETELEJ *GET,ETELEJ,ELEM,ELEJ,ETAB,ITNFJ !平均節(jié)點K的單元的I節(jié)點端的內(nèi)力和節(jié)點K的單元的J節(jié)點端的內(nèi)力 ETAVE=(ETELEI+ET