ansys實體單元含義的案例
基于ansys的梁單元、實體單元徐變精細化分析（含各參數解釋） ￥25
徐變應變可表達為：
其中， ?(t,τ)為徐變系數，需通過規范公式或實驗數據擬合確定
Ansys程序中內置金屬蠕變規律如下：
命令中詳細解釋了改公式的具體用法，以及參數意義。
二者除個別參數外形式具有異曲同工之妙，因此本案例給出用ansys精確分析混凝土徐變的方法，案例背景模擬了一個混凝土PK梁特定工況下的徐變發生過程。
案例文件中包含：
1. 00-ConcreteCreep-benchmark.mac【徐變標定文件，開箱即用，可以用來和手算對比是否正確】
2. 01-ConcreteCreep-solid.mac【分輸入模塊的參數化徐變計算文件【詳細解釋了各參數取值】。只需要改文件和計算邊界荷載即可計算實體徐變。】
3. ansa文件，用來生成網格
4. .cdb文件，網格文件
5. excel轉apdl命令流文件，用來輸入徐變系數。
進一步白話闡述一下：
1、什么是徐變？別看公式一大堆，理論一大推，簡單講就是：受力的結構，啥邊界條件、荷載不變的情況下，結構還是慢慢變形了。將這種慢慢變形的變形結果以及應力重分配準確分析出來就是徐變分析。機理一大堆，教科書上都比較詳盡，在此不做贅述，只講應用，而且是拿到案例開箱即用。
白話闡述要點：
1、案例是ansys apdl（命令流）分析的，給出了全套參數化命令流，材料模型定義、材料參數定義、求解，拿過來可以直接運行。
2、機理是用了ansys中關于金屬蠕變的材料模型。（細想蠕變和徐變的現象，表征都是一樣的。至于機理，各有各的理論，但不影響材料模型使用。）
具體使用：
1、，先跑一遍，看看到底徐變是怎么個事兒。
展開 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 APDL實體單元和殼單元（不共節點）之間的連接 ￥100
實體單元和殼單元之間的連接是ANSYS中常見的問題。即使兩種單元之間共節點，但單元之間不連續（實體單元每個節點有3個平動自由度，而殼單元每個節點有3個平動自由度和3個轉動自由度），對于兩種單元之間面面接觸，可直接定義剛域，本文主要采用MPC法對實體-殼單元的連接方法進行說明。
1 單元類型
算例模型中，實體單元采用SOLID45，殼單元采用SHELL63，接觸位置不共節點。對于兩種單元之間的連接，通過目標單元TARGE170和接觸單元CONTA175實現，定義約束為實體-殼約束，接觸單元為MPC算法，接觸類型為綁定接觸。
2 有限元模型和綁定接觸
圖1 底部固定約束，殼單元施加均布荷載
圖2 目標單元和接觸單元
3 計算結果
圖3 von Mises stress
圖4 X-Component of displacement
付費內容為相關命令流。
展開 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模擬該情況。
展開 
基于ANSYS經典界面的實體-板單元連接建模
1.8 創建第2個接觸
ASEL,S,,,47
NSLA,S,1
TYPE,3
REAL,3
ESURF
ALLS
LSEL,S,,,6
NSLL,S,1
TYPE,4
REAL,3
ESURF
ALLS
上述命令在此線和上述面之間建立接觸對
創建的接觸對如下
1.9 創建第3個接觸
ASEL,S,,,47
NSLA,S,1
TYPE,3
REAL,4
ESURF
ALLS
LSEL,S,,,7
NSLL,S,1
TYPE,4
REAL,4
ESURF
ALLS
上述命令在此線和上述面之間建立接觸對
創建的接觸對如下
1.10 創建第4個接觸
ASEL,S,,,47
NSLA,S,1
TYPE,3
REAL,5
ESURF
ALLS
LSEL,S,,,8
NSLL,S,1
TYPE,4
REAL,5
ESURF
ALLS
上述命令在此線和上述面之間建立接觸對
創建的接觸對如下
1.11 生成附加的虛擬殼單元
SHSD,2,CREATE
該命令是關鍵命令，它用于創建殼-實體裝配中的殼-實體接口。
其含義是對于2號實常數所定義的接觸對，創建其殼-實體接口。
2號實常數所定義的接觸對，就是下圖所示的接觸對
運行上述命令后，ANSYS在接觸線處創建了虛擬的殼單元
并在虛擬單元的每個節點處（如果此處沒有沒有接觸單元的話）創建了新的接觸單元。
通過這種方式，在實體和殼單元之間建立了關系。
SHSD,3,CREATE
SHSD,4,CREATE
SHSD,5,CREATE
FINISH
上述命令依次在其它三處創建這種連接關系，
然后退出前處理。
展開 ANSYS巧用殼單元給實體劃分六面體網格
圖10 實體網格
到了這一步還沒有結束，最后計算時我們不需要殼單元，只需要實體單元，但是表面一層是殼單元劃分的，怎么辦呢？
先刪除表面一層的單元，是的，可以直接刪除。
圖11 網格清除
清除之后，再刪除shell181。如圖12。
圖12 刪除shell181
ANSYS Classical 中如何獲取實體單元某截面的內力
ANSYS Classical 中如何獲取實體單元某截面的內力
相信很多童鞋在采用ANSYS進行實體單元進行分析的時候，對于如何輸出某截面的內力甚是困惑，由于實體單元的特性，ANSYS中沒有相應的集成命令來幫助我們輸出截面內力，唯一的方法只能是通過相關后處理得到我們想要的結果。
實體單元截面內力輸出，本人在這里分為兩類。
第一類：支座截面內力輸出
這種是最為簡單的內力輸出了，想要獲取支座的全部反力，我們只需輸入FSUM這個命令，即可列表顯示。如果在參數化過程中，需要提取支座反力，我們需要使用*Get命令。
例如：獲取支座X方向的反力
*get,X-force,fsum,0,item,fx
在這里我們也可以獲取一個提示，如果我們想要獲取部分支座反力，我們只需將這部分節點選取出來，然后使用上述相關命令就行了。
第二類：非支座截面的內力輸出
這類截面內力需要用到ANSYS后處理中一種比較高級的操作了，也即是面操作，核心思想在于定義結果面，將該面所包含的節點結果映射到該面上，在采用相應的積分即可得到結構內力。
下面以一個懸臂梁為例說明上述方法。
某懸臂梁，長2m，截面尺寸為300mmX500mm,混凝土等級為C30，端部固定，頂面受10KN/m的線荷載，試求端部截面和中間截面的剪力和彎矩。
展開 ANSYS輸出實體模型表面的節點信息 和單元拓撲關系
ANSYS輸出實體模型表面的節點信息
和單元拓撲關系
遇到一個問題，一個給定的實體模型，劃分了solid185的單元，假如實體模型單元劃分如下。需要提取實體模型外表面節點位置信息和單元拓撲關系（也就是每一個單元是由哪幾個節點組成的），目的是方便做其他分析，比如流體分析，提取外表面的節點可以施加溫度載荷。
圖1
對于此問題，在ansys里面很難直接提取所有外表面的節點和單元信息，因為外表面也是實體單元的一個單元面，不可能剝離出來。
因此，想要提取外表面的單元和節點，最好是需要外表面存在平面單元。
對于此，可以采用ansys里面的特殊單元mesh200，這個單元用于面網格的劃分，而且劃分后的單元不參與實際計算。
于是：
et,2,200 ！定義mesh200單元類型
asel,s,ext ！選擇所有的外表面
aatt,,,2 ! 設置劃分單元為mesh200
KEYOPT, 2, 1, 6 ! 4節點的四邊形單元
amesh,all ! 劃分所以的外表面
此時劃分的面網格和原來的實體網格的節點是一一對應的，這就保證了最后輸出的節點的坐標與原來實體模型的對應節點是一一對應的。
此時可以選擇刪除實體模型和實體單元。
展開 基于ANSYS經典界面的實體-板單元連接建模
1.8 創建第2個接觸
ASEL,S,,,47
NSLA,S,1
TYPE,3
REAL,3
ESURF
ALLS
LSEL,S,,,6
NSLL,S,1
TYPE,4
REAL,3
ESURF
ALLS
上述命令在此線和上述面之間建立接觸對
創建的接觸對如下
1.9 創建第3個接觸
ASEL,S,,,47
NSLA,S,1
TYPE,3
REAL,4
ESURF
ALLS
LSEL,S,,,7
NSLL,S,1
TYPE,4
REAL,4
ESURF
ALLS
上述命令在此線和上述面之間建立接觸對
創建的接觸對如下
1.10 創建第4個接觸
ASEL,S,,,47
NSLA,S,1
TYPE,3
REAL,5
ESURF
ALLS
LSEL,S,,,8
NSLL,S,1
TYPE,4
REAL,5
ESURF
ALLS
上述命令在此線和上述面之間建立接觸對
創建的接觸對如下
1.11 生成附加的虛擬殼單元
SHSD,2,CREATE
該命令是關鍵命令，它用于創建殼-實體裝配中的殼-實體接口。
其含義是對于2號實常數所定義的接觸對，創建其殼-實體接口。
2號實常數所定義的接觸對，就是下圖所示的接觸對
運行上述命令后，ANSYS在接觸線處創建了虛擬的殼單元
并在虛擬單元的每個節點處（如果此處沒有沒有接觸單元的話）創建了新的接觸單元。
通過這種方式，在實體和殼單元之間建立了關系。
SHSD,3,CREATE
SHSD,4,CREATE
SHSD,5,CREATE
FINISH
上述命令依次在其它三處創建這種連接關系，
然后退出前處理。
展開 ANSYS巧用殼單元給實體劃分六面體網格
圖10 實體網格
到了這一步還沒有結束，最后計算時我們不需要殼單元，只需要實體單元，但是表面一層是殼單元劃分的，怎么辦呢？
先刪除表面一層的單元，是的，可以直接刪除。
圖11 網格清除
清除之后，再刪除shell181。如圖12。
圖12 刪除shell181
ANSYS巧用殼單元給實體劃分六面體網格
圖10 實體網格
到了這一步還沒有結束，最后計算時我們不需要殼單元，只需要實體單元，但是表面一層是殼單元劃分的，怎么辦呢？
先刪除表面一層的單元，是的，可以直接刪除。
圖11 網格清除
清除之后，再刪除shell181。如圖12。
圖12 刪除shell181
Ansys中級認證窗口課程：LS-DYNA中殼體與實體單元連接技術應用
摘要：在LS-DYNA分析中經常會使用實體單元與殼體單元以滿足不同部位的分析要求，這就存在殼與實體單元連接時自由度不匹配的問題。本文詳述三種不同的連接方法案例。如果不需要傳遞轉動可以使用合并節點法和約束法，合并節點法要求節點重合，計算效率最高，約束法不要求節點重合。接觸法可以傳遞轉動，接觸法使用最為靈活，消耗的計算資源較多。
殼體單元的每個節點只有３個沿著x、y和z方向的平動自由度UX、UY、UZ；在實體單元中，每個節點具有六個自由度：沿x、y 和z方向的平動自由度UZ、UY、UZ以及繞X、Y和Z軸的轉動自由度TOTX、TOTY、ROTZ。當實體單元和殼單元連接在一起共同工作時，即存在自由度不協調問題。
案例部分分為四步，第一步建立沒有連接的模型，后三步都是在第一步模型的基礎上進行連接。具體操作視頻請在技術鄰搜索“李安民”，關注我，收看視頻。
1.1 模型建立
1.1.1 幾何模型
Geometry->Solid->Box，在Creat Box對話框或者圖形視口（Graphics Viewport）輸入實體單元尺寸，如果所示，點擊Apply關閉完成長方體。
Geometry->Surface->Plane，在Create Plane輸入如下圖所示的參數，點擊Appley生成平面。
1.1.2 網格劃分
FEM->Element and Mesh->Solid Mesher對實體網格劃分，填入Elem Size為0.5，點擊Try Meshing Automatically，若不滿意可以點擊Reject拒絕，再從新調整尺寸等參數，確認無誤，點擊Accept。
展開 ANSYS Workbench如何獲取實體單元某截面的剪力和軸力
其實在workbench中原理和classi一樣，也是通過建立結構面，投影結果并積分求解而來，依然以前文的懸臂梁為例來簡單說明。
第一步：建立該懸臂梁，并求解，如下。
第二步，在我們需要查看內力截面位置處建立局部坐標系，這里我們建立中間位置處截面，如下所示：
第三步，建立求解面。選擇Model,可以在工具欄中選中construction geometry，插入surface，坐標系選擇我們剛才建立的坐標系。
第四步，提取各個應力，也即是投影節點應力到我們的面上。選擇我們需要投影的節點力，點擊worksheet,然后在表格中右鍵 create user defined results.這兒我們提取SXZ和SZ，來獲取我們Y方向剪力和Z方向軸力。
第五步，觀察應力，并計算內力。
注意提取的時候要注意選擇SURFACE。
SXZ應力分布：
SZ應力分布：
我們這個截面的最終內力也即是 該截面的平均應力乘以我們的面積。
比如：
剪力
FY=66667*0.3*0.5N=10KN
這是與理論結果較為符合的。
細心的小伙伴可能會發現，為什么這里只說了WORKBENCH獲取軸力和剪力的方法，彎矩怎么獲取呢？
因為水哥也還不知道~~~場面一度十分尷尬。有興趣的歡迎可以一起研究討論哦~~~
展開 基于ANSYS的實體單元扭矩施加方法總結(原創帖子，轉載請注明出處，謝謝！技術鄰ID有限元中解人生) ￥1
基于ANSYS的實體單元扭矩施加方法總結
1、 引言
在實際工程問題中，扭矩無處不在。如攻絲的絲錐、車床的光桿、攪拌軸、汽車傳動軸等等，均為受扭構件，承受扭矩作用。為了更好的分析上述構件在扭（轉）矩作用下的變形、應力、應變等物理量，現代先進設計制造分析方法引入有限元來模擬結構在外載荷作用下的響應問題。對于很多工程模型，必須考慮結構的一些幾何特征，如軸的鍵槽、絲錐的螺紋面等。因此，實體模型上扭矩的施加就成為一個非常關鍵的問題。這包括扭矩施加的形式、位置，不同方式施加的扭矩會導致整體剛度矩陣的不同，最終會導致應力奇異，影響結果的評定。ANSYS作為全球最通用的大型有限元分析軟件之一，其強大的分析功能已為國內外一致認同，現已成為許多領域結果評定的行業標準。由于ANSYS中不能直接對實體單元施加力矩，傳統方法采用若干對力偶來代替扭矩，該方法容易導致局部應力集中；改進的方法引入一些特殊單元如rbe3單元、mpc184單元、mass21單元等，通過引入這些特殊單元，能夠比較好的實現扭矩的施加，但是特殊單元的引入又改變了整體剛度矩陣。為了解決由于引入特殊單元而導致影響整體剛度矩陣的問題，有學者等提出采用接觸單元能夠很好的解決扭矩的施加問題。
本文旨在綜合關于扭矩施加的各種方法，并對這些方法進行分析比較，從而找到關于實體單元扭矩施加有效、合理的方法，為結構有限元分析提供有益的參考。
2、 ANSYS中扭矩的施加
2.1 工程實例
現以長為0.2m直徑為100mm的實心鋼管為例說明扭矩的施加。鋼管材料視為線彈性，其彈性模量及泊松比分別為：E=2e11Pa，μ=0.3。 鋼管一端固定，另一端受1000N.m扭矩作用。
展開 