有限元計算注意事項

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

2、復 雜而且龐大時，隱式求解方法的計算量還是相當大的.與之相對應的，顯式求解方法.顯式求解方法是通過動態(tài)方法從 一個增量步前推到下一個增量步得到的.具體顯式求解方法和隱式求 解方法例子如下：（1）隱式求解圖1-3桁架間題圖1-4桁架問題的離散化模型圖1 5顯示了模型中每個節(jié)點的分離圖，在通常情況下，模型中的每個節(jié)點將承受 外部載荷P和內(nèi)部載荷L后者是由于節(jié)點相連單元的應力引起的。對于一個處于靜平 衡狀態(tài)的模型，每個節(jié)點上的合力必須為零，即每巳 : 個節(jié)點上的外部載荷和內(nèi)部載荷必須相互平，衡 一- 節(jié)占口 0 與對于節(jié)點八其平衡方程可建立如下.= j2 r假設圓桿的長度變化很小，則單元1的應變節(jié)點

3、3 L可表示為,圖1-5每個節(jié)點的分離圖u U11 = 7其中,必利分別是節(jié)點莊和&的位移工是單元的初始長度。假設材料是彈性的，圓桿的應力可以通過應變乘以楊氏模量E紿出:少1 Eeu作用在端部節(jié)點的軸向力等于圓桿的應力乘以其橫截面積A.這樣,內(nèi)力、材料性能 和位移的關系為7： = mA = Feu A = ua）在節(jié)點a處的平衡方程因此可以寫成巳+華（把一/） = 0節(jié)點&的平衡必須同時考慮與該節(jié)點相連的兩個單元的內(nèi)力.對于節(jié)點鼠單元1的 內(nèi)力作用在反方向上，因此變?yōu)樨撝怠Ｆ胶夥匠虨镻b 錚（盤一黃）+ 詈（if 一/） = 0節(jié)點的平衡方程為巳一譬3-）=0對于隱式方法，這些平衡方程需要同時

4、進行求解以獲取每個節(jié)點的位移，以上求解 最好采用矩陣算法，因此，將內(nèi)力和外力的貢獻寫成矩陣形式，如果兩個單元的性質(zhì)和 維數(shù)相同，平衡方程可以化簡成如下形式：0通常情況下，每個單元的剛度(本例中的EA/L項)是不同的*因此，可以將模型中 兩個單元的剛度分別寫成和我們感興趣的是獲得平衡方程的解答，使施加的外 力P與產(chǎn)生的內(nèi)力I達到平衡若考慮到問題的收斂性和非線性，可以將平衡方程寫成 一洱=o對于完整的兩單元三節(jié)點結構，經(jīng)過這種符號變換后，可以將平衡方程重寫為- M-kl-0在隱式方法中，例如應用在/XBAQUS/Standard中，由此方程組能夠解出3個未知變量的 值；護，討及在該問題中，/已指定

5、為0.0R 一旦位移求出后，就能利用位移返回計 算出桁架單元的應力當每一個求解增量步結束時，隱式的有限元方法要求解一組方 程組。顯式求解本節(jié)提供從概念上理解當應用顯式動力學算法時應力是如何在模型中傳播的在這 個演示例子中，考慮應力波沿著一個由3個單元構成的桿件模型傳播的過程，如圖1-6所 示，隨著時間增量的變化，我們將研究桿件的各個狀態(tài)。1 . 渤1-6自由端作用有集中力F的桿件的初始狀態(tài)在第1個時間增量段，施加在節(jié)點1的集中力尸使節(jié)點I產(chǎn)生了一個加速度成，這 個加速度使節(jié)點1產(chǎn)生速度W，在.單元內(nèi)引起應變速率鈕在第1個時間增量段內(nèi)對 應變速率進行積分就獲得了單元的應變增量&叫。總應變是初始應

6、變勤和應變增 量的和在這個問題中的初始應變?yōu)榱?一旦計算出了單元應變，就可以通過材料的本 構模型求出單元的應力如對于線彈性材料，應力是彈性模量與總應變的乘積。這個過 程如圖7所示。在第1個時間增量段里，節(jié)點2和節(jié)點3因為沒有力作用在其上，所以 沒有移動.在第2個時間增錄段，由單元的應力產(chǎn)生的單元內(nèi)力被施加到了與單元相連的所 有節(jié)點上，如圖8所示。然后應用這些單元應力計算節(jié)點1和節(jié)點2的動力平衡 方程，（D 234；p I I -IW/ /Hi =壽 = I Hi denial = =昏deg =的 +&=&cil圖1-7自由端作用有集中力的桿件在第1時間增履段結束時的旅態(tài) TOC o 1-5

