有限元計(jì)算注意事項(xiàng)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上1、計(jì)算原理任何有限元模擬的第一步都是利用一個(gè)有限單元的集合離散結(jié)構(gòu)的實(shí)際幾何形狀,每個(gè)單元（element）代表這個(gè)實(shí)際結(jié)構(gòu)的一個(gè)離散部分.這些單元通過公用的節(jié)點(diǎn)（node）來連接.節(jié)點(diǎn)和單元的集合成為網(wǎng)格（mesh）.在一個(gè)特定網(wǎng)格中的單元數(shù)目稱為網(wǎng)格密度（mesh density）.在ansys計(jì)算過程中,程序以每個(gè)節(jié)點(diǎn)的每個(gè)自由度建立平衡方程,以節(jié)點(diǎn)的位移作為未知量,利用矩陣求解節(jié)點(diǎn)的位移.一旦節(jié)點(diǎn)位移求出,整個(gè)結(jié)構(gòu)的應(yīng)力和應(yīng)變都很容易計(jì)算出來.這種計(jì)算的過程和方法,數(shù)學(xué)上稱之為隱式方法.從上敘述來看,整個(gè)計(jì)算過程中就是求解n個(gè)n元一次方程組（n表示節(jié)點(diǎn)數(shù)量）

2、,當(dāng)計(jì)算模型復(fù)雜而且龐大時(shí),隱式求解方法的計(jì)算量還是相當(dāng)大的.與之相對(duì)應(yīng)的,顯式求解方法.顯式求解方法是通過動(dòng)態(tài)方法從一個(gè)增量步前推到下一個(gè)增量步得到的.具體顯式求解方法和隱式求解方法例子如下：（1） 隱式求解（2） 顯式求解隱式求解中,計(jì)算的精度完全控制于計(jì)算步數(shù),在一般的計(jì)算軟件中（flac、abaqus）,軟件均是利用不平衡力來控制計(jì)算步數(shù)（當(dāng)不平衡力<10-5時(shí),停止計(jì)算）.不平衡力A+B.A表示施加在節(jié)點(diǎn)上的集中力；B表示：在n步數(shù)下,根據(jù)第n步計(jì)算出來的應(yīng)力,求出節(jié)點(diǎn)的內(nèi)力. Flac軟件中 B,以上公式是根據(jù)虛功原理推倒而得到.具體推倒過程見flac 原理.2、Ansys計(jì)

3、算注意事項(xiàng)：計(jì)算單位、參數(shù)、荷載、標(biāo)準(zhǔn)值、設(shè)計(jì)值,計(jì)算過程中系數(shù)的加入.（1） beam單元對(duì)于beam單元.Ansys軟件中我們常用的有兩種梁?jiǎn)卧篵eam188和beam4.這兩種單元均是三維的梁?jiǎn)卧?每個(gè)節(jié)點(diǎn)都具有6個(gè)自由度（ux、uy、uz、mx、my、mz）,并且單元坐標(biāo)系x軸是i點(diǎn)指向j點(diǎn).Beam188單元是基于鐵木辛哥理論的梁,beam4單元是我們常用的經(jīng)典結(jié)構(gòu)力學(xué)梁.（鐵木辛哥理論考慮了梁的剪切變形,而我們常用的經(jīng)典結(jié)構(gòu)力學(xué)梁只考慮了彎矩對(duì)結(jié)構(gòu)的變形影響）所以說,beam188可以更精確的計(jì)算梁?jiǎn)卧?因此我們結(jié)構(gòu)計(jì)算中,一般都采用beam188單元.當(dāng)然還有beam189單元

4、,189單元屬于三維二次的梁?jiǎn)卧╞eam188屬于三維一次梁?jiǎn)卧?精度比beam188更加高.定義beam188單元,一般采用如下形式：!定義單元/prep7 !進(jìn)入前處理et,1,beam188 !定義單元188號(hào)標(biāo)號(hào)為1!定義材料屬性mp,ex,1,2.55e7 !定義彈性模量(kn/m2)mp,nuxy,1,0.167 !定義泊松比mp,dens,1,2.5 !定義密度（KN/N*KG/M3）nummrg,all ！合并重合節(jié)點(diǎn)numcmp,all ！壓縮編號(hào)!定義梁截面SECTYPE,1,BEAM,RECT,A1,0 ! 1表示梁編號(hào) ; RECT表示是矩形梁（還有其他t型等等,具

