ABAQUS創(chuàng)建點單元的案例
基于tcl語言實現(xiàn)單元中心點創(chuàng)建坐標系 ￥15
<p class="ql-align-justify"> 本案例是基于tcl語言實現(xiàn)用戶自定義的單元，并獲取單元的中心點，并依據(jù)單元中心點及單元節(jié)點創(chuàng)建坐標系。具體實現(xiàn)過程見本案例的程序部分。</p><p class="ql-align-justify"> 詳情見收費的程序部分，凡購買本案例的朋友針對該案例有疑問，可私信，謝謝！</p><p class="ql-align-justify"><br></p><p><br></p>
展開 abaqus二次開發(fā)-參考點的創(chuàng)建 ￥20
<p>選擇coord_name,實現(xiàn)有name的坐標創(chuàng)建，選擇coord_name_non,參考點的名稱系統(tǒng)自動創(chuàng)建</p>
abaqus二次開發(fā)實例1:多個RP點一次創(chuàng)建
背景:在abaqus如果需要創(chuàng)建多個RP點，通常我們需要來回切換復制坐標黏貼并創(chuàng)建，較為重復性，
插件功能: 多個RP坐標一次創(chuàng)建，可以提前自定義RP點名稱（第一列）或者是否copy（第五列），如果這兩列都沒設置，按系統(tǒng)默認類型創(chuàng)建，同時也會創(chuàng)建相對于的set集，set名稱為自定義名稱或系統(tǒng)。
創(chuàng)建的方式:基于abaqus的Rsg模塊提前定義用戶界面，如圖所示，這個較為簡單，保存這個用戶界面會生成*DB.py和*_plugins.py兩個文件，這里不需要編輯這兩個文件。 內(nèi)核程序:需要設置函數(shù)def用于在rsg界面調用。
代碼只有100行以內(nèi)，適合初學者學習使用，以下為內(nèi)核程序，界面可以在RSG自己設置，選擇內(nèi)核程序的函數(shù)即可
createcoordinates.rar
展開 hypermesh_abaqus中fastener焊點單元和襯套BUSH單元創(chuàng)建流程 ￥1
hypermesh_abaqus中fastener焊點單元和襯套BUSH單元創(chuàng)建流程
HyperWorks（Hypermesh）+Abaqus彈簧單元（spring）創(chuàng)建及設置方法 ￥9.9
<p> 彈簧單元有3種類型：接地彈簧（spring1）、兩結點彈簧（spring2）、軸向彈簧（springA）。</p><p> <strong>spring1</strong>，接地彈簧，一個結點在大地上，只需定義另一個結點；需要定義彈簧力的方向。</p><p> <strong> spring2</strong>，兩結點彈簧，需要定義彈簧力的方向。</p><p> <strong> springA</strong>，軸向彈簧，不需要定義彈簧力的方向，由兩結點的連線方向確定。</p><p> 常使用springA彈簧單元。</p><p><br></p><p> 本案例分別介紹<u>HyperWorks（Hypermesh）</u>和<u>Abaqus</u>中彈簧單元springA和spring2的創(chuàng)建及設置方法（spring1可參照設置）。</p><p> 以圖文方式詳細描述每一步需要填寫的內(nèi)容及釋義，通過本使用教程，您將可以按教程詳細步驟一步步設置彈簧單元spring，即便是小白也能快速上手使用。</p><p><br></p><p><br></p><p><br></p>
展開 abaqus-python 利用getByBoundingCylinder(...)創(chuàng)建單元集合
在學習的時候被網(wǎng)上的帖子和幫助文檔誤導（其實是新新手的緣故），命令老寫成這個樣子：
elist = e.getByBoundingCylinder(7,6,0,7,6,20,5.8)
# 正確命令：
elist = e.getByBoundingCylinder(center1=(7,6,0),center2=(7,6,20),radius=5.8)
abaqus python 二次開發(fā)攻略 P195,P196
abaqus 用戶幫助文檔
2 簡單案例
案例簡介：選取多個圓柱體框中的單元并創(chuàng)建為集合。
具體命令：
################################
# 可以運行#######
from abaqus import*
from abaqusConstants import*
p = mdb.models['Model-3'].parts['Part-1']
e = p.elements
elist_I = e.getByBoundingCylinder(center1=(7,6,0),center2=(7,6,20),radius=5.8)
elist = elist_I
elist_I = e.getByBoundingCylinder(center1=(21,6,0),center2=(21,6,20),radius=5.8)
elist = elist+elist_I
p.Set(elements=elist,name='Set-matrix')
print(type(elist)) # <type 'Sequence'>
這就是運行結果。
展開 Abaqus創(chuàng)建零厚度cohesive單元
建立cohesive層的方法主要有：
方法一、建立完整的結構，然后在上面切割出一個薄層來模擬cohesive單元，用這種方法建立的cohesive單元與其他單元公用節(jié)點，并以此傳遞力和位移。
方法二、分別建立cohesive層和其他結構部件的實體模型，通過“tie”綁定約束，使得cohesive單元兩側的單元位移和應力協(xié)調。
有些情況下需要建立零厚度的cohesive單元以進行裂紋擴展的計算模擬，下面以第一種方法進行零厚度的cohesive單元的建模。
新建一個二維模型，如下圖所示：
假設在模型的中間有一層零厚度的cohesive單元層，做一個切割處理，如下圖所示：
切割出一個很小的縫隙之后進行網(wǎng)格劃分，如下所示：
完成這一步之后需要通過菜單欄Mesh-creat mesh part將單元變成孤立網(wǎng)格。
在此基礎上即可創(chuàng)建零厚度的單元。點擊菜單欄Mesh下面的Edit，如下所示：
然后框選出中間縫隙的上面一層節(jié)點，之后再選擇下面一層的一個單元面，將縫隙的上面的所有節(jié)點投影到縫隙的下層。
投影完的效果如下，中間的縫隙不見了。
打開節(jié)點號的顯示，如下所示，可以發(fā)現(xiàn)在中間的一層位置處相同位置有兩個節(jié)點，也即是該層單元為零厚度單元
再之后，通過Assign Element Type可修改中間這一層重合的單元，修改為cohesive單元。
長安CAE的博客
http://blog.sina.com.cn/zuoerninan
展開 基于ABAQUS彈簧單元（springa+spring1）的創(chuàng)建及靜力學分析 ￥40
彈簧單元是ABAQUS的特殊功能單元，可以直接定義結點受到的力與結點（相對）位移的關系，所以彈簧單元可以方便的表達一些界面接觸特性。彈簧的種類：按照本構分，彈簧單元可分為線性（linear）彈簧和非線性（nonlinear）彈簧。按照彈簧力的方向及彈簧幾何特點，又可以分為接地彈簧（spring1），兩結點彈簧（spring2），軸向彈簧（springa）。
線性（linear）彈簧可以通過CAE定義，非線性（nonlinear）彈簧則需要修改或者寫inp 文件。spring1，接地彈簧，該彈簧的一個結點（隱含的，不需要定義）是不動的，另一個結點定義在我們需要約束的節(jié)點上。彈簧力的方向，即被約束自由度方向，需要我們定義，既可以在整體坐標系下定義，又可以在結點局部坐標系下定義（查看orientation）。
spring2，兩結點彈簧，彈簧力的方向同上。
springa，軸向彈簧，彈簧力的方向由兩結點的連線方向確定。
另外注意：spring1，spring2可以約束轉角自由度，即抗扭彈簧，而springa不可以。
本案例主要講述如何在ABAQUS中創(chuàng)建接地彈簧（spring1）、軸向彈簧（springa）。購買本案例的朋友附件中同時贈送了非線性接地彈簧的創(chuàng)建模型。
彈簧變形動畫
ABAQUS中spinga彈簧及接地彈簧的創(chuàng)建
ABAQUS中有限元分析結果
展開 基于Hyperworks+ABAQUS彈簧單元（springa）的創(chuàng)建及靜力學分析 ￥35
彈簧單元是ABAQUS的特殊功能單元，可以直接定義結點受到的力與結點（相對）位移的關系，所以彈簧單元可以方便的表達一些界面接觸特性。彈簧的種類：按照本構分，彈簧單元可分為線性（linear）彈簧和非線性（nonlinear）彈簧。按照彈簧力的方向及彈簧幾何特點，又可以分為接地彈簧（spring1），兩結點彈簧（spring2），軸向彈簧（springa）。
線性（linear）彈簧可以通過CAE定義，非線性（nonlinear）彈簧則需要修改或者寫inp 文件。spring1，接地彈簧，該彈簧的一個結點（隱含的，不需要定義）是不動的，另一個結點定義在我們需要約束的節(jié)點上。彈簧力的方向，即被約束自由度方向，需要我們定義，既可以在整體坐標系下定義，又可以在結點局部坐標系下定義（查看orientation）。
spring2，兩結點彈簧，彈簧力的方向同上。
springa，軸向彈簧，彈簧力的方向由兩結點的連線方向確定。
另外注意：spring1，spring2可以約束轉角自由度，即抗扭彈簧，而springa不可以。
本案例講述的是如何在Hyperworks的ABAQUS模塊中創(chuàng)建springa（軸向彈簧）。
彈簧變形動畫
ABAQUS中有限元分析結果
展開 Abaqus 中創(chuàng)建零厚度cohesive單元的幾種方法
基于0厚度cohesive單元的冰雹隨機開裂分析
本帖主要介紹四種咋abaqus 中創(chuàng)建零厚度cohesive單元的方法，四種方法的詳細說明可參見技術鄰學院教學視頻，鏈接如下：
http://www.yqgqt.org.cn/college/video/c10020
現(xiàn)列舉四種方法如下：
1.網(wǎng)格節(jié)點偏移
先生成孤立網(wǎng)格，在Edit mesh工具中，選擇node →edit ，將需要偏置的節(jié)點按照偏移距離進行設置即可，如下圖所示。
2.網(wǎng)格節(jié)點投影
先生成孤立網(wǎng)格，在Edit mesh工具中，選擇node →project ，將需要投影的節(jié)點投影到指定的單元面上。如下圖所示
以上兩種方法都是基于孤立網(wǎng)格進行操作，適合簡單的平板結構。
3.高版本Abaqus直接生成
如果你有高版本的abaqus 比如2016版，那么很幸運，這個版本自帶了插入cohesive的工具，并且同時支持基于幾何的網(wǎng)格和孤立網(wǎng)格。該功能同樣位于Edit mesh工具欄，mesh→insert cohesive seams按鍵，如下圖所示。
該功能模塊，較為強大，能對復雜，曲面結構插入cohesive，但是不能實現(xiàn)任意兩個單元面之間插入cohesive。
4.自定義腳本
如果你想模擬隨你開裂，想在任意兩個單元之間插入零厚度的cohesive,只能自己編寫腳本，去修改模型CAE或者修改inp文件，現(xiàn)在提供一下這種方法的思路，如下：
即先將原單元離散，然后提取原相鄰單元共用面上的節(jié)點，復制該節(jié)點，進行網(wǎng)格重構。以一個實心圓球為例，下圖左為六面體單元組成的網(wǎng)格模型，右圖為插入的零厚度cohesive單元，該方法可以實現(xiàn)任意結構任意單元面之間插入0厚度的cohesive單元。
展開 基于Hyperworks+ABAQUS彈簧單元（spring2）的創(chuàng)建及靜力學分析 ￥35
彈簧單元是ABAQUS的特殊功能單元，可以直接定義結點受到的力與結點（相對）位移的關系，所以彈簧單元可以方便的表達一些界面接觸特性。彈簧的種類：按照本構分，彈簧單元可分為線性（linear）彈簧和非線性（nonlinear）彈簧。按照彈簧力的方向及彈簧幾何特點，又可以分為接地彈簧（spring1），兩結點彈簧（spring2），軸向彈簧（springa）。
線性（linear）彈簧可以通過CAE定義，非線性（nonlinear）彈簧則需要修改或者寫inp 文件。spring1，接地彈簧，該彈簧的一個結點（隱含的，不需要定義）是不動的，另一個結點定義在我們需要約束的節(jié)點上。彈簧力的方向，即被約束自由度方向，需要我們定義，既可以在整體坐標系下定義，又可以在結點局部坐標系下定義（查看orientation）。
spring2，兩結點彈簧，彈簧力的方向同上。
springa，軸向彈簧，彈簧力的方向由兩結點的連線方向確定。
另外注意：spring1，spring2可以約束轉角自由度，即抗扭彈簧，而springa不可以。
本案例講述的是如何在Hyperworks的ABAQUS模塊中創(chuàng)建spring2（兩結點彈簧，彈簧力的方向同上。），后面有時間將陸續(xù)在后續(xù)案例中講述如何在Hyperworks+ABAQUS中創(chuàng)建spring1（接地彈簧）及springa（軸向彈簧）。
彈簧變形動畫
ABAQUS中有限元分析結果
展開 
ABAQUS輸出單元積分點坐標
方法
在ABAQUS CAE的場輸出中選擇的坐標點是節(jié)點的坐標，而節(jié)點是從積分點插值出來的，單元積分點的信息相對真實。所以最好是獲取積分點的信息，其中積分點的坐標無法在CAE中獲取，需要在關鍵字中添加。具體在每個分析步的單元輸出下面添加COORD，如果需要輸出節(jié)點的坐標也可以在節(jié)點場輸出下面添加COORD（這和CAE中場輸出選擇節(jié)點坐標的效果是一致的）。具體如下圖：
2.注意
在ODB結果中創(chuàng)建場輸出時會附帶著一份XYZ坐標，這個應該也可以當做單元的坐標，，但是我比較過這個附帶的坐標和單元的COORD輸出的坐標，有時候有點差別，可能是數(shù)據(jù)精度的問題。
展開 abaqus怎么獲取單元積分點上的坐標呢
如題 為什么我查詢的時候坐標都是0
abaqus C3D8 單元 計算中采用了多少個積分點？
按照正常的理解，毫無.疑問，abaqus 全積分一定是采用了2x2x2=8個積分點。
從后處理結果來看，似乎也是如此，每個單元存在8個積分點。
然而，如果自己動手跑一遍程序，就會發(fā)現(xiàn)事實遠非如此，采用全積分計算得到的結果與abaqus 存在差異，原因何在？
事實賞，abaqus C3D8 采用的選擇積分方式（selective intergation schema），即對于偏應變，采用8個積分，對于球應變，采用中心點積分。這樣計算得到的結果才能與abaqus 完全對標，亦可從abaqus 幫助文檔得到答案。
展開 Abaqus中平面應力單元高斯積分點的順序
可以輸出umat接口中的變量coords進行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結果為：
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
coords = 5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 3
coords = -5.77350269E-01 5.77350269E-01 1.00000000E-02
npt = 4
coords = 5.77350269E-01 5.77350269E-01 1.00000000E-02
因此Abaqus中平面應力單元高斯積分點的順序為：
展開 