采用有限元的非線性屈曲分析就是要尋找上述過程中的A點時的壓力，有兩種方式： （1） 加強制位移約束，然后輸出反力，做出反力和位移的曲線圖，直接在圖上查看馬鞍點位置對應的反力大小。

現在Abaqus、LS-DYNA、Ansys等結構商軟都說可以處理復雜的上萬零部件接觸的整車、整機等模型仿真，沒做過實際的這種仿真分析，很好奇，接觸分析算法往往涉及大變形、邊界不連續，只要輸入條件或者算法稍微變化一些，兩個零部件算出來的接觸結果就可能差異很大，更不用說上萬個零部件的接觸結果了，對這種大規模組裝模型的仿真結果不知如何來判斷它的可靠性，像普通的只校核一下材料的應力還是看一下動畫是否和試驗一致

現在Abaqus、LS-DYNA、Ansys等結構商軟都說可以處理復雜的上萬零部件接觸的整車、整機等模型仿真，沒做過實際的這種仿真分析，很好奇，接觸分析算法往往涉及大變形、邊界不連續，只要輸入條件或者算法稍微變化一些，兩個零部件算出來的接觸結果就可能差異很大，更不用說上萬個零部件的接觸結果了，對這種大規模組裝模型的仿真結果不知如何來判斷它的可靠性，像普通的只校核一下材料的應力還是看一下動畫是否和試驗一致

我們關注CAE中的結構有限元，所以主要選擇了商用結構有限元軟件中文檔相對較完備的Abaqus來研究內部實現方式，同時對某些問題也會涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限可能有許多的理論錯誤，歡迎交流討論，也期待有更多的合作機會。

讀者對象為初步接觸求解器的軟件研發人員，文中提到的很多用顏色標記的內容可以作為關鍵字單獨搜索。

3）畢竟現在同屬于Ansys的模塊，相比EDEM，Rocky和Fluent的耦合更加絲滑。而EDEM與Fluent需要通過UDF耦合，并且有版本限制，配置較為繁瑣。關鍵是這個耦合UDF的更新并不是Ansys的人員執行的，而是Altair的工程師，而且更新得并不勤。

（5） 團隊協作困難：前后處理和求解器不同，你可以發現MSC和Abaqus的前身HKS公司名都是由兩三個人名的簡稱MSC（Mac Neal（下圖右方）和Robert Schwendler（下圖左方））、HKS（David Hibbitt、Bengt Karlsson、Paul Sorensen）組成，而他們公司成立一開始都是只有求解器，也就是求解器可以靠兩三個理論大牛開發，可能一個月也改不了幾行代碼

，可以使用光譜過濾對燈光效果進行后處理，得到RGB三種不同顏色下和不同RGB權重的白燈的效果。

個人沒做過實際的物理試驗，只是查了一下力矩的加載方法，一種加載方式是在梁的兩個邊緣加力，然后用力乘以離彎曲軸的距離當做力矩。但試驗室截面到底是觀測到過節點還是過形心個人理解很難說的清楚，如果問一個小朋友，看上面的彎曲圖，轉動角度θ是繞哪個軸旋轉，也許很多人就會說是底部，而不是中性軸，也就是說從外形上很難確定哪個軸旋轉，或者你可以認為是繞任意一個軸旋轉。

如果是在純彎矩情況下，由試驗得到彎曲時依然保持平面，符合實際： 如果一個橡皮條做成的懸臂梁，那么原來在鉛直線上畫的直線受力后將會變成如下圖所示的m1n1的曲線。 很明顯，上面的彎曲下原來的平面mn已經變為了曲面m1n1了。如果是這樣，不就和梁彎曲后依然為平面的假設矛盾了嗎？