5、體見ansys幫助）; A1 表示梁的名稱 ; 0表示薄壁梁?jiǎn)卧W(wǎng)格劃分精細(xì)程度（05）.SECDATA,1,3,4,12 !1表示梁b ; 3表示梁h ; 4和12定義對(duì)應(yīng)寬長(zhǎng)等分份數(shù).SECOFFSET,CENT !cent質(zhì)心 ; shrc剪切中心 ; origin原始中心 ; user用戶定義;!注意：當(dāng)梁?jiǎn)卧蜌卧黄鹗褂脮r(shí),可以設(shè)置梁?jiǎn)卧钠?使梁的一面和殼的一面共面.（secoffset,user,offsety,offsetz）,如下圖：!劃分網(wǎng)格LSEL,S,1 !選中編號(hào)為1的線.LATT,1,1,1,1 !mp,r,et,方向點(diǎn),SECTYPE截面號(hào).LESIZE,A

6、LL,0.2,1 !0.2是單元大小,1是確認(rèn)細(xì)分規(guī)則.LMESH,ALL !用beam單元離散模型,形成網(wǎng)格.!對(duì)于劃分網(wǎng)格,空間的beam單元,由于需要確定b、h的方向,ansys軟件利用方向點(diǎn)來控制b、h的方向.方向點(diǎn)的編號(hào)最好定義的很大,如果定義太小,會(huì)影響后面的加載.具體方向點(diǎn)如何控制見上面的latt命令和ansys幫助.自己試兩下就知道怎么用了.AllsFINISH!加載加約束/SOLUACEL,9.8 !重力加速度.注意方向,數(shù)值和整體坐標(biāo)相反,比如重力指向z軸負(fù)向,則為正值.SFCUM,ALL,ADD !設(shè)置單元荷載是疊加還是替代,只對(duì)加在單元和節(jié)點(diǎn)上的荷載有效,對(duì)于加在面、線

7、上的荷載,都只有替代作用（對(duì)同一個(gè)面,第二次加的荷載替代第一次加的荷載）!對(duì)于beam單元,只能根據(jù)sfbeam命令增加均布荷載 等大小的均布荷載.Lsel,s,1ESLL,S,1sfbeam,ALL,1,PRES,-161.5 !1表示作用在beam單元的面上(如下圖,面表示beam單元的軸向,面表示單元側(cè)面,面表示beam單元頂面),-161.5表示均布荷載大小,正負(fù)號(hào)可以控制作用力的方向. 梯形均布荷載Sfbeam命令是對(duì)每個(gè)單元進(jìn)行加載.如果一根梁承受10100的梯形均布荷載,而且這根梁被分成了10個(gè)beam單元,這樣施加荷載就非常困難.因此我將這種加載過程寫成命令流,讓軟件自動(dòng)進(jìn)行加

8、載.命令流如下：LSEL,S,1 !選中要加載的那根梁（線）ESLL,S,1 !選中屬于這根梁（線）的beam單元*GET,Nelem,ELEM,COUNT, , , , !獲得當(dāng)前所選單元個(gè)數(shù),賦予參數(shù)Nelem*GET,Ne,ELEM,NUM,MIN, , , , !獲得當(dāng)前所選單元最小編號(hào),賦予參數(shù)Ne*DO,I,1,Nelem !循環(huán)加載,循環(huán)次數(shù)單元個(gè)數(shù)ESEL,S,NeNSLE,S,1*GET,Nnode,NODE,COUNT, , , , !獲得當(dāng)前所選節(jié)點(diǎn)個(gè)數(shù),賦予參數(shù)Nnode*GET,Nn,NODE,NUM,MIN, , , , !獲得當(dāng)前所選節(jié)點(diǎn)最小編號(hào),賦予參數(shù)NnNN