7、h z 1W 34，P - 4* K-M S-4ii=hr. F Lllf- old. _L （-M2 M1 Jr 街.fl = H +“由 陽-七二=理小=jEaJ r=fU =E1 +血1緝H金寸臨=J*由FBg圈1-8桿件在第2個時間增量段開始時的狀態(tài)這個過程繼續(xù)下去，到第3個時間增量段開始時，單元和單元中都已存在了應 力，而節(jié)點1,2和3則存在了作用力，如圖1.-9所示。這個過程將-一直繼續(xù)下去，直到1總 的分析時間結束- TOC o 1-5 h z | 2 1 A；p - I y扇1 412F圖9桿件在第3個時間增量段開始時的狀態(tài)隱式求解中，計算的精度完全控制于計算步數(shù)，在一般的計算

8、軟 件中（flac、abaqus），軟件均是利用不平衡力來控制計算步數(shù)（當 不平衡力10-5時，停止計算）.不平衡力=A+B.A表示施加在節(jié)點上的 集中力；B表示：在n步數(shù)下，根據(jù)第n步計算出來的應力，求出節(jié)點 的內(nèi)力.pby mFlac軟件中B 4 3，以上公式是根據(jù)虛功原理推倒而得到.具體推倒過程見flac原理.2、Ansys計算注意事項：計算單位、參數(shù)、荷載、標準值、設計值，計算過程中系數(shù)的加入.（1） beam 單元對于beam單元.Ansys軟件中我們常用的有兩種梁單元:beam188和beam4.這兩種單元均是三維的梁單元，每個節(jié)點都具有6個自由度（ux、uy、uz、mx、my、m

9、z），并且單元坐標系x軸是i點指向j點.Beam188單元是基于鐵木辛哥理論的梁,beam4單元是我們常用的經(jīng)典結構力學梁.（鐵木辛哥理論考慮了梁的剪切變形,而我們常用的經(jīng)典結構力學梁只考慮了彎矩對結構的變形影響）所以說,beam188可以更精確的計算梁單元,因此我們結構計算中，一般都采用beam188單元.當然還有beam189單元，189單元屬于三維二次的梁單元（beam188屬于三維一次梁單元），精度比beam188更加高.定義beam188單元,一般采用如下形式：!定義單元!進入前處理!定義單元188號標號為1/prep7 et,1,beam188!定義材料屬性mp,ex,1,2.55

10、e7!定義彈性模量（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型等 等,具體見ansys幫助）；A1表示梁的名稱；0表示 薄壁梁單元網(wǎng)格劃分精細程度（05）.SECDATA,1,3,4,12!1表示梁b ; 3表示梁h ; 4和12定義對應寬長等分份數(shù).SECOFFSET,CENT!cent 質(zhì)心；shrc 剪切中心；origin 原始中心；use

11、r用戶定義；!注意：當梁單元和殼單元一起使用時，可以設置梁單元的偏心，使梁的一面和殼的一面共面.(secoffset,user,offsety,offset z),如下圖：W2!劃分網(wǎng)格LSEL,S,1LATT,1,1,1,1LESIZE,ALL,0.2,1LMESH,ALL!加載加約束/SOLU ACEL,9.8SFCUM,ALL,ADD!選中編號為1的線.!mp,r,et,，方向點，SECTYPE 截面號.!0.2是單元大小，1是確認細分規(guī)則.!用beam單元離散模型，形成網(wǎng)格.!對于劃分網(wǎng)格，空間的beam單元，由于需要確定b、h的方向，ansys軟件利用方向點來控 制b、h的方向.方向

