ANSYS中彎矩、剪力圖的繪制_第1頁
ANSYS中彎矩、剪力圖的繪制_第2頁
ANSYS中彎矩、剪力圖的繪制_第3頁
ANSYS中彎矩、剪力圖的繪制_第4頁
ANSYS中彎矩、剪力圖的繪制_第5頁
已閱讀5頁,還剩13頁未讀 繼續免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

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

2、#160;   !單元J點軸力ETABLE,QI,SMISC,2     !單元I點剪力ETABLE,QJ,SMISC,8     !單元J點剪力ETABLE,MI,SMISC,6     !單元I點彎矩ETABLE,MJ,SMISC,12 !單元J點彎矩     ! 更新單元表ETABLE,REFL     ! 畫軸力分布圖/TITLE,Axial   force  

3、; 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 

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中彎矩圖,云圖繪制總結在回答別人問題時,利用前人的回復和總結,自己進行了總結改正,發表在這里,供各位參考 (1)ANSYS彎矩等可以直接標注在圖上嗎?如何實現? 如果三維問題,在剖面上標出某一結構的軸心力、彎矩等,如何實現 (2)后處理圖形,其等值線的數值能否直接標注在圖上,而不是采

8、用圖例的形式 后處理結果往往用云圖表示,下跟一圖例表示數值大小,能夠實現等值線直接標注在圖上 回答 (1) 1.繪制彎矩圖 建立彎矩單元表。例如梁單元 i節點單元表名稱為imom,j節點單元表名稱為jmom, ETABLE,NI,SMISC,1     !單元I點軸力ETABLE,NJ,SMISC,7     !單元J點軸力ETABLE,QI,SMISC,2     !單元I點剪力ETABLE,QJ,SMI

9、SC,8     !單元J點剪力ETABLE,MI,SMISC,6     !單元I點彎矩ETABLE,MJ,SMISC,12 !單元J點彎矩plls,imom,jmom 2.標注彎矩圖 PLOTCTRLS>>NUMBERING>>SVAL ON即可在畫出彎矩圖的同時在圖上標出彎矩值的大小 3.調整彎矩圖 如果彎矩圖方向錯誤,則繪制彎矩圖命令為 plls,imom,jmom,-1 同一個節點處兩邊的單元內力有細微差別, 導

10、致內力數字標注出現重影。觀察上面整體軸力圖也可以發現, 一段一段的,好像馬賽克,其實上面整體彎矩圖也是,不過不是 很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成 的(詳細原因就不展開了,看看梁理論的書和ANSYS的說明吧)。為了修正重影和節點兩邊內力值不一樣的問題,遍制了宏文件ITFAVG.MAC 命令檔內容如下:!-!宏:ITFAVG.MAC(INTERNAL FORCE AVERAGE MACRO) !獲取線性單元內力,并對單元邊界處的內力進行平衡!輸入信息 !內力類型: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'!需處理的節點包 *ASK,NASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF NODE

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

13、0;!NASSEMBLY='NOUTER'!按內力類型確定ANSYS輸出信息SMISC的編號 *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!對不需平均的節點進行處理 *IF,UNASSEMBLY,NE,'NONE',THEN!選出不進行處理的節點包并獲取不進行處理節點的數目 CMSEL,S,UNASSEMBLY *GET,UNNODNUM,NODE,0,COUNT!定義長度為UNNODNUM的

15、數組(UNNOD),以存放選中單元的單元編號 *DIM,UNNOD,ARRAY,UNNODNUM!將選中單元的編號按順序存入數組UNNOD *DO,I,0,UNNODNUM-1,1 UNNOD(I+1)=NDNEXT(I) *ENDDO *ELSE UNNODNUM=0 *ENDIF!選出所需的單元和節點包 CMSEL,S,EASSEMBLY CMSEL,S,NASSEMBLY!獲得當前選中單元總數(存入變數SELELENUM) *GET,SELELENUM,ELEM,0,COUNT!定義長度為