9、1X=NX(NN) !將nn節(jié)點(diǎn)的x坐標(biāo)賦予NN1X(NX表示x坐標(biāo),NY表示y坐標(biāo))NN=NDNEXT(NN) !NN=當(dāng)前所選節(jié)點(diǎn)的下一個(gè)編號(hào)NN2X=NX(NN) !將nn節(jié)點(diǎn)的x坐標(biāo)賦予NN2Xsfbeam,ALL,1,PRES,-1630.76/3.23*(3.23-NN1X),-1630.76/3.23*(3.23-NN2X)!以上荷載公式應(yīng)根據(jù)實(shí)際情況進(jìn)行調(diào)整LSEL,S,1ESLL,S,1NE=ELnext(NE) !NE=當(dāng)前所選單元的下一個(gè)編號(hào) *ENDDO!對(duì)于此命令流,根據(jù)不同的實(shí)際情況,ABC部分需要修改,其他不需要修改.!后處理 (XY平面) （大拇指指向y，就是m

10、y）etable,ImY,smisc,2 !顯示彎距etable,JmY,smisc,15PLLS,IMY,JMYETABLE,IFX,SMISC,1 !顯示軸力ETABLE,JFX,SMISC,14PLLS,IFX,JFXETABLE,IFY,SMISC,5 !顯示剪力ETABLE,JFY,SMISC,18PLLS,IFY,JFY!注意：beam單元的結(jié)果輸出都是以單元坐標(biāo)系輸出的,且拉為正、壓為負(fù).前面我們已經(jīng)知道,單元坐標(biāo)系x軸就是i點(diǎn)指向j點(diǎn),其他坐標(biāo)可以根據(jù)整體坐標(biāo)系推出.詳細(xì)內(nèi)容見ansys 幫助.（2） Shell單元對(duì)于shell單元應(yīng)用的范圍,ansys軟件并沒有強(qiáng)制規(guī)定,只

11、是從字面上區(qū)分了薄殼和厚殼.我以前看過一本電子教案仿真在線,里面說一般規(guī)定殼體的主尺寸是厚度的10倍左右,都是可以用殼體來模擬的.一般高度與跨度之比（非與單元尺寸比較）<1/15,可以當(dāng)作薄殼處理,>1/15 & <1/10,可以當(dāng)作厚殼來處理.shell63是薄殼單元,他包含彎曲和薄膜效應(yīng),但是忽略橫向剪切變形；shell43,shell143,shell181,shell91,shell93和shell99,都屬于厚殼單元,不僅有彎曲、薄膜效應(yīng),他也包含了橫向剪切效應(yīng)橫向剪力所引起的變形和彎曲變形屬同一量級(jí),中厚板應(yīng)該考慮.橫向剪切被表示為整個(gè)厚度上的常剪切應(yīng)變.

12、這種一階近似只適用于中等厚度殼體. 線形分析時(shí),如果不包含橫向剪切應(yīng)變,使用63,163單元;如果橫向剪切變形重要,則遵守以下原則:均勻材料,使用43,93,143單元,復(fù)合材料使用91,99,181.我們土木工程中,一般利用shell43計(jì)算.!定義單元/prep7 !進(jìn)入前處理et,1,shell43 !定義單元43號(hào)標(biāo)號(hào)為1!定義材料屬性mp,ex,1,2.55e7 !定義彈性模量(kn/m2)mp,nuxy,1,0.167 !定義泊松比mp,dens,1,2.5 !定義密度（KN/N*KG/M3）!定義墻體厚度!等厚度板R,1,2 !1表示編號(hào),2表示厚度（m）R,2,3Asel,s,

13、1 !選中1號(hào)面Aatt,1,1,1 !mp,real,typeESIZE,0.2 !定義單元大小為0.2左右MSHAPE,0,2D !規(guī)定劃分單元形狀,0表示四邊形（1表示三角形）,2d表示劃分面（3d表示劃分體）MSHKEY,2 !指定是自由劃分還是映射劃分,2表示：盡量用映射劃分,不符合要求就自動(dòng)使用自由劃分,具體參見ansys幫助的eshkey命令.AMESH,ALL !劃分面單元.！注意：在網(wǎng)格剖面方面，最好全部用四邊形，而且形狀盡量規(guī)則、均勻！因?yàn)閷砗筇幚韮?nèi)力提取的時(shí)候，提取出來的力和單元的大小有直接的關(guān)系。后處理部分將詳細(xì)講解。!漸變厚度的板對(duì)于漸變的板,首先選中要?jiǎng)澐值拿?然