12、點的編號最好定義的很大，如果定義太小,會影響后面的加載.具體方向點如何控制見上面的latt命令和ansys幫助.自己試兩下就知道怎么用了.AllsFINISH!重力加速度.注意方向，數(shù)值和整體坐標相反,比如重力指向z軸 負向，則為正值.!設置單元荷載是疊加還是替代，只對加在單元和節(jié)點上的荷載有 效，對于加在面、線上的荷載，都只有替代作用(對同一個面，第 二次加的荷載替代第一次加的荷載)!對于beam單元，只能根據(jù)sfbeam命令增加均布荷載等大小的均布荷載.Lsel,s,1ESLL,S,1sfbeam,ALL,1,PRES,-161.5!1表示作用在beam單元的面上(如下圖，面表示beam單

13、元的軸向，面表示單元側面，面表示beam單元頂 面)，-161.5表示均布荷載大小，正負號可以控制作用力的 方向.梯形均布荷載Sfbeam命令是對每個單元進行加載.如果一根梁承受10100的梯形均布荷載，而且這根梁 被分成了 10個beam單元，這樣施加荷載就非常困難.因此我將這種加載過程寫成命令流，讓 軟件自動進行加載.命令流如下：!選中要加載的那根梁(線)LSEL,S,1ESLL,S,1!選中屬于這根梁(線)的beam單元*GET,Nelem,ELEM,COUNT, , , ,!獲得當前所選單元個數(shù)，賦予參數(shù)Nelem*GET,Ne,ELEM,NUM,MIN, , , ,!獲得當前所選單元

14、最小編號，賦予參數(shù)Ne*DO,I,1,Nelem!循環(huán)加載，循環(huán)次數(shù)=單元個數(shù)ESEL,S,NeNSLE,S,1*GET,Nnode,NODE,COUNT, , , , !獲得當前所選節(jié)點個數(shù)，賦予參數(shù)Nnode*GET,Nn,NODE,NUM,MIN,NN1X=NX(NN)NN=NDNEXT(NN)NN2X=NX(NN)sfbeam,ALL,1,PRES,!獲得當前所選節(jié)點最小編號，賦予參數(shù)Nn!將nn節(jié)點的x坐標賦予NN1X(NX表示x坐標,NY表示y坐標)!NN=當前所選節(jié)點的下一個編號!將nn節(jié)點的x坐標賦予NN2X-1630.76/3.23*(3.23-NN1X),-1630.76/