16、SELELENUM的數組(ELENUM),以存放選中單元的單元編號 *DIM,ELENUM,ARRAY,SELELENUM!將選中單元的編號按順序存入數組ELENUM *DO,I,0,SELELENUM-1,1 ELENUM(I+1)=ELNEXT(I) *ENDDO!獲得當前選中節點總數(存入變數SELNODNUM) *GET,SELNODNUM,NODE,0,COUNT!定義長度為SELNODNUM的數組(NODNUM),以存放選中單元的單元編號 *DIM,NODNUM,ARRAY,SELNODNUM!將選中單元的編號按順序存入數

17、組NODNUM *DO,I,0,SELNODNUM-1,1 NODNUM(I+1)=NDNEXT(I) *ENDDO!定義所需的線性單元內力ETABLE,節點I的內力存入數組ITNFI, !節點J的內力存入數組ITNFJ ETABLE,ITNFI,SMISC,ITFINUM ETABLE,ITNFJ,SMISC,ITFJNUM!定義所需的結果數組,并將其置零 ETABLE,ITNFINEO,SMISC,5 SADD,ITNFINEO,ITNFI,1 ETABLE,ITNFJNEO,SMISC,11

18、0;SADD,ITNFJNEO,ITNFJ,1*DO,K,1,SELNODNUM,1!處理不需平均的節點 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!選出和節點K相連的線性單元中,I節點(對線性單元而言)為節點K的單元編號 *IF,NELEM(ELENUM(J),1),EQ,NO

19、DNUM(K),THEN ELEI=ELENUM(J) *EXIT *ELSE *ENDIF*ENDDO*DO,J,1,SELELENUM,1!選出和節點K相連的線性單元中,J節點(對線性單元而言)為節點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 !有可能出現ELEJ為0的情況!取出I節點為節點K的單元

20、的I節點端的內力放入參數ETELEI *GET,ETELEI,ELEM,ELEI,ETAB,ITNFI !取出J節點為節點K的單元的J節點端的內力放入參數ETELEJ *GET,ETELEJ,ELEM,ELEJ,ETAB,ITNFJ!平均節點K的單元的I節點端的內力和節點K的單元的J節點端的內力 ETAVE=(ETELEI+ETELEJ)/2!將平均后的內力存入結果數組中 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) 對體和面來說,ANSYS默認的結果輸出格式是云圖格式,而這種彩色云圖打印為黑白圖像時對比很不明顯,無法表達清楚,對于發表文章非常不便。發文章所用的結果圖最好是等值線圖,并且最好是黑白的等值線圖。一般借用photoshop等第三方軟件,很麻煩,效果也不好。ansys本身也能實現這項功能。 1.將要輸出的結果調出,這時為彩色云圖; 2.將云圖轉換為等值線圖的形式 GUI:plotCtrls>D

22、evice Options>/DEVI中的vector mode 選為on 命令:/DEVICE,VECTOR,1 這時結果為彩色等值線,若直接輸出,打印為黑白圖像時仍然不清晰,為此需進行以下幾步將圖像轉換為黑白形式; 3. 將背景變為白色 命令:jpgprf,500,100,1 /rep 4.對等值線中的等值線符號(圖中為A,B,C等)的疏密進行調整 GUI:plotCtrls>Style>Contours> Contours Labeling 在Key Vector mode contour la

23、bel 中選中on every Nth elem,然后在N= 輸入框中輸入合適的數值,例如5,多試幾次,直到疏密合適 命令:/clabel,1,5 5.將彩色等值線變為黑色 GUI:plotCtrls>Style>Colors>Contours Colors 將Items Numbered 1,Items Numbered 2等復選框中的顏色均選為黑色,圖像即可變為黑白等值線圖像 命令:/color,cntr,whit,1 等等 6.修改 在顯示等值線時,圖形中太多的標識字母(A、B等)使得等值線圖完全看不清楚。請問

24、有沒有什么方法可以去掉圖形的這些標識字母呢 在plotCtrls下面的style下面的contours下面的contour Labeling,在彈出的contour labeling options里面的vector mode contour labels選off 后處理中畫出了結構的等值線圖,可否將等值線邊的字母改為相應的數值標注 /PNUM,SVAL,1 /replot 如果還想將等值線的字母去掉 /PNUM,SVAL,1 /CLABEL,1,-1 /replot 7.出圖 GUI:plot