14、后進(jìn)行網(wǎng)格劃分.但是在網(wǎng)格劃分的時(shí)候,賦予的厚度編號(hào)開始選用100（前面沒有定義的就可以）,然后利用以下命令進(jìn)行厚度賦值.ASEL,S,1ESLAMXNODE = NDINQR(0,14)*DIM,THICK,MXNODE*DO,NODE,1,MXNODE*IF,NSEL(NODE),EQ,1,THENTHICK(node) = 3-(NZ(NODE)-7.5)/25.5*1.5 !厚度和點(diǎn)坐標(biāo)系的關(guān)系,NZ(NODE)表示編號(hào)為node的點(diǎn)的Z坐標(biāo)值.*ENDIF*ENDDONODE = $ MXNODE =RTHICK,THICK(1),1,2,3,4!注意:根據(jù)不同的實(shí)際情況,以上命令流

15、只需要修改ABC部分內(nèi)容.!劃分好網(wǎng)格后,有一個(gè)很重要的事情：重新調(diào)整單元坐標(biāo)系,因?yàn)楹枚鄦卧膯卧鴺?biāo)系是亂七八糟的,但是將來結(jié)果文件輸出值的方向是單元坐標(biāo)系的.因此必須把所有單元的單元坐標(biāo)系調(diào)整的和整體坐標(biāo)系一樣(只是方向一樣,不一定x、y、z對(duì)應(yīng),因?yàn)闅んw的z軸永遠(yuǎn)是殼體的法向).命令如下：Asel,s,1 !選中1號(hào)面Esla,s,1 !選中屬于1號(hào)面的所有單元Local,11,0,1,0,0 !11局部坐標(biāo)系編號(hào);0表示笛卡爾坐標(biāo)系;1表示x值Emodif,all,esys,11 !修改單元坐標(biāo)系為局部坐標(biāo)系11注意:修改單元坐標(biāo)系的時(shí)候,必須先定義一個(gè)局部坐標(biāo)系.因?yàn)閑modif命

16、令只能利用局部坐標(biāo)系修改單元坐標(biāo)系.!將殼體法向反方向：有時(shí)候需要調(diào)整殼體法向的時(shí)候,可以利用如下命令：Asel,s,1AREVERSE,ALL !將面的方向反過來,ansys會(huì)提示你：是否也要將屬于面的單元法向也反過來?AllsFINISH!加載加約束/SOLUACEL,9.8 !施加重力加速度,詳細(xì)解釋見上.等大小的均布荷載 對(duì)殼體施加等大小的均布荷載比較簡(jiǎn)單,以下命令均可以實(shí)現(xiàn),只要你選中相應(yīng)的面、點(diǎn)或者單元.Sfa,ALL,1,PRES,100 Sf,ALL,PRES,100 !個(gè)人不推薦此命令,因?yàn)榇嗣顚?duì)施加在實(shí)體面上的節(jié)點(diǎn)有限制.Sfe,ALL,1,PRES,0,100 !主力推

