ansys鋼板用實體單元的案例
ANSYS巧用殼單元給實體劃分六面體網格
ANSYS巧用殼單元給實體劃分六面體網格
1 概述
眾所周知，ANSYS經典劃分網格的功能比較弱，映射劃分（Map）和掃掠劃分(Sweep)對幾何形狀的要求都十分高。而四面體網格一方面導致單元數目多余六面體，一方面給計算后處理帶來一定的不便。
有些情況下，幾何模型的結構導致即使再怎么對模型進行切分都不可能掃掠出六面體網格，這種情況下，可以巧妙地利用殼單元。
ANSYS經典里對于一個平面，劃分網格非常簡單，而且幾何形狀約束很少，即使是自由劃分的網格，一般情況下都比較規整。利用這個特點，用殼單元對面進行網格劃分，然后再對整個實體模型進行網格劃分。本次以一個例子示意此過程。
2 過程
首先在ANSYS經典界面定義兩個單元類型，分別是shell181和solid185。如圖1所示。
圖1 單元類型
建立幾何模型，采用block命令，建立100x40x10的長方體：
block,-50,50,0,10,-20,20
如圖2所示。
圖2 幾何模型
之后為了演示網格劃分，將模型切分成幾塊，如圖3所示。
圖3 切分模型
再然后選擇殼單元shell181，如圖4所示：
圖4 選擇shell181單元
然后設置模型最上層的面各個線條的分數：
圖5 操作
份數分別如圖6.
圖6 線條份數
之后點擊MeshTool，如圖7所示。
圖7 劃分面網格設置
如圖7設置，點擊Mesh，選中模型的最上一層表面劃分，得到圖8的結果。
圖8 面網格劃分
再重復前面的選擇單元的操作，選擇單元類型為solid185，并且在MeshTool里選擇Volumes 的掃掠（sweep）劃分，如圖9所示。
圖9 操作
點擊Sweep選中所有的體，即得到如圖10所示的網格。
展開 ANSYS巧用殼單元給實體劃分六面體網格
ANSYS巧用殼單元給實體劃分六面體網格
1 概述
眾所周知，ANSYS經典劃分網格的功能比較弱，映射劃分（Map）和掃掠劃分(Sweep)對幾何形狀的要求都十分高。而四面體網格一方面導致單元數目多余六面體，一方面給計算后處理帶來一定的不便。
有些情況下，幾何模型的結構導致即使再怎么對模型進行切分都不可能掃掠出六面體網格，這種情況下，可以巧妙地利用殼單元。
ANSYS經典里對于一個平面，劃分網格非常簡單，而且幾何形狀約束很少，即使是自由劃分的網格，一般情況下都比較規整。利用這個特點，用殼單元對面進行網格劃分，然后再對整個實體模型進行網格劃分。本次以一個例子示意此過程。
2 過程
首先在ANSYS經典界面定義兩個單元類型，分別是shell181和solid185。如圖1所示。
圖1 單元類型
建立幾何模型，采用block命令，建立100x40x10的長方體：
block,-50,50,0,10,-20,20
如圖2所示。
圖2 幾何模型
之后為了演示網格劃分，將模型切分成幾塊，如圖3所示。
圖3 切分模型
再然后選擇殼單元shell181，如圖4所示：
圖4 選擇shell181單元
然后設置模型最上層的面各個線條的分數：
圖5 操作
份數分別如圖6.
圖6 線條份數
之后點擊MeshTool，如圖7所示。
圖7 劃分面網格設置
如圖7設置，點擊Mesh，選中模型的最上一層表面劃分，得到圖8的結果。
圖8 面網格劃分
再重復前面的選擇單元的操作，選擇單元類型為solid185，并且在MeshTool里選擇Volumes 的掃掠（sweep）劃分，如圖9所示。
圖9 操作
點擊Sweep選中所有的體，即得到如圖10所示的網格。
展開 ANSYS巧用殼單元給實體劃分六面體網格
ANSYS巧用殼單元給實體劃分六面體網格
1 概述
眾所周知，ANSYS經典劃分網格的功能比較弱，映射劃分（Map）和掃掠劃分(Sweep)對幾何形狀的要求都十分高。而四面體網格一方面導致單元數目多余六面體，一方面給計算后處理帶來一定的不便。
有些情況下，幾何模型的結構導致即使再怎么對模型進行切分都不可能掃掠出六面體網格，這種情況下，可以巧妙地利用殼單元。
ANSYS經典里對于一個平面，劃分網格非常簡單，而且幾何形狀約束很少，即使是自由劃分的網格，一般情況下都比較規整。利用這個特點，用殼單元對面進行網格劃分，然后再對整個實體模型進行網格劃分。本次以一個例子示意此過程。
2 過程
首先在ANSYS經典界面定義兩個單元類型，分別是shell181和solid185。如圖1所示。
建立幾何模型，采用block命令，建立100x40x10的長方體：
block,-50,50,0,10,-20,20
如圖2所示。
圖2 幾何模型
之后為了演示網格劃分，將模型切分成幾塊，如圖3所示。
圖3 切分模型
再然后選擇殼單元shell181，如圖4所示：
圖4 選擇shell181單元
然后設置模型最上層的面各個線條的分數：
圖5 操作
份數分別如圖6.
圖6 線條份數
之后點擊MeshTool，如圖7所示。
圖7 劃分面網格設置
如圖7設置，點擊Mesh，選中模型的最上一層表面劃分，得到圖8的結果。
圖8 面網格劃分
再重復前面的選擇單元的操作，選擇單元類型為solid185，并且在MeshTool里選擇Volumes 的掃掠（sweep）劃分，如圖9所示。
圖9 操作
點擊Sweep選中所有的體，即得到如圖10所示的網格。
展開 基于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經典界面的實體-板單元連接建模
（2）中間空心部分使用殼單元，邊上實心部分使用實體單元。
（3）上述兩種單元需要建立連接關系。實心單元每個節點有3個自由度，而殼單元每個節點有6個自由度，如何建立連接關系呢？ANSYS提供了SHSD命令來建立這種連接。要使用該命令，首先需要創建接觸對，并且要對目標-接觸單元的關鍵字進行設置。下面的絕大多數操作都是圍繞該命令進行的。
【求解步驟】
1.前處理
1.1 創建單元
/PREP7
ET,1,SOLID187
ET,2,SHELL181
ET,3,TARGE170
KEYOPT,3,5,1
ET,4,CONTA175
KEYOPT,4,2,2
KEYOPT,4,12,5
上述命令分別定義了4種單元。
第1種是實體單元，第2種是殼單元，他們分別用于建模上述梁的實體部分和空心部分。
第3-4種則是用于模擬接觸部分，就是實體與空心的接觸部分。
這里對于這兩種單元均設置了關鍵字，這些關鍵字的設置是使用后面的命令“SHSD”所必須的。
1.2 創建實常數
R,1,0.02
R,2
R,3
R,4
R,5
這里創建了5個實常數。
第1個實常數用于定義空心梁的厚度
第2-5個實常數分別用于定義4個接觸對。
1.2 創建材料類型
MP,EX,1,2e11
MP,PRXY,1,0.3
上述命令定義了材料的彈性模量和泊松比。
1.3 創建中間的空心梁
/VIEW,1,1,1
BLOCK,-0.14,0.14,-0.14,0.14,0,0.98
VDELE,1,,,0
ADELE,1,2,1,1
上述命令首先創建了一個長方體，
然后刪除了體本身，留下構成長方體的面，線和關鍵點。
最后又刪除了兩端的面。
結果如下圖。
展開 ANSYS Classical 中如何獲取實體單元某截面的內力
ANSYS Classical 中如何獲取實體單元某截面的內力
相信很多童鞋在采用ANSYS進行實體單元進行分析的時候，對于如何輸出某截面的內力甚是困惑，由于實體單元的特性，ANSYS中沒有相應的集成命令來幫助我們輸出截面內力，唯一的方法只能是通過相關后處理得到我們想要的結果。
實體單元截面內力輸出，本人在這里分為兩類。
第一類：支座截面內力輸出
這種是最為簡單的內力輸出了，想要獲取支座的全部反力，我們只需輸入FSUM這個命令，即可列表顯示。如果在參數化過程中，需要提取支座反力，我們需要使用*Get命令。
例如：獲取支座X方向的反力
*get,X-force,fsum,0,item,fx
在這里我們也可以獲取一個提示，如果我們想要獲取部分支座反力，我們只需將這部分節點選取出來，然后使用上述相關命令就行了。
第二類：非支座截面的內力輸出
這類截面內力需要用到ANSYS后處理中一種比較高級的操作了，也即是面操作，核心思想在于定義結果面，將該面所包含的節點結果映射到該面上，在采用相應的積分即可得到結構內力。
下面以一個懸臂梁為例說明上述方法。
某懸臂梁，長2m，截面尺寸為300mmX500mm,混凝土等級為C30，端部固定，頂面受10KN/m的線荷載，試求端部截面和中間截面的剪力和彎矩。
展開 ANSYS輸出實體模型表面的節點信息 和單元拓撲關系
刪除后，顯示的就只有外表面的單元和節點，通過菜單欄List下面的Nodes，可選擇列表輸出節點坐標，如下所示：
圖2
同樣，可通過List輸出單元拓撲信息，如下：
先用命令：
numcmp,all
壓縮所有的單元和節點的編號，也就是重新編號，否則可能不是從1開始。
最后顯示如下：
圖3
這兩個文件都可以保存，之后需要可以再寫程序提取。
這里有一個問題，現在得到的表面的節點號和原來實體模型對于位置的節點號不是對應的，處理這個問題需要重新寫程序，用什么語言都可以，Python，C++等等，目的是讀取ansys輸出的節點信息文件，讀出固定坐標處對應的節點號，通過對比可以找到所以的節點對應關系。
展開 
基于ANSYS經典界面的實體-板單元連接建模
（2）中間空心部分使用殼單元，邊上實心部分使用實體單元。
（3）上述兩種單元需要建立連接關系。實心單元每個節點有3個自由度，而殼單元每個節點有6個自由度，如何建立連接關系呢？ANSYS提供了SHSD命令來建立這種連接。要使用該命令，首先需要創建接觸對，并且要對目標-接觸單元的關鍵字進行設置。下面的絕大多數操作都是圍繞該命令進行的。
【求解步驟】
1.前處理
1.1 創建單元
/PREP7
ET,1,SOLID187
ET,2,SHELL181
ET,3,TARGE170
KEYOPT,3,5,1
ET,4,CONTA175
KEYOPT,4,2,2
KEYOPT,4,12,5
上述命令分別定義了4種單元。
第1種是實體單元，第2種是殼單元，他們分別用于建模上述梁的實體部分和空心部分。
第3-4種則是用于模擬接觸部分，就是實體與空心的接觸部分。
這里對于這兩種單元均設置了關鍵字，這些關鍵字的設置是使用后面的命令“SHSD”所必須的。
1.2 創建實常數
R,1,0.02
R,2
R,3
R,4
R,5
這里創建了5個實常數。
第1個實常數用于定義空心梁的厚度
第2-5個實常數分別用于定義4個接觸對。
1.2 創建材料類型
MP,EX,1,2e11
MP,PRXY,1,0.3
上述命令定義了材料的彈性模量和泊松比。
1.3 創建中間的空心梁
/VIEW,1,1,1
BLOCK,-0.14,0.14,-0.14,0.14,0,0.98
VDELE,1,,,0
ADELE,1,2,1,1
上述命令首先創建了一個長方體，
然后刪除了體本身，留下構成長方體的面，線和關鍵點。
最后又刪除了兩端的面。
結果如下圖。
展開 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中級認證窗口課程：LS-DYNA中殼體與實體單元連接技術應用
1.2 合并節點法
合并節點法只傳遞平動自由度，不傳遞轉動自由度
FEM->Element Tools->Duplicate Nodes，點擊Show Dup Nodes，如圖所示，重合節點用黃色方框顯示
點擊Merge Dup Nodes進行合并。
進行計算，可以用LS-RUN，也可以用ANSYS Mechanical APDL Product Launcher，計算結果如下圖所示，從圖中明顯看出殼單元并未被拉開。
修改受力方向為Z向，變形很大，并且板沒有變彎，說明沒有傳遞彎矩。
1.3 約束法
FEM->Model and Part->Create Entity打開Entity Creation窗口選擇Set Data->*SET_NODE，點擊右側的Cre單選按鈕，即新建一個節點組，在Title中輸入名稱為by hand，如下左圖所示：
分別選擇同一位置殼單元和實體單元上的點，選擇殼單元時在左上角關閉實體單元，選擇實體單元時打開殼單元，關閉殼單元，最后恢復選擇。
在Entity Creation中點擊Apply，在下側列表中出現新生成的節點對。依次生成七個節點對，如果想看批量處理法，請查看視頻。最后點擊Done關閉Entity Creation。
在Keyword Manager中選擇Constrained->NODE_SET，雙擊或者選擇后點擊Edit按鈕，打開Keyword Input Form，輸入剛剛創建的NSID，每一對節點分別進行三個平動方向的約束，即DOF分別選擇1、2、3，也就是x、y、z方向。
依次生成七個約束，如果想看批量處理法，請查看視頻。最后點擊Done關閉Entity Creation。
展開 基于ANSYS的實體單元扭矩施加方法總結(原創帖子，轉載請注明出處，謝謝！技術鄰ID有限元中解人生) ￥1
基于ANSYS的實體單元扭矩施加方法總結
1、 引言
在實際工程問題中，扭矩無處不在。如攻絲的絲錐、車床的光桿、攪拌軸、汽車傳動軸等等，均為受扭構件，承受扭矩作用。為了更好的分析上述構件在扭（轉）矩作用下的變形、應力、應變等物理量，現代先進設計制造分析方法引入有限元來模擬結構在外載荷作用下的響應問題。對于很多工程模型，必須考慮結構的一些幾何特征，如軸的鍵槽、絲錐的螺紋面等。因此，實體模型上扭矩的施加就成為一個非常關鍵的問題。這包括扭矩施加的形式、位置，不同方式施加的扭矩會導致整體剛度矩陣的不同，最終會導致應力奇異，影響結果的評定。ANSYS作為全球最通用的大型有限元分析軟件之一，其強大的分析功能已為國內外一致認同，現已成為許多領域結果評定的行業標準。由于ANSYS中不能直接對實體單元施加力矩，傳統方法采用若干對力偶來代替扭矩，該方法容易導致局部應力集中；改進的方法引入一些特殊單元如rbe3單元、mpc184單元、mass21單元等，通過引入這些特殊單元，能夠比較好的實現扭矩的施加，但是特殊單元的引入又改變了整體剛度矩陣。為了解決由于引入特殊單元而導致影響整體剛度矩陣的問題，有學者等提出采用接觸單元能夠很好的解決扭矩的施加問題。
本文旨在綜合關于扭矩施加的各種方法，并對這些方法進行分析比較，從而找到關于實體單元扭矩施加有效、合理的方法，為結構有限元分析提供有益的參考。
2、 ANSYS中扭矩的施加
2.1 工程實例
現以長為0.2m直徑為100mm的實心鋼管為例說明扭矩的施加。鋼管材料視為線彈性，其彈性模量及泊松比分別為：E=2e11Pa，μ=0.3。 鋼管一端固定，另一端受1000N.m扭矩作用。
展開 