25、Ctrls>Capture ImageANSYS中如何加彎矩或扭矩A轉矩一般有三種施加的方法: 第一種,將矩轉換成一對一對的力偶,直接施加在對應的節點上面.第二種,在構件中心部位建立一個節點,定義為MASS21單元,然后跟其它受力節點藕荷,形成剛性區域,就是用CERIG命令.然后直接加轉矩到主節點,即中心節點上面第三種,使用MPC184單元.是在構件中心部位建立一個節點,跟其它受力節點分別形成多根剛性梁,從而形成剛性面.最后也是直接加載荷到中心節點上面,通過剛性梁來傳遞載荷.上面三種方法計算的結果基本一致,我做過實驗的.只不過是后兩種情況都是形成剛性區域,但是CERIG命令是要

26、在小變形或者小旋轉才能用,只支持靜力,線形分析.而第三種方法適用多種情況,不僅支持大應變,還支持非線形情況.如果你需要例子,我下次在發給你看A思路1:矩或扭矩說白了就是矩,所謂矩就是力和力臂的乘積。         施加矩可以等效為施加力;思路2:直接施加彎矩或扭矩,此時需要引入一個具有旋轉自由度的節點;可以選擇單元21,或者1841將矩轉換成一對的力偶,直接施加在對應的節點上面。2在構件中心部位建立一個節點,定義為mass21單元,然后跟其它受力節點耦合,形成剛性區域,就是用cerig命令。然后直接加轉矩到主節點,即中心節點上面。3

27、使用mpc184單元。是在構件中心部位建立一個節點,跟其它受力節點分別形成多根剛性梁,從而形成剛性面。最后也是直接加載荷到中心節點上面,通過剛性梁來傳遞載荷。4通過rbe3命令。該方法與方法2很接近。5基于表面邊界法:主要通過定義一個接觸表面和一個目標節點接觸來實現,彎矩荷載可以通過在目標節點上用“F”命令施加。對于方法1,通過轉換為集中力或均布力,比如施加扭矩,把端面節點改成柱坐標,然后等效為施加環向的節點力;而施加彎矩,可以將力矩轉化為端面的剪切均布力;但這種方法比較容易出現應力集中現象;方法2,定義局部剛性區域,施加過程venture講的很詳細,這里就不在贅述。根據他的例子,我在下面給出

28、了一段命令流。該方法有個不足,它在端面額外的增加了一定的剛度,只能適用于小變形分析。方法3,相對方法2來說,采用剛性梁單元,適用范圍更廣一些,對于大應變分析也能很好的適用。但在小應變分析下,方法2和方法3沒有什么區別。方法4,定義一個主節點,施加了分布力面,應該說跟實際比較接近一點,但端面的結果好像不是很理想,結果有點偏大,在遠離端面處的位置跟實際很符合。方法5,它具體的受力形式有如下兩種:剛性表面邊界(Rigid surface constraint)認為接觸面是剛性的,沒有變形,和通過節點耦合命令CERIG比較相似; 分布力邊界(Force-distributed constraint)允

29、許接觸面的變形,和邊界定義命令RBE3相似。使用這種方法,需要用KEYOPT(2) = 2打開接觸單元的MPC(多點接觸邊界)算法ANSYS繪制彎矩、剪力、軸力圖1.繪制彎矩圖 建立彎矩單元表。例如梁單元 i節點單元表名稱為imom,j節點單元表名稱為jmom, ETABLE,NI,SMISC,1     !單元I點軸力ETABLE,NJ,SMISC,7     !單元J點軸力ETABLE,QI,SMISC,2     !單元I點剪力ETABLE,QJ,SMISC,8 