17、薦此命令.1表示作用殼體的哪個(gè)面（如下圖所示 .）;0表示確實(shí)存在的壓力（SURF153 or SURF154單元選 用2）;100表示壓力(kn/m2).梯形荷載SFGRAD,PRES,0,Z,5,-10*1.2 !為下面的加載形成梯度。pres表示施加壓力；0表示 坐標(biāo)系代號(hào)(笛卡爾坐標(biāo)系）；z表示梯度沿z方向；5 表示：在0坐標(biāo)系下z5的位置(下面命令的300就加 在這個(gè)位置）；-10*1.2表示施加的梯度，負(fù)號(hào)表示隨 著z的增大，壓力在減小.SFe,ALL,2,PRES,0,300 !見上sfe命令解釋.!后處理（如下圖）!殼體單元輸出的是每個(gè)單元在單元坐標(biāo)系下，單位長(zhǎng)度的值。(這句話

18、是ansys幫助里翻譯而來的）有點(diǎn)含糊哈，大家湊合理解。我舉個(gè)列子給你們聽聽。如下例題： 計(jì)算圖 選取單元圖求基礎(chǔ)和板之間的剪力？用屁股想都想的出來，精確解是100kn。如此簡(jiǎn)單的問題，我們選用shell63做計(jì)算（當(dāng)然我們也可以選用平面單元，這里只是為了試驗(yàn)shell63求解結(jié)果的含義）對(duì)這個(gè)模型，我們分兩種情況計(jì)算，兩種情況唯一不同之處就是網(wǎng)格大小，第一種情況網(wǎng)格大小1*1,第二種情況網(wǎng)格大小0.5*0.5。取出靠近地面的最近的一排單元（如上圖），利用單元table輸出每個(gè)單元的剪力值（單元坐標(biāo)系下），然后相加。第一種情況，相加后A99.8 kn第二種情況，相加后A198.9 kn由上可得

19、:剪力值A(chǔ)*單元大小這里可以看出，劃分單元盡量規(guī)則、均勻的重要性，特別是將來你要提取內(nèi)力的位置，在劃分網(wǎng)格的時(shí)候就要做好準(zhǔn)備！而且計(jì)算前單元坐標(biāo)系要統(tǒng)一也非常重要！本人編寫了一個(gè)提取內(nèi)力的通用程序，各位可以參考參考。只要選中你要提取內(nèi)力的那排單元（如同上圖的那排），然后輸入以下命令即可：JLzonghe=0 !剪力ZLZONGHEx=0 !X方向軸力ZLZONGHEy=0 !Y方向軸力mzonghex=0 !X方向彎矩mzonghey=0 !Y方向彎矩etab,JL,smisc,3ETAB,ZL2,SMISC,2ETAB,ZL1,SMISC,1etab,m1,smisc,4etab,m2,sm

20、isc,5*GET,MAXN,ELEM,0,NUM,MAX !所選單元的最大編號(hào)給maxn*GET,MINN,ELEM,0,NUM,MIN*DO,I,MINN,MAXN *IF,ESEL(I),EQ,1,THEN !判斷單元號(hào)I是否屬于所選范圍 *get,jianli,elem,i,etab,JL *get,ZHOULIx,elem,i,etab,ZL1 *get,ZHOULIy,elem,i,etab,ZL2 *get,wanjux,elem,i,etab,m1 *get,wanjuy,elem,i,etab,m2 JLzonghe=JLzonghe+jianli ZLzonghex=ZLz

21、onghex+ZHOULIx ZLzonghey=ZLzonghey+ZHOULIy mzonghex=mzonghex+wanjux mzonghey=mzonghey+wanjuy *ENDIF*ENDDO注意：（1）以上命令求解出來的 Jlzonghe 等等還沒有乘以單元大小。（2）以上求解出來的X方向、Y方向的軸力和彎矩都是單元坐標(biāo)系下的值（具體方向如下圖,SX、SY就是單元坐標(biāo)系），單元坐標(biāo)系的X、Y方向可以通過視圖顯示出來，命令如下：/psymb,esys,1藍(lán)線表示Z軸，白線表示X軸，青線表示Y軸（3） Solid單元實(shí)體單元其實(shí)是ansys最簡(jiǎn)單的單元，但是也是ansys最難操

22、作的單元。實(shí)體單元從建模到形成計(jì)算書，70的時(shí)間用在劃分網(wǎng)格上，20的時(shí)間用在后處理上，還有10就是一些小雜事了。前期操作方面，有兩種方法可以完成：（1）利用ansys建模同時(shí)劃分網(wǎng)格。我有一個(gè)很牛的同事，做窩殼方面很有研究，他曾經(jīng)做過一個(gè)工程的窩殼，寫了整整2萬行的命令流，那窩殼真是漂亮，不得不讓人佩服。不過這種方法需要非常非常熟悉ansys，不是一般人可以的哦（汗顏！我為什么不行，55555）。（2）利用三維軟件建模(個(gè)人推薦solidwork,和ansys的兼容性非常ok)，然后導(dǎo)入ansys。這種方法比較通用，但是一般solidwork導(dǎo)入的都是一個(gè)整體塊，需要人工不停的用wp（工作面

23、）切割，切到可以映射劃分為止。這個(gè)切割技術(shù)還是比較難的，需要自己多去試試。一定要自己動(dòng)手去切割、去劃分，只要認(rèn)真做兩個(gè)例題，你就能自學(xué)成才了。上面廢話一堆，實(shí)體單元的前處理我就不說了，重點(diǎn)說一下我們土建專業(yè)比較關(guān)注的后處理：內(nèi)力（剪力、軸力、彎矩）提取。由于solid單元最后輸出的都是應(yīng)力，無法在工程中利用（當(dāng)然也可以根據(jù)規(guī)范提到的應(yīng)力配筋，但個(gè)人感覺不太穩(wěn)妥，我曾經(jīng)比較過：同一截面，應(yīng)力配筋和提取內(nèi)力的結(jié)果，內(nèi)力配筋結(jié)果都大于應(yīng)力配筋），但是我們可以利用這些應(yīng)力積分得到我們所需要的內(nèi)力。積分的原理就不說了，哪個(gè)兄弟要是忘了就去再看看結(jié)構(gòu)力學(xué)書吧。首先將工作平面的xy面移到你要計(jì)算內(nèi)力的位置，

24、然后你要判斷：（1） 你準(zhǔn)備計(jì)算內(nèi)力的截面屬于整體坐標(biāo)系的xy平面SUCR,SUZ2,CPLANE,3 !移動(dòng)工作平面、創(chuàng)建面SUZ2 SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y !映射X和Y方向應(yīng)力 SUMAP,MYSZ,S,Z$SUMAP,MYSXY,S,XY !映射Z和XY方向應(yīng)力 SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ !映射YZ和XZ方向應(yīng)力 SUEVAL,XFORCE,MYSXZ,INTG !求截面上FXSUEVAL,YFORCE,MYSYZ,INTG !求截面上FYSUEVAL,ZFORCE,MYSZ,INTG !求截面上FZSUEVA

25、L,MYA,DA,SUM !求截面面積并賦給MYA變量 SUEVAL,MYYA,GCY,INTG !求關(guān)于X軸的面積矩并賦給變量MYYA MYYA=MYYA/MYA !得到面積重心到X軸的距離=面積矩/面積 SUEVAL,MYXA,GCX,INTG !求關(guān)于Y軸的面積矩并賦給變量MYXA MYXA=MYXA/MYA !得到面積重心到Y(jié)軸的距離=面積矩/面積 SUCALC,SZGCY,MYSZ,MULT,GCY !計(jì)算MYSZ×GCY，并賦給面項(xiàng)SZGCY SUEVAL,MX1,SZGCY,INTG !對(duì)面項(xiàng)SZGCY在面上積分得到MX1 SUCALC,SZGCX,MYSZ,MULT,

26、GCX !計(jì)算MYSZ×GCZ，并賦給面項(xiàng)SZGCX SUEVAL,MY1,SZGCX,INTG !對(duì)面項(xiàng)SZGCX在面上積分得到MY1 !上述彎矩基于總體直角坐標(biāo)系原點(diǎn)而言的，應(yīng)對(duì)面積重心取矩，將內(nèi)力簡(jiǎn)化到面積重心上 MX1=MX1-ZFORCE*MYYA !MY1=MY1-ZFORCE*MYXA !（2） 你準(zhǔn)備計(jì)算內(nèi)力的截面屬于整體坐標(biāo)系的yz平面SUCR,SUZ2,CPLANE,3 !移動(dòng)工作平面、創(chuàng)建面SUZ2 SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y !映射X和Y方向應(yīng)力 SUMAP,MYSZ,S,Z$SUMAP,MYSXY,S,XY !映射Z和XY方向

27、應(yīng)力 SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ !映射YZ和XZ方向應(yīng)力 SUEVAL,XFORCE,MYSx,INTG !求截面上FXSUEVAL,zFORCE,MYSXZ,INTG !求截面上FzSUEVAL,yFORCE,MYSXY,INTG !求截面上FySUEVAL,MYA,DA,SUM !求截面面積并賦給MYA變量 SUEVAL,MYZA,GCZ,INTG !求關(guān)于X軸的面積矩并賦給變量MYYA MYZA=MYZA/MYA !得到面積重心到Y(jié)軸的距離=面積矩/面積 SUEVAL,MYYA,GCY,INTG !求關(guān)于Z軸的面積矩并賦給變量MYYA MYYA=

28、MYYA/MYA !得到面積重心到Z軸的距離=面積矩/面積 SUCALC,SZGCZ,MYSX,MULT,GCZ !計(jì)算MYSX×GCZ，并賦給面項(xiàng)SZGCZ SUEVAL,MY1,SZGCZ,INTG !對(duì)面項(xiàng)SZGCZ在面上積分得到MY1 SUCALC,SZGCY,MYSX,MULT,GCY !計(jì)算MYSX×GCY，并賦給面項(xiàng)SZGCY SUEVAL,MZ1,SZGCY,INTG !對(duì)面項(xiàng)SZGCY在面上積分得到MZ1 !上述彎矩基于總體直角坐標(biāo)系原點(diǎn)而言的，應(yīng)對(duì)面積重心取矩，將內(nèi)力簡(jiǎn)化到面積重心上 MY1=MY1-XFORCE*MYZA !MZ1=MZ1-XFORCE

29、*MYYA !（3） 你準(zhǔn)備計(jì)算內(nèi)力的截面屬于整體坐標(biāo)系的xz平面SUCR,SUZ2,CPLANE,3 !移動(dòng)工作平面、創(chuàng)建面SUZ2 SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y !映射X和Y方向應(yīng)力 SUMAP,MYSZ,S,Z$SUMAP,MYSXY,S,XY !映射Z和XY方向應(yīng)力 SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ !映射YZ和XZ方向應(yīng)力 SUEVAL,XFORCE,MYSxy,INTG !求截面上FXSUEVAL,zFORCE,MYSYZ,INTG !求截面上FzSUEVAL,yFORCE,MYSy,INTG !求截面上FySUEVA

30、L,MYA,DA,SUM !求截面面積并賦給MYA變量 SUEVAL,MYZA,GCZ,INTG !求關(guān)于X軸的面積矩并賦給變量MYYA MYZA=MYZA/MYA !得到面積重心到X軸的距離=面積矩/面積 SUEVAL,MYXA,GCX,INTG !求關(guān)于Z軸的面積矩并賦給變量MYXA MYXA=MYXA/MYA !得到面積重心到Z軸的距離=面積矩/面積 SUCALC,SZGCZ,MYSy,MULT,GCZ !計(jì)算MYSY×GCZ，并賦給面項(xiàng)SZGCZ SUEVAL,MX1,SZGCZ,INTG !對(duì)面項(xiàng)SZGCZ在面上積分得到MX1 SUCALC,SZGCX,MYSy,MULT,

31、GCX !計(jì)算MYSY×GCX，并賦給面項(xiàng)SZGCX SUEVAL,MZ1,SZGCX,INTG !對(duì)面項(xiàng)SZGCX在面上積分得到MY1 !上述彎矩基于總體直角坐標(biāo)系原點(diǎn)而言的，應(yīng)對(duì)面積重心取矩，將內(nèi)力簡(jiǎn)化到面積重心上 MX1=MX1-yFORCE*MYZA !MZ1=MZ1-yFORCE*MYXA !注意：（1） XFORCE、YFORCE、ZFORCE、MX1、MY1、MZ1為最后的計(jì)算結(jié)果，而且是在整體坐標(biāo)系下的結(jié)果。（2） 對(duì)于xy平面，ZFORCE為該斷面的軸力；XFORCE、YFORCE分別表示x、y方向的剪力；MX1、MY1分別表示右手準(zhǔn)則，大拇指分別指向整體坐標(biāo)系下的

32、x、y軸。其他yz、xz平面類似。（3） 這里強(qiáng)調(diào)一下求得內(nèi)力的方向問題。我們求解內(nèi)力的時(shí)候，其實(shí)就是對(duì)各應(yīng)力進(jìn)行面積積分，所得內(nèi)力正負(fù)值的意義和各應(yīng)力一樣。在ansys中各應(yīng)力的正負(fù)的意義是和有限元定義正負(fù)號(hào)的方法是一致的，和單純的整體坐標(biāo)系的正負(fù)號(hào)是不一樣的。舉例如下： 比如你求一個(gè)屬于xy平面的截面的內(nèi)力，求出剪力100kn。那么對(duì)于這個(gè)剪力方向是什么呢？根據(jù)有限元正負(fù)號(hào)定義原則！可參照有限元原理母校教材哈，徐之倫寫的經(jīng)典的。實(shí)在不想看，就q我，我qq：。（4） 以前我一個(gè)同事發(fā)現(xiàn)了一個(gè)命令，很簡(jiǎn)單的一個(gè)命令輸入，然后就可以輸出以上結(jié)果。但是對(duì)于此命令得到結(jié)果，本人不太明白他輸出的意思，