ansys中實體單元的案例
ANSYS Classical 中如何獲取實體單元某截面的內力
ANSYS Classical 中如何獲取實體單元某截面的內力
相信很多童鞋在采用ANSYS進行實體單元進行分析的時候，對于如何輸出某截面的內力甚是困惑，由于實體單元的特性，ANSYS中沒有相應的集成命令來幫助我們輸出截面內力，唯一的方法只能是通過相關后處理得到我們想要的結果。
實體單元截面內力輸出，本人在這里分為兩類。
第一類：支座截面內力輸出
這種是最為簡單的內力輸出了，想要獲取支座的全部反力，我們只需輸入FSUM這個命令，即可列表顯示。如果在參數化過程中，需要提取支座反力，我們需要使用*Get命令。
例如：獲取支座X方向的反力
*get,X-force,fsum,0,item,fx
在這里我們也可以獲取一個提示，如果我們想要獲取部分支座反力，我們只需將這部分節點選取出來，然后使用上述相關命令就行了。
第二類：非支座截面的內力輸出
這類截面內力需要用到ANSYS后處理中一種比較高級的操作了，也即是面操作，核心思想在于定義結果面，將該面所包含的節點結果映射到該面上，在采用相應的積分即可得到結構內力。
下面以一個懸臂梁為例說明上述方法。
某懸臂梁，長2m，截面尺寸為300mmX500mm,混凝土等級為C30，端部固定，頂面受10KN/m的線荷載，試求端部截面和中間截面的剪力和彎矩。
展開 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的梁單元、實體單元徐變精細化分析（含各參數解釋） ￥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的實體單元扭矩施加方法總結(原創帖子，轉載請注明出處，謝謝！技術鄰ID有限元中解人生) ￥1
基于ANSYS的實體單元扭矩施加方法總結
1、 引言
在實際工程問題中，扭矩無處不在。如攻絲的絲錐、車床的光桿、攪拌軸、汽車傳動軸等等，均為受扭構件，承受扭矩作用。為了更好的分析上述構件在扭（轉）矩作用下的變形、應力、應變等物理量，現代先進設計制造分析方法引入有限元來模擬結構在外載荷作用下的響應問題。對于很多工程模型，必須考慮結構的一些幾何特征，如軸的鍵槽、絲錐的螺紋面等。因此，實體模型上扭矩的施加就成為一個非常關鍵的問題。這包括扭矩施加的形式、位置，不同方式施加的扭矩會導致整體剛度矩陣的不同，最終會導致應力奇異，影響結果的評定。ANSYS作為全球最通用的大型有限元分析軟件之一，其強大的分析功能已為國內外一致認同，現已成為許多領域結果評定的行業標準。由于ANSYS中不能直接對實體單元施加力矩，傳統方法采用若干對力偶來代替扭矩，該方法容易導致局部應力集中；改進的方法引入一些特殊單元如rbe3單元、mpc184單元、mass21單元等，通過引入這些特殊單元，能夠比較好的實現扭矩的施加，但是特殊單元的引入又改變了整體剛度矩陣。為了解決由于引入特殊單元而導致影響整體剛度矩陣的問題，有學者等提出采用接觸單元能夠很好的解決扭矩的施加問題。
本文旨在綜合關于扭矩施加的各種方法，并對這些方法進行分析比較，從而找到關于實體單元扭矩施加有效、合理的方法，為結構有限元分析提供有益的參考。
2、 ANSYS中扭矩的施加
2.1 工程實例
現以長為0.2m直徑為100mm的實心鋼管為例說明扭矩的施加。鋼管材料視為線彈性，其彈性模量及泊松比分別為：E=2e11Pa，μ=0.3。 鋼管一端固定，另一端受1000N.m扭矩作用。
展開 
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模擬該情況。
展開 ABAQUS中實體單元的應用
ABAQUS中實體單元的應用（完整可以看附件）
（此為學習時看到的一個很好的講解。 學習交流群：1063594113）
在ABAQUS的單元庫中，應用最廣泛的是應力/位移實體單元族。對三維單元，可以選擇六面體、四面體和楔形體；對二維單元則可在三角形與四邊形之間進行選擇。這些基本的單元形狀，每一種都有線性和二次的兩類選擇。對六面體和四邊形，還可選擇完全積分或減縮積分。最后，還可選用標準元或雜交元列式。另外對線性六面體或四邊形單元，還有個附加的功能，可選擇非協調模式，而對二次的三角形或四面體單元可以應用修正列式。
若列出所有種類的單元，所面臨的實體單元的總數目是相當大的，僅三維單元而言就超過20種。模擬的精度將強烈地依賴于所采用的單元類型。特別是在初次使用時，在這些單元中選擇哪一個最為合適很可能是一件令人苦惱的事情。然而，用戶會逐漸把這個工作看作是從一個20多件的工具組中，有能力選擇最恰當的工具或單元來完成的一個有價值的工作。
這一章討論了不同的單元列式和積分水平對一個特定分析的精度的影響。同時也討論了一些選擇實體單元的一般性原則。這些討論提供了獲得更多應用ABAQUS經驗和知識的基礎。在本節末的例子將允許用戶應用這些知識建立和分析一個連接柄構件的模型。
4.1 單元列式和積分
通過圖4-1所示的懸臂梁，可闡明單元階數（線性或二次），單元列式及積分水平等因素對結構模擬精度的影響。這是評估一個給定單元的性能的經典測試。因為該構件相對是細長的，我們通常用梁單元來對它建立模型。但在這里我們用這個測試來幫助評估各種實體單元的效率。
梁長150mm，寬2.5mm，高5mm；一端固定；自由端承受5N的荷載。材料的楊氏模量E為70GPa，泊松比為0.0。
展開 ABAQUS中實體單元的應用
然后，對于接觸問題，應采用細化網格的線性、減縮積分單元或者非協調單元(CAX4I，C3D8I等)。
ABAQUS中實體單元的應用.pdf
基于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巧用殼單元給實體劃分六面體網格
圖10 實體網格
到了這一步還沒有結束，最后計算時我們不需要殼單元，只需要實體單元，但是表面一層是殼單元劃分的，怎么辦呢？
先刪除表面一層的單元，是的，可以直接刪除。
圖11 網格清除
清除之后，再刪除shell181。如圖12。
圖12 刪除shell181
基于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中桿單元和殼單元的單元耦合問題
在比較復雜的結構的有限元分析中，不同的結構部件通常使用不同類型的單元來模擬。
通常情況下，不同類型的單元的各個節點的自由度數目是不同的，不同類型單元的連接節點處的自由度的耦合問題，是一個比較令人頭疼的問題。
在ANSYS中通常可以用耦合命令CP來耦合不同類型單元在連接節點處的自由度(DOF)。
也可以用CE命令來認為添加自由度之間的約束方程來達到耦合的目的。
下面是一個簡單的算例，使用了CE命令來耦合連接節點處的自由度。
模型是航天器的機翼的一個Section的某一個隔框。上下表皮是薄殼結構，用Shell63單元來模擬，在上下表皮之間有起支撐作用的桿件，用link8單元來模擬。
建模的時候，link8單元和shell63單元在連接有各自獨立的節點。即：link8單元和shell63單元的節點在連接處是重合的，但是，節點編號是各自獨立的。
link8單元在每個節點有 ux,uy,uz3個平動自由度；
shell63在每個節點有ux,uy,uz這3個平動自由度和rotx,roty,rotz這3個轉個自由，共6個自由度。
在耦合節點處，兩個耦合節點的ux,uy,uz自由度應該是相等的。
這個等式可以用CE命令來描述。
完整的命令流如下：
finish
/clear,start
/prep7
!定義第一種材料屬性；
mp,ex,1,30e6
mp,prxy,1,0.3
!定義shell63單元和實常數；
et,1,shell63
r,1,1e-3
!建立幾何模型；
rectng,31.8,33.2,0,0.3556
agen,2,1,1,1,0,0,1
a,1,4,8,5
a,6,7,3,2
KL,7,0.5, ,
KL,3,0.5, ,
在關鍵點處生成節點；
nkpt,100,4 !與編號為117的節點耦合
nkpt,101,9 !
展開 