ansys 梁單元應力的案例
基于ansys的梁單元、實體單元徐變精細化分析（含各參數解釋） ￥25
徐變是混凝土在長期恒定應力作用下產生的時變不可逆變形，其發展規律呈現前期快速增長、后期漸趨穩定的特征。主要受應力水平、材料配比、環境濕度、構件尺寸及加載齡期等因素影響。
常用方法包括有效模量法、疊加法和老化理論。國內規范（如JTG3362-2018）推薦基于線性疊加原理的徐變系數法。徐變應變可表達為：
其中， ?(t,τ)為徐變系數，需通過規范公式或實驗數據擬合確定
Ansys程序中內置金屬蠕變規律如下：
命令中詳細解釋了改公式的具體用法，以及參數意義。
二者除個別參數外形式具有異曲同工之妙，因此本案例給出用ansys精確分析混凝土徐變的方法，案例背景模擬了一個混凝土PK梁特定工況下的徐變發生過程。
案例文件中包含：
1. 00-ConcreteCreep-benchmark.mac【徐變標定文件，開箱即用，可以用來和手算對比是否正確】
2. 01-ConcreteCreep-solid.mac【分輸入模塊的參數化徐變計算文件【詳細解釋了各參數取值】。只需要改文件和計算邊界荷載即可計算實體徐變。】
3. ansa文件，用來生成網格
4. .cdb文件，網格文件
5. excel轉apdl命令流文件，用來輸入徐變系數。
進一步白話闡述一下：
1、什么是徐變？別看公式一大堆，理論一大推，簡單講就是：受力的結構，啥邊界條件、荷載不變的情況下，結構還是慢慢變形了。將這種慢慢變形的變形結果以及應力重分配準確分析出來就是徐變分析。機理一大堆，教科書上都比較詳盡，在此不做贅述，只講應用，而且是拿到案例開箱即用。
展開 WB14.0水上滑道結構應力分析（殼單元，梁單元在WB中應用）
復件 水上滑道結構應力分析報告.doc
這是我做的一個游樂設施的 仿真分析的報告，由于涉及企業的隱私，刪除了數據和名稱。
報告的關鍵點：殼單元和梁單元在workbench中的應用，殼單元和梁單元是DM建模的顯著特征，在WB中建模，分析，方便，快捷。
梁單元結構建模optistruct求解查看應力，沒有Von mises、normal stress? ￥20
本帖子是關于：整體以梁單元結構建模進行預應力模態分析，optistruct求解后查看應力結果，沒有von mises stress、normal/shear stress應力信息的原因，以及如何解決這個問題的方法。
前段時間接觸到桁架橋的結構分析，桿件橫截面主要為BOX和C型槽，C型槽的剪切中心和中性軸不重合，前處理采用梁單元cbeam建模，單元類型選擇cbar還是cbeam，可以參考：【HyperMesh寶典】之梁單元 (qq.com)。建立梁單元截面類型選擇HYPER BEAM庫下的thinwalled box和standard channel，屬性卡片選擇pbeam，求解后，hyperview查看應力結果發現只有element stress1D(s)下的CBAR/CBEAM Axial stress和long stress，沒有von mises stress、normal stress等應力。
網上搜索了一圈都沒有找到相關的問題的解決方法，也可能是我沒找全面，只能老老實實啃幫助文件，找到了關于Stress Result Written in HyperView，附上鏈接以及截圖：Stress Results Written in HyperView .h3d Format (altair.com)
展開 Ls-Dyna顯示梁單元應力
在Ls-Dyna中，如果需要在后處理中查看梁單元應力應變，不能設置為Belytschko- Schwer無應力梁（Belytschko- Schwer resultant formulation），并且需要在前處理時添加*DATABASE_EXTENT-BINARY關鍵字，將BEAMIP設置為大于0的數字，如下圖所示。
BEAMIP為梁單元積分點個數，只有在大于零的時候才才輸出梁單元應力應變。
計算完成后，可在Ls-Prepost中看到梁單元應力應變云圖。具體位置為：Fringe Component - Beam - Von Mises Stress。
展開 
關于OPTISTRUCT求解梁單元應力顯示問題
梁單元位移顯示沒有問題，應力顯示幾乎為零？請問怎樣正確顯示梁單元的應力問題，謝謝
2.jpg
1.jpg
ANSYS梁單元與實體單元的耦合與約束方程
ANSYS梁單元與實體單元的耦合與約束方程
By長安CAE
1 概述
在ANSYS計算過程中，有時候會遇到不同單元之間進行連接，由于不同的單元自由度不同，連接時通常需要通過耦合和約束方程建立節點自由度的關系，保證結果的準確性。
耦合可以理解成是將耦合的對象某個自由度作相等處理，而約束方程則不局限于相等這個關系，其可以描述具有某種關系的自由度。如圖1所示，為梁單元與平面單元的連接。如果不采用約束方程，力矩的傳遞無法完成，因為平面單元沒有轉動自由度。
圖1 梁單元與平面單元連接
為使節點2具有力矩傳遞的能力，要求1、2、3節點之間的自由度滿足以下關系：
ROTZ2 = (UY3 - UY1)/10
再通過CE命令，即可將此關系通過約束方程的形式施加給1、2、3節點。
2 命令
查看ANSYS的幫助文檔，查詢CE命令的解釋，如圖2所示。
圖2 ANSYS的CE命令解釋
CE, NEQN, CONST, NODE1, Lab1, C1, NODE2, Lab2, C2, NODE3, Lab3, C3
其中，NEQT表示常數，用于區別約束方程，一般可以用數字1、2、3表示即可，表示第幾個約束方程；
CONST表示方程的常數項，一般為0；
NODE1，表示第一個節點；
Lab1，表示自由度標簽，對于結構而言，就是三個平移和三個轉動自由度；
C1，表示該自由度的系數；
同理，后面的也一樣。
展開 ANSYS各類型單元連接專題講解(五)之3D梁單元與殼單元剛接
例如采用ANSYS模擬一個多層混凝土框架結構，一般除計算整體指標外，我們在計算具體荷載作用時（如風荷載、地震作用、恒載、活載等），樓板一般采用彈性版，此時可用殼單元模擬，主梁、次梁采用梁單元模擬，此時變為梁單元包含在殼面內的情況，當然此類情況是否需要考慮截面偏置，可根據具體工程而定。
對這中梁單元包含在殼單元面內的情況，只需要將梁單元與殼單元共用節點即可，而無須格外建立約束方程。
三、梁單元在殼單元內但不包含
此種情況為梁與殼位于同一面內，但其中面不包含梁線，適用于多尺度建模分析（如下圖）。梁單元與殼單元的連接在端部可以通過剛性梁和剛性區域兩種方式連接。剛性梁采用MPC184單元，剛性區域采用Cerig命令，具體使用方法下期文章討論。
展開 ANSYS各類型單元連接專題講解（四）之2D梁與2D實體單元剛接
前面文章主要講解了梁單元與其他類型單元鉸接的情況，從本篇文章開始，將主要講解梁單元與各類單元剛接的情況，而這也是我們日常工程中比較常見的一種連接方式。
首先從2D平面單元單元開始說起。
盡管現在的ANSYS版本已經摒棄了很古老的2D梁單元，改用Beam18x系列單元代替，但為究其連接方法，這類方面仍具有一定的講解價值，例如我們計算一榀框架的時候多數時候是采用2D平面單元的。
2D梁單元包括：beam3、beam23、beam54
2D實體單元：plane單元
一般來講，2D梁單元與2D實體單元剛接一般分為三種方法：
1）約束方程法；2）偽梁法；3）MPC法。
三種方法的連接原理無非是建立自由度之間的關系方程，但值得注意的是由于采用了局部區域的節點，因而在建立關系的局部區域內可能會有應力集中的情況，后處理當中應格外注意。
約束方程法后續講解3D梁單元連接時會詳細說明，此處簡單說下偽梁法與MPC法。
其實偽梁法與MPC法原理基本一樣，構造一個虛擬梁單元，使虛擬梁單元與外部梁單元剛接，虛擬梁單元與內部實體單元強制剛接，從而間接實現外部梁單元與實體單元的剛接效果。
使用偽梁法需注意的是，在建立虛擬梁單元時，虛擬梁單元應至少與實體單元的兩個節點相連，剛度可取為無窮大或者實際梁單元的10^5倍。
下面以一個小案例來演示。
如上圖所示，兩塊小鋼板中間靠一小鋼梁連接，小鋼梁上有均布荷載，尺寸如上所示，均以mm計，中間鋼梁所受均布荷載為10KN/m,采用ANSYS模擬該情況。
展開 Hypermesh為ansys創建梁單元（一） ￥1
Hypermesh與ansys聯合仿真系列之Hypermesh為ansys創建梁單元（一）。
本文介紹ansys梁單元中的beam188和beam189及它們之間的本質區別，以及仿真時對兩種梁單元的選擇建議。簡介梁單元的關鍵字設置及截面設置，主要介紹通過Hypermesh在ansys求解器下兩種創建梁單元的詳細步驟及效果對比。
ANSYS各類型單元連接專題講解（二）之桿與梁殼體單元的連接
前一篇文章主要介紹了單元之間連接的主要原則，今天開始主要從具體方面講解連接方法。
按照桿、梁、殼、實體的順序，先說說桿單元與各單元的連接方法。
那么什么時候需要用到桿單元與各種單元的連接呢？水哥稍微列舉下實際工程中需要考慮此類連接的例子。
案例一：工業廠房
此類結構一般橫向跨度較大，屋頂采用鋼結構形式，在具體模擬屋架時，此時各個桿件可看成鉸接，采用桿單元模擬。而下方框架柱則采用梁單元進行模擬，在相交部位則需要用到桿單元與梁單元的連接。
案例二：門廳鋼結構雨棚
在具體模擬該結構時，雨棚上方拉桿采用桿單元模擬，而下方的鋼梁采用梁單元模擬，混凝土框架柱可采用實體單元模擬。
一直以來，桿單元一般用于模擬桁架結構的時候比較多，其特點是桿件兩端不考慮承受彎矩作用，節點只有平動自由度，是所有單元中最為簡單的一種。
桿單元分為2D桿單元和3D桿單元，2D桿單元節點只有Ux和Uy兩個平動自由度，而3D桿單元除了這兩個，還有Uz。其他單元，梁單元、殼單元、體單元都包含了這三個自由度，且具有相同的物理意義，按照前面一篇文章所介紹的連接總則，桿單元與其他單元連接時只需要共用節點即可，無需建立約束方程。
下面是一個簡單的類似雨棚案例，注意本案例各構件尺寸僅為演示操作需要所擬，未經仔細推敲，各工程大佬可忽略。
某屋外雨棚平面簡化模型如上，長度為4m，折算荷載為10 KN/m，雨棚梁采用工字型鋼I40，系桿截面面積為238.64mm^2，材料均為Q235，采用ANSYS模擬該結構。
下面為建模過程
!
展開 天方地圓結構-梁殼單元建模實例！再次驗證應力奇異的可怕性！
應力集中和應力奇異的驗證
▲▲▲
因此模型存在矩形截面，即有90°夾角部分，從理論上來說，此部分位置必然會產生應力集中，筆者通過不斷細化網格求解驗證發現，網格越細化，應力值則越大，如果僅僅出現應力集中的話，即使網格不斷細化，那么最終的應力值也會趨于一個穩定值，但此模型則是隨著網格的不斷細化，應力值似乎有無限變大的趨勢，所以此時，應力集中的部位已然產生了應力奇異，但此應力奇異是實際模型本身的幾何缺陷造成的，因而是無可避免的。
因而，對于有幾何結構突變部位并非網格越細化越好的，網格細化卻帶來了無法判斷結果準確性的一個可怕后果，建議：
1. 在建模的時候如果存在幾何突變部位，務必要采用倒圓角進行圓滑過渡，以避免應力集中變成了應力奇異（即使真實模型中并未有圓角過渡）。
2. 梁殼單元在突變部位隨網格細化產生應力奇異是無可避免的，所以盡量避免在有幾何突變的部位采用梁殼單元，最好采用實體單元。
3. 實體建模時，四面體網格也會隨網格細化產生無可避免的應力奇異，所以最好采用六面體網格，這一點壓力容器行業基本上都采用六面體網格。
應力分析僅僅有漂亮的計算云圖和計算結果是遠遠不夠的，分析設計人必須具有判斷計算結果準確性的能力，試想如果以并非準確的應力結果作為評定依據，那將是多么可怕的一件事情！
有從事壓力容器分析設計行業的朋友如感興趣可掃描上方二維碼關注一下公眾號！
展開 
ANSYS梁單元自定義截面
梁單元作為一種簡單且高效的計算單元，在結構分析尤其是建筑結構中得到廣泛的應用。使用梁單元可以避免將結構中梁柱全部轉換為實體單元，從而降低了計算量，且梁單元結構形式簡單，求解精度也相對較高。在ANSYS中，梁單元基本上可以分為線性單元和二次單元，二者之間計算理論不同，經典的二次單元即BEAM189單元的積分點如下圖所示：
在ANSYS中可以為BEAM單元定義截面，其中大部分經典的截面形式都包含在ANSYS的截面庫中，但是經典的梁單元計算時截面方向分為四個單元，這對于一般計算來說是足夠的，但如果需要仔細分析截面方向的內力，可能就略顯的粗糙了。除此之外，鋼管混凝土、組合梁之類也都是異形梁截面，此時標準截面庫中的數據也沒什么用。針對這個問題存在兩種解決方式，一種是使用ASEC自定義截面參數，這個命令不管截面如何，只需要給出截面相關的信息即可，截面的信息輸入如下圖所示：
至于這些截面的參數可以使用簡單的截面計算工具得到，如果是鋼筋混凝土梁這種比較復雜的復合梁，那么需要使用Xtract之類的截面有限元軟件進行計算。將截面信息填入。采用ASEC的截面輸入方式計算效率高，截面信息準確的話，精度也不差，但缺點是不能輸出截面積分點和柵點的數據。
另一種方式就是自定義截面，其基本思路如下：
1.設定MESH200單元，建立截面幾何形狀；
2.用MESH200單元劃分截面，并保存截面數據；
3.建立計算幾何模型，讀取截面數據；
4.賦予模型截面，施加邊界條件計算；
5.后處理。
展開 ANSYS梁單元自定義截面
ANSYS梁單元自定義截面
梁單元作為一種簡單且高效的計算單元，在結構分析尤其是建筑結構中得到廣泛的應用。使用梁單元可以避免將結構中梁柱全部轉換為實體單元，從而降低了計算量，且梁單元結構形式簡單，求解精度也相對較高。在ANSYS中，梁單元基本上可以分為線性單元和二次單元，二者之間計算理論不同，經典的二次單元即BEAM189單元的積分點如下圖所示：
在ANSYS中可以為BEAM單元定義截面，其中大部分經典的截面形式都包含在ANSYS的截面庫中，但是經典的梁單元計算時截面方向分為四個單元，這對于一般計算來說是足夠的，但如果需要仔細分析截面方向的內力，可能就略顯的粗糙了。除此之外，鋼管混凝土、組合梁之類也都是異形梁截面，此時標準截面庫中的數據也沒什么用。針對這個問題存在兩種解決方式，一種是使用ASEC自定義截面參數，這個命令不管截面如何，只需要給出截面相關的信息即可，截面的信息輸入如下圖所示：
至于這些截面的參數可以使用簡單的截面計算工具得到，如果是鋼筋混凝土梁這種比較復雜的復合梁，那么需要使用Xtract之類的截面有限元軟件進行計算。將截面信息填入。采用ASEC的截面輸入方式計算效率高，截面信息準確的話，精度也不差，但缺點是不能輸出截面積分點和柵點的數據。
另一種方式就是自定義截面，其基本思路如下：
1.設定MESH200單元，建立截面幾何形狀；
2.用MESH200單元劃分截面，并保存截面數據；
3.建立計算幾何模型，讀取截面數據；
4.賦予模型截面，施加邊界條件計算；
5.后處理。
展開 hypermesh-ansys聯合仿真-《梁單元4》 ￥1
在《hypermesh-ansys聯合仿真-梁單元3》中對比了梁單元和實體單元的結果，表明梁單元計算結果更容易接近理論計算值，且付出的計算資源是很小的。但并非所有情況都是這樣，下面介紹一種情況實例來說明問題。
如圖兩端固支的C型薄壁梁，在梁中心位置作用一個F=100N的集中力，具體作用點是C型截面的上邊沿（上右圖），下面分別采用梁單元和殼單元分別計算該結構工況下梁的變形梁，讀者可以自行計算嘗試并分析哪種結算結果更可靠？造成這個結果的原因是什么？我們如何在梁單元與殼單元之間做選擇
梁截面尺寸
展開 hypermesh-ansys聯合仿真-《梁單元1》
前文已經通過《hypermesh-ansys聯合仿真-基本步驟1》系列詳細說明了hypermesh與ansys聯合仿真時的基本過程，下面通過一列文章按照單元類型分別介紹不同單元類型的建模方法以及使用這些單元時需要注意的問題，當忽略這些問題時往往會造成較大的誤差甚至錯誤。
梁單元簡介
當結構長度方向尺寸明顯大于截面尺寸時（常常設定為10:1），我們可以將結構簡化為一維的梁單元，相比于三維的實體單元可以在保證求解精度的情況下大大降低計算量。ANSYS中3D分析中的梁單元主要采用Beam188和Beam189，均被默認為鐵摩辛柯（Timoshenko）梁單元，該假設認為變形后橫截面保持平面且不發生扭曲，可以計算彎曲、拉壓、扭轉并且考慮剪切變形的影響，兩種梁單元可以滿足大多數工程問題。
Beam188 單元具有I和J兩個節點，每個節點有6-7個自由度。在高版本中可以通過KEYOPT(3)選擇插值函數，默認KEYOPT(3)=0，單元在梁長度方向只有一個積分點，因此在采用SMISC獲取節點I和J的結果時，以重心的結果表示兩個節點的結果，導致彎矩圖呈現鋸齒狀；如果KEYOPT(3)=2，ANSYS采用增加一個內部節點（用戶無法訪問）沿著梁長設置兩個積分點，這樣彎矩圖沿著梁線性變化不再呈鋸齒狀。當KEYOPT(3)=2時除了初始幾何為直線（不管是否采用二次插值）和不能訪問增加的節點外其他和Beam189相同。目前2019版本的Hypermesh可以為ANSYS求解器設置線性插值、二次插值和三次插值。Beam188單元支持線性分析，大轉動、大應變等非線性分析，支持彈性材料定義， 塑性、蠕變以及其余非線性材料定義。
展開 