ansys 單元生成的案例
ansys模塊化仿真系列文章（一）梁單元截面特性標準生成
開篇點題，不說廢話，直接給出生成梁單元的手動操作方式和模塊化命令流。
手動操作
介紹一下標準化生產梁單元截面特性，便于后續的梁單元建模和仿真。
1，CAD做成sat文件：首先生成面域
2，file導入ACIS
3，定義單元，劃分網格
ET,1,plane82 !添加單元類型plane82
LSEL,all !選擇所有線段
LESIZE,all,10 !設定網格尺寸，根據具體圖形尺寸進行調整
MSHAPE,0,2D !采用四邊形網格單元
MSHKEY,0 !采用自由網格
AMESH,ALL !劃分網格
4，截面寫出-界面操作
section->beam->write
5，截面寫入-界面操作
section->beam->read->plot
模塊化命令流
! 模塊化寫出截面命令流
finish
/clear
/prep7
str1 = 'name'
~SATIN,'name','sat',,SURFACES,0
*get,a_count,area,,count ! 獲得面號
/facet,normal ! 面顯示正常
allsel
ET,1,plane82 !添加單元類型plane82
LSEL,all !選擇所有線段
LESIZE,all,12 !設定網格尺寸，根據具體圖形尺寸進行調整
MSHAPE,0,2D !采用四邊形網格單元
MSHKEY,0 !采用自由網格
AMESH,ALL !
展開 cohesive單元生成插件
4.插件文件夾中還有一個簡單的2D例子（Job-2.inp），驗證此方法生成的cohesive單元的可計算性（可開裂）
課程地址：詳見本人的視頻課程，有詳細講解代碼和python入門
空氣彈簧單元生成方法
彈簧單元生成
從inp文件中，將*Element, type=S4R下面的單元與節點復制到excel表格中。每行有5個數字，第一個數字代表單元編號，后面四個代表這個單元所包含的節點編號。我們要做的就是把單元編號改掉，而節點編號不變。這樣，改動過的新單元與舊單元就擁有了共同的節點。新單元type=F3D4，這樣我們就有了與S4R共節點的但是單元編號不一樣的F3D4空氣彈簧單元。我這里將S4R單元統一移動了100個編號，變成F3D4單元。注意，這兩種單元的編號一定不能重復。如圖2所示。
圖2
將F3D4單元編號和節點復制到INP文件中，并在開頭寫上關鍵字
*Element, type=F3D4
還要增加一個空氣彈簧單元集合，
*Elset, elset=SET-AIRSKIN, generate
101, 152, 1
2. 預定義場（壓力、溫度）
打開cae，重新讀入上面改動過的inp模型。
添加一個與節點不重合的參考點集合：Set-AirRP(0,1,0)
在model-edit keywords里面編輯predified field。加入
** PREDEFINED FIELDS
**
** Name: Predefined Field-InitialTemperature Type: Temperature
*Initial Conditions, type=TEMPERATURE
Part-2-1.Set-AirRP, 300.
*Initial Conditions, TYPE=PORE PRESSURE
Part-2-1.Set-AirRP, 414., 0, 414. , 500
預定義初始溫度為300k，初始內壓為414Pa。
注意，這里要檢查一下shell的方向。
展開 cohesive單元批量生成插件 ￥60
cohesive單元批量生成插件
Samcef field超單元生成
在大型復雜的機構中，為了節省計算資源，保證計算精度，常需要利用超單元對部件進行柔性化。Samcef field中超單元生成的步驟比較簡單，計算完成后也沒有可顯示的計算結果。有些重要的步驟需要考慮，否則生成的超單元難以與其他部件連接。兩個比較容易忽略的操作步驟是：將連接幾何點設置為“not a datum”；將連接點在assembly中設置為“retained nodes”這樣才能與其他部件進行連接。
以直驅風機的底座為例，進行了超單元生成的操作，視頻網址：
http://v.youku.com/v_show/id_XOTQ2MDk4MzI0.html
http://pan.baidu.com/s/1ntOfKTF
展開 單元和網格生成(Zone and Grid Generation)方法
1 引言
數值模擬的第一步是對模型進行單元和網格劃分，選擇正確的網格生成方法以及有效地部署該方法是數值模擬過程中的關鍵步驟。由于網格生成需要考慮很多因素，因此一個項目使用的方法不能照搬到另一個項目。這個筆記簡要討論了FLAC3D(也適用于3DEC)網格劃分的一些注意事項。
2 方法分類
(1) FLAC3D內置了三種方法來劃分單元和網格，它們是：基元生成(primitives)、擠壓生成(extrusions)和積木生成(building blocks)。每種方法都有明顯的優點和缺點。對于網格的初始生成，它們是相互排斥的(即不能同時用基元和擠壓來創建一組區域)。然而，程序提供的命令允許根據需要對來自任何方法的區段進行進一步修改、組合等等。實現模型幾何形狀的正確途徑可能涉及到這些區域生成方法的巧妙組合。
(2) 此外，FLAC3D可以從導入的文件中創建單元，例如ANSYS和Abaqus。也可以使用Griddle或其它第三方工具來生成網格,例如CUBIT, HyperMesh, TrueGrid等來劃分網格。
(3) FLAC3D也有一些附加的工具來生成和修改單元。對于那些幾何形狀復雜，但單元面與表面的精確一致性在物理上并不重要的情況，使用幾何描述和網格密集化來創建八叉(octree )網格，當有復雜的不相交的地層或不均勻的表面地形時，可以使用zone generate from-topography從地形生成區域的命令。
3 劃分策略
沒有哪一種網格生成方法能夠適應所有情況，最佳的方法極大依賴于模型的幾何形狀和分析目標。下面所列的是一個總的策略。
(1) 如果問題具有非常簡單、規則的幾何形狀，或者恰好符合zone create創建命令可用的形狀之一，那么使用基元生成網格理論上是最快的方法。
展開 空氣彈簧單元生成方法
彈簧單元生成
從inp文件中，將*Element, type=S4R下面的單元與節點復制到excel表格中。每行有5個數字，第一個數字代表單元編號，后面四個代表這個單元所包含的節點編號。我們要做的就是把單元編號改掉，而節點編號不變。這樣，改動過的新單元與舊單元就擁有了共同的節點。新單元type=F3D4，這樣我們就有了與S4R共節點的但是單元編號不一樣的F3D4空氣彈簧單元。我這里將S4R單元統一移動了100個編號，變成F3D4單元。注意，這兩種單元的編號一定不能重復。如圖2所示。
圖2
將F3D4單元編號和節點復制到INP文件中，并在開頭寫上關鍵字
*Element, type=F3D4
還要增加一個空氣彈簧單元集合，
*Elset, elset=SET-AIRSKIN, generate
101, 152, 1
2. 預定義場（壓力、溫度）
打開cae，重新讀入上面改動過的inp模型。
添加一個與節點不重合的參考點集合：Set-AirRP(0,1,0)
在model-edit keywords里面編輯predified field。加入
** PREDEFINED FIELDS
**
** Name: Predefined Field-InitialTemperature Type: Temperature
*Initial Conditions, type=TEMPERATURE
Part-2-1.Set-AirRP, 300.
*Initial Conditions, TYPE=PORE PRESSURE
Part-2-1.Set-AirRP, 414., 0, 414. , 500
預定義初始溫度為300k，初始內壓為414Pa。
注意，這里要檢查一下shell的方向。
展開 機器學習 |使用門控循環單元網絡生成文本 ￥2
本文將演示如何通過構建 Gated Recurrent Unit Network 來構建文本生成器。訓練網絡的概念過程是首先向網絡提供網絡正在訓練的文本中存在的每個字符的映射到唯一數字。然后將每個字符熱編碼為一個向量，這是網絡所需的格式。
所述程序的數據是著名詩人的短詩和著名詩集，格式為 .txt。它可以從 kaggle 下載。
samcef 二維轉子模型的超單元生成
二維轉子動力學模型不僅可以用于系統的簡化分析，也可以對其進行超單元生成，然后作為次級部件導入其他系統模型中。需要指出的是，超單元生成的計算，沒有任何可以顯示的三維結果。
通過對本案例的介紹，可以了解到samcef在超單元生成方面的功能，主要包括：
· Howto define 2D rotor geometry in Field (vertex, lines, face)
· Howto define an elastic material
· Howto create a superelement from this 2D geometry
· Howto create a model, and import the superelement in it
· Howdo define a Lumped mass attached to the superelement
· Howdo define a Ground bearing linked to retainednodes of the superelement
· Howdo define locking on translational degrees of freedom of the superelement
· Howto perform a “critical speed” analysis (ROTOR) of the model containing thesuperelement
· Howto examine the results and report them.
How to perform recovery ofthe superelement and see internal modal displacements
附件為本案例的具體操作介紹。
展開 一鍵生成非線性彈簧單元！！！
概述
在使用ABAQUS中的非線性彈簧單元研究鋼筋混凝土粘結滑移、土體和樁的非線性剛度等問題時，需要在樁基和土體間建立彈簧單元。手動操作不太現實，因此本文使用python開發了腳本，可用于快速生成彈簧單元。
2.效果演示
3.核心代碼
給出核心代碼如下供大家參考，如想快速獲取需完整代碼可聯系小編（扣q1871858827）。
4.非線性彈簧單元
ABAQUS/CAE中暫時僅支持定剛度彈簧單元，如需創建非線性彈簧單元，需要在inp文件中修改關鍵字。
inp文件中修改示意（僅供參考）。
5.參考文獻
Abaqus Example Problems Guide (6.14)
Abaqus Analysis User's Guide (6.14)
展開 自研批量生成0厚度內聚力單元插件 ￥400
該插件基于python語言批量生成內聚力單元，使用該插件需要配置python語言環境。使用者可隨時反饋使用過程中遇到的任何問題，發送至郵箱shenz1hao@126.com，承諾修正一切問題。
插件功能具體包含：
（1）三種插入方式：兩種材料內部+邊界、僅邊界、邊界+一種材料內部
（2）任意單元形式：三角形單元、四邊形單元、六面形單元、楔形單元、四面體單元以及這些單元之間的任意組合
（3）單元數量及INP生成方式無限制
展開 
基于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、，先跑一遍，看看到底徐變是怎么個事兒。
展開 samcef wind turbines 的超單元生成視頻
這里上傳了SWT高精度建模時，超單元模型的建模方法視頻。主要是在前處理軟件samcef filed中完成的。
步驟文檔見附件。優酷視頻：http://v.youku.com/v_show/id_XNzMwMTU4MDE2.html
百度網盤中的高清視頻：http://pan.baidu.com/s/1o6FNSaI
Samcef超單元創建.pdf
Abaqus前處理插件-生成Voronoi單元集 ￥100
通過本款插件可方便生成類Voronoi的單元集，二維圖形和三維結構通用，適于任意幾何構型。
閉源程序保證支持2016及以下版本，高版本建議采用開源程序：
https://www.yqgqt.org.cn/content/post/f6b5ef4e-5878-43cb-bb9a-9cfd4dcb05b8
更新日志：
2020.7.19 ： 改進對早于2016的舊版本兼容性。
Abaqus前處理插件-生成Voronoi單元集(開源程序) ￥200
通過本款插件可方便生成類Voronoi的單元集，通用于二維圖形和三維結構，適于任意幾何構型。
包含程序的開源 .py 文件，可直接作為腳本文件執行，相信能夠幫你打開Abaqus二次開發的思路。
更新日志：
2020.7.19 ： 改進對早于2016的舊版本兼容性；
2020.9.15 ：增加隨機程度控制參數 (0~1)。
二維實例，隨機程度 = 0：
三維實例，隨機程度 = 0：
二維實例，隨機程度 = 1：
三維實例，隨機程度 = 1：
程序界面：
如有需求可私信詳聊。