15、3.23*(3.23-NN2X!以上荷載公式應根據(jù)實際情況進行調(diào)整LSEL,S,1ESLL,S,1NE=ELnext(NE)*ENDDO!對于此命令流，根據(jù)不同的實際情況,!NE =當前所選單元的下一個編號部分需要修改，其他不需要修改.!后處理（XY平面）（大拇指指向y,就是my）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,

16、JFY!注意：beam單元的結果輸出都是以單元坐標系輸出的，且拉為正、壓為負前面我們已經(jīng)知 道，單元坐標系x軸就是i點指向j點,其他坐標可以根據(jù)整體坐標系推出.詳細內(nèi)容見ansys 幫助.（2）Shell 單元對于shell單元應用的范圍,ansys軟件并沒有強制規(guī)定，只是從 字面上區(qū)分了薄殼和厚殼.我以前看過一本電子教案仿真在線， 里面說一般規(guī)定殼體的主尺寸是厚度的10倍左右，都是可以用殼體 來模擬的.一般高度與跨度之比（非與單元尺寸比較）1/15 & 1/10,可以當作厚殼來處理.shell63是薄殼單元，他包含彎曲和薄膜效應,但是忽略橫向剪切變 形;shell43,shell143,sh

17、ell181,shell91,shell93 和 shell99，都屬于厚殼單元, 不僅有彎曲、薄膜效應，他也包含了橫向剪切效應.橫向剪切被表示為 整個厚度上的常剪切應變.這種一階近似只適用于中等厚度殼體.線形分析時，如果不包含橫向剪切應變，使用63,163單元;如果橫 向剪切變形重要,則遵守以下原則:均勻材料，使用43,93,143單元，復合 材料使用91,99,181.我們土木工程中，一般利用shell43計算.!定義單元!進入前處理!定義單元43號標號為1!定義彈性模量(kn/m2)!定義泊松比!定義密度(KN/N*KG/M3)!1表示編號,2表示厚度(m)/prep7 et,1,she

18、ll43 !定義材料屬性 mp,ex,1,2.55e7 mp,nuxy,1,0.167 mp,dens,1,2.5 !定義墻體厚度!等厚度板R,1,2R,2,3!選中1號面!mp,real,type!定義單元大小為0.2左右!規(guī)定劃分單元形狀,0表示四邊形(1表示三角形),2d表示劃分面(3d 表示劃分體)!指定是自由劃分還是映射劃分,2表示:盡量用映射劃分，不符合要求 就自動使用自由劃分，具體參見ansys幫助的eshkey命令.!劃分面單元.Asel,s,1Aatt,1,1,1ESIZE,0.2MSHKEY,2MSHAPE,0,2DAMESH,ALL！注意：在網(wǎng)格剖面方面，最好全部用四邊形

19、，而且形狀盡量規(guī)則、均勻！因為將來后處 理內(nèi)力提取的時候，提取出來的力和單元的大小有直接的關系。后處理部分將詳細講解。漸變厚度的板對于漸變的板，首先選中要劃分的面，然后進行網(wǎng)格劃分.但是在網(wǎng)格劃分的 時候，賦予的厚度編號開始選用100 (前面沒有定義的就可以)，然后利用以下命 令進行厚度賦值.ASEL,S,1ESLAMXNODE = NDINQR(0,14)*DIM,THICK,MXNODE*DO,NODE,1,MXNODE3-(NZ(NODE)-7.5)/25.5*1.5 !厚度和點坐標系的關系,NZ(NODE)表示 編號為node的點的Z坐標值.*IF,NSEL(NODE),EQ,1,TH

20、ENTHICK(node)= *ENDIF*ENDDONODE = $ MXNODE =RTHICK,THICK(1),1,2,3,4!注意:根據(jù)不同的實際情況,以上命令流只需要修改ABC部分內(nèi)容.x、y、zAsel,s,1Esla,s,1Local,11,0,1,0,0Emodif,all,esys,11!劃分好網(wǎng)格后，有一個很重要的事情：重新調(diào)整單元坐標系，因為好多單元的單 元坐標系是亂七八糟的,但是將來結果文件輸出值的方向是單元坐標系的.因此 必須把所有單元的單元坐標系調(diào)整的和整體坐標系一樣（只是方向一樣,不一定對應，因為殼體的z軸永遠是殼體的法向）.命令如下：!選中1號面!選中屬于1號

21、面的所有單元!11局部坐標系編號;0表示笛卡爾坐標系;1表示x值!修改單元坐標系為局部坐標系11!將殼體法向反方向：有時候需要調(diào)整殼體法向的時候,可以利用如下命令:Asel,s,1AREVERSE,ALL注意:修改單元坐標系的時候,必須先定義一個局部坐標系.因為emodif命令只能利用局部 坐標系修改單元坐標系.!將面的方向反過來,ansys會提示你:是否也要將屬于面的單元法向 也反過來？AllsFINISH!加載加約束/SOLU!施加重力加速度，詳細解釋見上.ACEL,9.8等大小的均布荷載對殼體施加等大小的均布荷載比較簡單,以下命令均可以實現(xiàn)，只要你選中相應的面、點 或者單元.Sfa,AL

22、L,1,PRES,100Sf,ALL,PRES,100!個人不推薦此命令，因為此命令對施加在實體面上的節(jié)點有限制.Sfe,ALL,1,PRES,0,100!主力推薦此命令.1表示作用殼體的哪個面（如下圖所示.）;0表示確實存在的壓力（SURF153 or SURF154單元選 用2） ;100表示壓力（kn/m2）.梯形荷載SFGRAD,PRES,0,Z,5,-10*1.2!為下面的加載形成梯度。pres表示施加壓力；0表示坐標系代號（笛卡爾坐標系）；z表示梯度沿z方向；5 表示:在0坐標系下z = 5的位置（下面命令的300就加 在這個位置）；-10*1.2表示施加的梯度，負號表示隨 著z的

23、增大，壓力在減小.SFe,ALL,2,PRES,0,300!見上 sfe 命令解釋.!后處理（如下圖）!殼體單元輸出的是每個單元在單元坐標系下，單位長度的值。（這句話是 ansys幫助里翻譯而來的）有點含糊哈，大家湊合理解。我舉個列子給你們聽聽。 如下例題：PI I計算圖選取單元圖求基礎和板之間的剪力？用屁股想都想的出來，精確解是100kn。如此簡單的問題，我們選用shell63做計算（當然我們也可以選用平面單元， 這里只是為了試驗shell63求解結果的含義）對這個模型，我們分兩種情況計算，兩種情況唯一不同之處就是網(wǎng)格大小， 第一種情況網(wǎng)格大小1*1,第二種情況網(wǎng)格大小0.5*0.5。取出靠

24、近地面的最近的一排單元（如上圖），利用單元table輸出每個單元 的剪力值（單元坐標系下），然后相加。第一種情況，相加后A=99.8 kn第二種情況，相加后A= 198.9 kn由上可得:剪力值=A*單元大小這里可以看出，劃分單元盡量規(guī)則、均勻的重要性，特別是將來你要提取內(nèi)力的位置， 在劃分網(wǎng)格的時候就要做好準備！而且計算前單元坐標系要統(tǒng)一也非常重要!！ 本人編寫了一個提取內(nèi)力的通用程序，各位可以參考參考。只要選中你要提取內(nèi)力的那排單元（如同上圖的那排），然后輸入以下命令即可:JLzonghe=0!剪力ZLZONGHEx=0!X方向軸力ZLZONGHEy=0!Y方向軸力mzonghex=0!X

25、方向彎矩mzonghey=0!Y方向彎矩etab,JL,smisc,3ETAB,ZL2,SMISC,2ETAB,ZL1,SMISC,1etab,m1,smisc,4etab,m2,smisc,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,wa

26、njux,elem,i,etab,m1*get,wanjuy,elem,i,etab,m2JLzonghe二JLzonghe+jianliZLzonghex二ZLzonghex+ZHOULIxZLzonghey二ZLzonghey+ZHOULIymzonghex二mzonghex+wanjuxmzonghey二mzonghey+wanjuy*ENDIF*ENDDO注意：(1)以上命令求解出來的Jlzonghe等等還沒有乘以單元大小。(2)以上求解出來的X方向、Y方向的軸力和彎矩都是單元坐標系下的值(具體方向 如下圖,SX、SY就是單元坐標系)，單元坐標系的X、Y方向可以通過視圖顯示出來，命令如

27、 下：/psymb,esys,1藍線表示Z軸，白線表示X軸，青線表示Y軸,翊TQP3SX表j = Element x：-axis if ESY；is not supplied.k = Element if ESYS is supplied.Solid 單元實體單元其實是ansys最簡單的單元，但是也是ansys最難操作的單元。實體單元從建模到形成計算書，70%的時間用在劃分網(wǎng)格上，20%的時間用 在后處理上，還有10%就是一些小雜事了。前期操作方面，有兩種方法可以完成：（1）利用ansys建模同時劃分網(wǎng)格。我有一個很牛的同事，做窩殼方面很有研究，他曾經(jīng)做過一個工程的窩殼， 寫了整整2萬行的命令

28、流，那窩殼真是漂亮，不得不讓人佩服。不過這種方法需 要非常非常熟悉ansys，不是一般人可以的哦（汗顏！我為什么不行，55555）。（2）利用三維軟件建模（個人推薦solidwork,和ansys的兼容性非常ok）， 然后導入ansys。這種方法比較通用，但是一般solidwork導入的都是一個整體塊，需要人工 不停的用wp （工作面）切割，切到可以映射劃分為止。這個切割技術還是比較 難的，需要自己多去試試。一定要自己動手去切割、去劃分，只要認真做兩個例 題，你就能自學成才了。上面廢話一堆，實體單元的前處理我就不說了，重點說一下我們土建專業(yè)比 較關注的后處理：內(nèi)力（剪力、軸力、彎矩）提取。由于

29、solid單元最后輸出的都是應力，無法在工程中利用（當然也可以根據(jù)規(guī) 范提到的應力配筋，但個人感覺不太穩(wěn)妥，我曾經(jīng)比較過：同一截面，應力配筋和提取內(nèi)力 的結果，內(nèi)力配筋結果都大于應力配筋），但是我們可以利用這些應力積分得到我們 所需要的內(nèi)力。積分的原理就不說了，哪個兄弟要是忘了就去再看看結構力學書 吧。首先將工作平面的xy面移到你要計算內(nèi)力的位置，然后你要判斷:（1）你準備計算內(nèi)力的截面屬于整體坐標系的xy平面SUCR,SUZ2,CPLANE,3!移動工作平面、創(chuàng)建面SUZ2SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y !映射 X 和 Y 方向應力 SUMAP,MYSZ,S,Z$

30、SUMAP,MYSXY,S,XY !映射 Z 和 XY 方向應力 SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ !映射 YZ 和 XZ 方向應力 SUEVAL,XFORCE,MYSXZ,INTG !求截面上 FX SUEVAL,YFORCE,MYSYZ,INTG !求截面上 FY SUEVAL,ZFORCE,MYSZ,INTG!求截面上 FZSUEVAL,MYA,DA,SUM!求截面面積并賦給MYA變量SUEVAL,MYYA,GCY,INTG !求關于X軸的面積矩并賦給變量MYYA MYYA=MYYA/MYA!得到面積重心到X軸的距離二面積矩/面積SUEVAL,MYXA,

31、GCX,INTG!求關于Y軸的面積矩并賦給變量MYXAMYXA=MYXA/MYA!得到面積重心到Y軸的距離二面積矩/面積SUCALC,SZGCY,MYSZ,MULT,GCY !計算 MYSZ X GCY,并賦給面項 SZGCY SUEVAL,MX1,SZGCY,INTG !對面項SZGCY在面上積分得到MX1 SUCALC,SZGCX,MYSZ,MULT,GCX !計算 MYSZXGCZ,并賦給面項 SZGCX SUEVAL,MY1,SZGCX,INTG!對面項SZGCX在面上積分得到MY1!上述彎矩基于總體直角坐標系原點而言的，應對面積重心取矩，將內(nèi)力簡化到面積重心上 MX1=MX1-ZFO

32、RCE*MYYA!MY1=MY1-ZFORCE*MYXA!（2）你準備計算內(nèi)力的截面屬于整體坐標系的yz平面SUCR,SUZ2,CPLANE,3!移動工作平面、創(chuàng)建面SUZ2SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y !映射 X 和 Y 方向應力 SUMAP,MYSZ,S,Z$SUMAP,MYSXY,S,XY !映射 Z 和 XY 方向應力 SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ !映射 YZ 和 XZ 方向應力 SUEVAL,XFORCE,MYSx,INTG !求截面上 FX SUEVAL,zFORCE,MYSXZ,INTG !求截面上 Fz SU

33、EVAL,yFORCE,MYSXY,INTG!求截面上 FySUEVAL,MYA,DA,SUM!求截面面積并賦給MYA變量SUEVAL,MYZA,GCZ,INTG !求關于X軸的面積矩并賦給變量MYYA MYZA=MYZA/MYA!得到面積重心到Y軸的距離二面積矩/面積SUEVAL,MYYA,GCY,INTG!求關于Z軸的面積矩并賦給變量MYYAMYYA=MYYA/MYA!得到面積重心到Z軸的距離二面積矩/面積SUCALC,SZGCZ,MYSX,MULT,GCZ !計算 MYSXXGCZ,并賦給面項 SZGCZ SUEVAL,MY1,SZGCZ,INTG !對面項SZGCZ在面上積分得到MY1

34、 SUCALC,SZGCY,MYSX,MULT,GCY !計算 MYSX X GCY，并賦給面項 SZGCY SUEVAL,MZ1,SZGCY,INTG!對面項SZGCY在面上積分得到MZ1!上述彎矩基于總體直角坐標系原點而言的，應對面積重心取矩，將內(nèi)力簡化到面積重心上 MY1=MY1-XFORCE*MYZA!MZ1=MZ1-XFORCE*MYYA!你準備計算內(nèi)力的截面屬于整體坐標系的xz平面SUCR,SUZ2,CPLANE,3!移動工作平面、創(chuàng)建面SUZ2SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y !映射 X 和 Y 方向應力 SUMAP,MYSZ,S,Z$SUMAP,MYSXY,S,XY !映射 Z 和 XY 方向應力 SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ !映射 YZ 和 XZ 方向應力 SUEVAL,XFORCE,MYSxy,INTG !求截面上 FXSUEVAL,zFORCE,MYSYZ,INTG !求截面上 FzSUEVAL,y