有限元計算注意事項

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

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

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

4、,189單元屬于三維二次的梁單元（beam188屬于三維一次梁單元）,精度比beam188更加高.定義beam188單元,一般采用如下形式：!定義單元/prep7 !進入前處理et,1,beam188 !定義單元188號標(biā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é)點numcmp,all ！壓縮編號!定義梁截面SECTYPE,1,BEAM,RECT,A1,0 ! 1表示梁編號 ; RECT表示是矩形梁（還有其他t型等等,具

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

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

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

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

9、1X=NX(NN) !將nn節(jié)點的x坐標(biāo)賦予NN1X(NX表示x坐標(biāo),NY表示y坐標(biāo))NN=NDNEXT(NN) !NN=當(dāng)前所選節(jié)點的下一個編號NN2X=NX(NN) !將nn節(jié)點的x坐標(biāo)賦予NN2Xsfbeam,ALL,1,PRES,-1630.76/3.23*(3.23-NN1X),-1630.76/3.23*(3.23-NN2X)!以上荷載公式應(yīng)根據(jù)實際情況進行調(diào)整LSEL,S,1ESLL,S,1NE=ELnext(NE) !NE=當(dāng)前所選單元的下一個編號 *ENDDO!對于此命令流,根據(jù)不同的實際情況,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點指向j點,其他坐標(biāo)可以根據(jù)整體坐標(biāo)系推出.詳細(xì)內(nèi)容見ansys 幫助.（2） Shell單元對于shell單元應(yīng)用的范圍,ansys軟件并沒有強制規(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)橫向剪力所引起的變形和彎曲變形屬同一量級,中厚板應(yīng)該考慮.橫向剪切被表示為整個厚度上的常剪切應(yīng)變.

12、這種一階近似只適用于中等厚度殼體. 線形分析時,如果不包含橫向剪切應(yīng)變,使用63,163單元;如果橫向剪切變形重要,則遵守以下原則:均勻材料,使用43,93,143單元,復(fù)合材料使用91,99,181.我們土木工程中,一般利用shell43計算.!定義單元/prep7 !進入前處理et,1,shell43 !定義單元43號標(biā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表示編號,2表示厚度（m）R,2,3Asel,s,

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

14、后進行網(wǎng)格劃分.但是在網(wǎng)格劃分的時候,賦予的厚度編號開始選用100（前面沒有定義的就可以）,然后利用以下命令進行厚度賦值.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 !厚度和點坐標(biāo)系的關(guān)系,NZ(NODE)表示編號為node的點的Z坐標(biāo)值.*ENDIF*ENDDONODE = $ MXNODE =RTHICK,THICK(1),1,2,3,4!注意:根據(jù)不同的實際情況,以上命令流

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

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

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

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

19、:剪力值A(chǔ)*單元大小這里可以看出，劃分單元盡量規(guī)則、均勻的重要性，特別是將來你要提取內(nèi)力的位置，在劃分網(wǎng)格的時候就要做好準(zhǔn)備！而且計算前單元坐標(biāo)系要統(tǒng)一也非常重要！本人編寫了一個提取內(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 !所選單元的最大編號給maxn*GET,MINN,ELEM,0,NUM,MIN*DO,I,MINN,MAXN *IF,ESEL(I),EQ,1,THEN !判斷單元號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藍線表示Z軸，白線表示X軸，青線表示Y軸（3） Solid單元實體單元其實是ansys最簡單的單元，但是也是ansys最難操

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

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

24、然后你要判斷：（1） 你準(zhǔn)備計算內(nèi)力的截面屬于整體坐標(biāo)系的xy平面SUCR,SUZ2,CPLANE,3 !移動工作平面、創(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 !計算MYSZ×GCY，并賦給面項SZGCY SUEVAL,MX1,SZGCY,INTG !對面項SZGCY在面上積分得到MX1 SUCALC,SZGCX,MYSZ,MULT,

26、GCX !計算MYSZ×GCZ，并賦給面項SZGCX SUEVAL,MY1,SZGCX,INTG !對面項SZGCX在面上積分得到MY1 !上述彎矩基于總體直角坐標(biāo)系原點而言的，應(yīng)對面積重心取矩，將內(nèi)力簡化到面積重心上 MX1=MX1-ZFORCE*MYYA !MY1=MY1-ZFORCE*MYXA !（2） 你準(zhǔn)備計算內(nèi)力的截面屬于整體坐標(biāo)系的yz平面SUCR,SUZ2,CPLANE,3 !移動工作平面、創(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 !計算MYSX×GCZ，并賦給面項SZGCZ SUEVAL,MY1,SZGCZ,INTG !對面項SZGCZ在面上積分得到MY1 SUCALC,SZGCY,MYSX,MULT,GCY !計算MYSX×GCY，并賦給面項SZGCY SUEVAL,MZ1,SZGCY,INTG !對面項SZGCY在面上積分得到MZ1 !上述彎矩基于總體直角坐標(biāo)系原點而言的，應(yīng)對面積重心取矩，將內(nèi)力簡化到面積重心上 MY1=MY1-XFORCE*MYZA !MZ1=MZ1-XFORCE

29、*MYYA !（3） 你準(zhǔn)備計算內(nèi)力的截面屬于整體坐標(biāo)系的xz平面SUCR,SUZ2,CPLANE,3 !移動工作平面、創(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 !計算MYSY×GCZ，并賦給面項SZGCZ SUEVAL,MX1,SZGCZ,INTG !對面項SZGCZ在面上積分得到MX1 SUCALC,SZGCX,MYSy,MULT,

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

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