30、60;   !單元J點剪力ETABLE,MI,SMISC,6     !單元I點彎矩ETABLE,MJ,SMISC,12 !單元J點彎矩plls,MI,MJ 2.標注彎矩圖 PLOTCTRLS>>NUMBERING>>SVAL ON即可在畫出彎矩圖的同時在圖上標出彎矩值的大小 3.調整彎矩圖 如果彎矩圖方向錯誤,則繪制彎矩圖命令為 plls,imom,jmom,-1 同一個節點處兩邊的單元內力有細微差別, 導致內力數字標注出現重影。觀察上面整體軸力圖也

31、可以發現, 一段一段的,好像馬賽克,其實上面整體彎矩圖也是,不過不是 很明顯罷了。這是EULER-BEONOULI梁理論以及ANSYS輸出定義造成 的(詳細原因就不展開了,看看梁理論的書和ANSYS的說明吧)。 為了修正重影和節點兩邊內力值不一樣的問題,遍制了宏文件ITFAVG.MAC 命令檔內容如下: !- !宏:ITFAVG.MAC(INTERNAL FORCE AVERAGE MACRO) !獲取線性單元內力,并對單元邊界處的內力進行平衡 !輸入信息 !內力類型: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' !需處理的節點包 *ASK,NASSEMBLY,'PLEASE INPUT THE COMPONENT

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

34、;!EASSEMBLY='EOUTER' !NASSEMBLY='NOUTER' !按內力類型確定ANSYS輸出信息SMISC的編號 *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 !對不需平均的節點進行處理 *IF,UNASSEM

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

37、NUM=0 *ENDIF !選出所需的單元和節點包 CMSEL,S,EASSEMBLY CMSEL,S,NASSEMBLY !獲得當前選中單元總數(存入變數SELELENUM) *GET,SELELENUM,ELEM,0,COUNT !定義長度為SELELENUM的數組(ELENUM),以存放選中單元的單元編號 *DIM,ELENUM,ARRAY,SELELENUM !將選中單元的編號按順序存入數組ELENUM *DO,I,0,SELELENUM-1,1 ELENUM(I+1)=ELN

38、EXT(I) *ENDDO !獲得當前選中節點總數(存入變數SELNODNUM) *GET,SELNODNUM,NODE,0,COUNT !定義長度為SELNODNUM的數組(NODNUM),以存放選中單元的單元編號 *DIM,NODNUM,ARRAY,SELNODNUM !將選中單元的編號按順序存入數組NODNUM *DO,I,0,SELNODNUM-1,1 NODNUM(I+1)=NDNEXT(I) *ENDDO !定義所需的線性單元內力ETABLE,節點I的內力存入數組ITNFI,

39、0;!節點J的內力存入數組ITNFJ ETABLE,ITNFI,SMISC,ITFINUM ETABLE,ITNFJ,SMISC,ITFJNUM !定義所需的結果數組,并將其置零 ETABLE,ITNFINEO,SMISC,5 SADD,ITNFINEO,ITNFI,1 ETABLE,ITNFJNEO,SMISC,11 SADD,ITNFJNEO,ITNFJ,1 *DO,K,1,SELNODNUM,1 !處理不需平均的節點 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 !選出和節點K相連的線性單元中,I節點(對線性單元而言)為節點K的單元編號 *IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THEN ELEI=ELENUM(J) *EXIT *ELSE *ENDIF

41、0;*ENDDO *DO,J,1,SELELENUM,1 !選出和節點K相連的線性單元中,J節點(對線性單元而言)為節點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 !有可能出現ELEJ為0的情況 !取出I節點為節點K的單元的I節點端的內力放入參數ETELEI *GET,ETELEI,ELEM,ELEI,ETAB,ITNFI !取出J節點為節點K的單元的J節點端的內力放入參數ETELEJ *GET,ETELEJ,ELEM,ELEJ,ETAB,ITNFJ !平均節點K的單元的I節點端的內力和節點K的單元的J節點端的內力 ETAVE=(ETELEI+ETELEJ)/2 !將平均后的內力存入結果數組中 DETAB,ELEI,ITNFINEO,ETAVE DETAB,ELEJ,ITNFJNEO,ETAVE *EL

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論