abaqus平面接觸的案例
平面應變單元CPE4R齒輪傳動接觸應力計算 ￥49.9
厚齒輪的應力符合平面應變狀態，可以采用平面應變單元CPE4R來進行快速接觸應力計算。
在sketch模塊建立非對稱結構齒輪的草圖，然后建立part，并在assembly模塊進行裝配。
非對稱齒輪草圖
齒輪裝配體
通過適當的結構設計，非對稱齒輪可以在定速轉動的情況下獲得按某規律的變化轉速，在工程上經常會用到。
非對稱齒輪傳動分析結果
非對稱齒輪應力云圖
非對稱齒輪齒合區域局部應力云圖
Abaqus+PyQt+Python平面變形歐拉角計算
5 相對歐拉角計算
利用計算絕對歐拉角時得到的坐標系文件，計算平面變形相對歐拉角，如下圖所示，計算平面2相當于平面1、平面3相對與平面1的相對歐拉角。
計算結果如下圖所示。
6 小結
上述軟件用的算法申請了發明專利，軟件申請了軟著。CAE工程師，也可以自制軟件工具，解決重復性、復雜性數據處理等工作痛點。
2021年8月24日于西昌衛星發射中心
SolidWorks平面模型導入ABAQUS建立軸對稱模型
（3）有了螺栓截面的草圖，接下來應用曲面工具中的平面工具按鈕，為螺栓零件區域建立截面模型。
圖5
如圖5所示，利用平面工具，根據草圖2生成了螺栓零件的截面模型，這時，在曲面實體下有了相應的截面實體列表。
圖6
重復新建草圖→平面工具過程完成全部6個零件截面建模，結果如圖6。最后可以隱藏草圖1，使得圖形區的圖面顯得較為干凈。
然后可以另存為Parasolid格式的文件，以供ABAQUS導入使用。
（4）如圖7所示，在ABAQUS中作為裝配導入Parasolid文件。在ABAQUS中自動創建了6個零件實例，這樣就可以為每個零件實例劃分網格和賦予材料、建立零件之間的接觸關系，然后加載分析。
ABAQUS導入的面模型默認是在三維空間中，為了分析軸對稱模型，需要回到部件位置對每一個零件編輯，改為軸對稱模型。如圖7所示。
圖7
為了在螺栓上施加預緊力，需要在螺栓桿部適當位置進行一次切分。如圖8所示。
圖8
后面在ABAQUS中的操作都是ABAQUS使用者所熟悉的（賦予材料、建立接觸、添加約束、添加螺栓預緊力等），完善模型后進行分析，結果如圖9所示。
圖9
上述過程還是比較簡單的，ABAQUS使用者有的可能不熟悉SolidWorks的草圖繪制和特征工具的操作，SolidWorks是公認學習曲線非常平緩的軟件，簡單的摸索就能用起來。需要注意的是：要找到SolidWorks里的曲面工具欄；在ABAQUS中導入時注意，要進行接觸分析需要從裝配位置右鍵導入；還需注意默認導入時三維空間（的曲面），要進行平面或者軸對稱分析，需要回到部件位置對每一個部件修改為二維平面或軸對稱，以使得模型的空間維度是正確的。
展開 ABAQUS UEL - 損傷材料本構簡單應用于4節點平面單元 ￥300
利用ABAQUS自定義單元子程序，既可以開發新的單元，同時也可以定義新的材料本構模型。本文以損傷模型簡單應用于4節點平面單元為案例，介紹ABAQUS UEL的開發和使用。
如上圖所示，該單元包含4個節點，每個節點有兩個自由度，分別在水平（X）和垂直（Y）方向運動。節點1的兩個自由度被固定，節點4的水平自由度被固定，節點2的垂直自由度被固定。節點3和節點4在垂直方向上向上運動，位移為0.1mm。該正方形單元的邊長為100mm。在input文件里，坐標表示為，
定義節點組合與邊界條件為，
為了讓模型收斂性更好，采用quasi-newton 求解器。時間步設置為，
在文件夾中通過Powershell提交job和子程序，
單個單元的變形為，
采用不同的 ??
，在后處理中得到損傷因子的變化，
相對應的力-時間關系為，
對于多個單元的情況，比如9單元組成的模型，
具體介紹見知乎：ABAQUS UEL - 損傷材料本構簡單應用于4節點平面單元 - 知乎 (zhihu.com)
相對應的UEL代碼和input文件在付費內容中，
展開 
Abaqus帶螺紋螺栓接觸應力分析淺析 Abaqus帶螺紋螺栓接觸應力分析淺析
目前的常規做法通常有兩種：1.簡化，用RBE2和beam梁來代替螺栓，這樣不能反映連接螺栓真實應力，圖1為某結構連接螺栓簡化的beam梁應力云圖，沒有接觸應力：
.直接做出來螺栓螺紋采用接觸分析，雖然得出的結果很精確，但這樣前處理工作量大（螺栓和螺紋用六面體網格建模）、計算量大（接觸收斂困難），如圖為某結構帶螺紋螺栓和連接件模型（圖2）和計算得出的結果（圖3）：
圖3 計算結果
那么，有什么好辦法可以不用簡化帶螺紋螺栓，不用直接做出帶螺紋螺栓，又能得到足夠精確的結果？
運用大型通用非線性有限元分析軟件Abaqus，只需要在接觸定義中設置跟實際螺紋形狀有關聯的參數，如牙角、螺距、螺栓小徑等，就可以模擬真實的連接螺栓接觸狀況。既可以得到足夠精確的分析結果，又節省了時間專注進行其他的分析設置。如圖4，為連接螺栓接觸來定義帶螺紋螺栓：
圖4 連接螺栓接觸來定義帶螺紋螺栓
圖5為某結構直徑10MM的帶螺紋的連接螺栓接觸壓力分布云圖：
圖5 某結構直徑10mm帶螺紋的連接螺栓接觸壓力分部云圖
展開 ABAQUS UEL-梯度損傷模型應用于4節點平面單元 ￥600
本文詳細介紹了如何將梯度損傷模型應用于4節點平面單元，并在有限元模型中進行模擬。
ABAQUS提供了UEL（user defined element）給使用者進行開發。筆者利用UEL開發4節點平面單元，其邊界條件如下圖所示。其中，節點1的X、Y方向被限制住，節點2的Y方向被限制，節點4的X方向被限制，節點3、4的Y方向有豎向位移0.1mm。單元為100*100mm的二維正方形。
每個節點除了X和Y方向的位移，還帶有非局部應變（nonlocal strain）。
單個單元模型，
多個單元模型，
具體內容可參見知乎文章：
ABAQUS UEL-梯度損傷模型應用于4節點平面單元 - 知乎 (zhihu.com)
相應的input文件和uel代碼付費可見，
展開 abaqus中平面應力應變厚度對切削力的影響 ￥5
在鋁合金的二位正交切削仿真中，不同的平面應力應變厚度的對切削力的影響結果
以上為設定值為1的情況
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中平面應力單元高斯積分點的順序為：
展開 基于abaqus的鋼筋混凝土平面框架倒塌性能分析 ￥100
<p>結構在遭遇偶然突發事件后， 不可避免的會導致結構局部破壞或者損傷， 如果剩余結構不能有效的承擔結構初始破壞和損傷造成的內力變化， 剩余結構就會發生進一步破壞， 造成多米諾骨牌式的連鎖反應，從而造成大范圍嚴重破壞乃至倒塌，這就是通常所說的連續倒塌。附件中只有一個cae有限元模型。</p><p><img src="https://img.jishulink.com/upload/201905/597ca43812cb414e98ab1fd96e276a82.jpg" alt="2019-05-07_105121.jpg"></p><p><img src="https://img.jishulink.com/upload/201905/77bbfce9cdc84985b02d22088e6933bc.jpg" alt="2019-05-07_105131.jpg"></p><p><img src="https://img.jishulink.com/upload/201905/5dd8c1f0916b47cab1cbf1df2992706f.jpg" alt="2019-05-07_105149.jpg"></p><p><img src="https://img.jishulink.com/upload/201905/984ad184d1ba41209ee701b4d0aec1de.jpg" alt="2019-05-07_105208.jpg"></p><p><br></p><p><img src="https://img.jishulink.com/upload/201905/48d8ca94aa6e42768f67ad19803e150b.jpg" alt="2019-05-07_105234.jpg"></p><p><img src="https://img.jishulink.com
展開 abaqus2020-三維-顯示分析-通用接觸或接觸對接觸-單元刪除法模擬裂紋，單元穿透問題！！
1 abaqus2020-三維-顯示分析-單元刪除法模擬裂紋，僅采用通用接觸時，模型中出現明顯穿透，結果不合理！
2 abaqus2020-三維-顯示分析-單元刪除法模擬裂紋，僅采用接觸對接觸時，模型中出現少許穿透，結果相對合理，但不是最理想狀態！
3 abaqus2020-三維-顯示分析-單元刪除法模擬裂紋，同時采用通用接觸+接觸對接觸時，模型中無明顯穿透，結果合理！
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題
摘要：
采用基于ABAQUS的UEL子程序開發4節點平面應變等參單元，采用雙線性形函數，4點高斯積分，本構關系為線彈性各向同性材料，得到的單元剛度矩陣和ABABUS自帶的CPE4單元的單元剛度矩陣（剛度矩陣輸出方式為*element matrix output, elset= ALLE, stiffness=yes, OUTPUT FILE=USER DEFINED）不同；對比ANSYS的單元剛度矩陣，結果顯示兩者也不相同。問題出在哪里呢？本文檔將對此問題進行回答。
本文可以作為ABAQUS高級子程序UEL的入門級教程，做UEL的應該關注下！
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題（技術鄰 藍牙）.pdf
展開 
【接觸分析】詳述ABAQUS接觸分析（1/2）
基于節點（node-based）的表面：定義了在節點集和表面之間的接觸，可以應用并取得與基于邊界的表面同樣的效果。
2.3 剛性表面
剛性體的表面，可以定義為一個解析形狀，或者是基于與剛體相關的單元的表面。
解析剛性體表面的三種基本形式：在二維中：一個二維的分段剛性表面，可以在二維平面上應用直線、圓弧和拋物線弧定義表面的橫截面；在三維中：在用戶指定的平面上應用對于二維問題相同的方式定義，然后由這個橫截面繞一個軸掃掠形成一個旋轉面，或沿一個矢量拉伸形成一個場的三維表面。
優點：只用少量的幾何點便可以定義、計算效率很高
缺點：在三維情況下，能夠創建的形狀范圍有限
離散形式的剛性表面：基于構成剛性體的單元面，可以創建比解析剛性表面幾何上更為復雜的剛性面，定義方法與定義可變形體表面的方法完全相同。
在ABAQUS/Explicit中，解析剛性表面只能應用于接觸對算法。
三、接觸面間的相互作用
接觸面間的相互作用包含兩部分：接觸面間的 法向作用、接觸面間的 切向作用。切向作用包括：接觸面間的 相對運動（滑動）、可能存在的 摩擦剪應力。ABAQUS中接觸相互作用的默認模型： 沒有粘結的無摩擦模型
3.1 接觸面的法向行為
間隙（clearance）：兩個表面分開的距離。
當兩個表面間的間隙為零時，在ABAQUS中施加了接觸約束。在接觸問題的公式中，對接觸面之間能夠傳遞的 接觸壓力的量值未作任何限制。當接觸面之間的接觸壓力變為零或負值時，兩個接觸面分離，并且約束被移開。這種行為代表了“硬”接觸。
當接觸條件從“開”（間隙值為正）到“閉”（間隙值為零）時，接觸壓力會發生劇烈的變化，有時可能會使得在ABAQUS/Standard中的接觸模擬難以完成。但在ABAQUS/Explicit中不會這樣，因為其顯式算法不需要迭代。
展開 【接觸分析】詳述ABAQUS接觸分析（2/2）
八、ABAQUS/Standard和ABAQUS/Explicit的比較
在ABAQUS/Standard和ABAQUS/Explicit中的力學接觸算法具有本質區別，體現在如何定義接觸條件，主要區別如下：
· ABAQUS/Standard在施加接觸約束時應用
嚴格的主從權重，約束從屬表面的節點不能侵入主控表面，而主控表面上的節點原則上可以侵入從屬表面。ABAQUS/Explicit包括這個公式，但是典型地它默認應用
平衡主從權重。
· ABAQUS/Standard和ABAQUS/Explicit都提供了
有限滑動接觸公式，但是，在ABAQUS/Standard中的二維有限滑動公式要求主控表面是光滑的，而在ABAQUS/Explicit的主控表面是由面元構成的，除非是光滑的解析剛性表面。
· ABAQUS/Standard和ABAQUS/Explicit都提供了
小滑移接觸公式，但是在ABAQUS/Standard中的小滑移公式根據從屬節點的當前位置向主控節點傳遞載荷，ABAQUS/Explicit總是通過固定點（anchor point）傳遞載荷。
· ABAQUS/Explicit在接觸邏輯中可以考慮殼和膜的當前厚度和中面偏移，而ABAQUS/Standard不能夠做到。
· ABAQUS/Explicit
通用接觸算法的許多優勢在ABAQUS/Standard中是不具備的。
由于以上差異，所以在一個ABAQUS/Standard分析中定義的接觸不能導入一個ABAQUS/Explicit分析中，反之亦然。
九、小結
· 接觸分析需要一個謹慎的、邏輯的方法。
展開 【技巧】abaqus輸出通用接觸的某個面的接觸力
INP關鍵字
*OUTPUT, HISTORY, TIME INTERVAL = 0.1 ##0.1為輸出頻率，如計算時長為1s，需要輸出10步
*INTEGRATED OUTPUT, SURFACE =FACE_NAME ##FACE_NAME是通用接觸設置中要輸出的接觸面的名稱
SOF ##輸出面接觸力
Abaqus中通用接觸(General contact)和接觸對（Contact Pairs）的區別。
對于大多數的接觸問題，在ABAQUS中有通用接觸（General Contact）和接觸對（Contact Pair）兩種算法處理，它們的異同主要體現在用戶交互、默認設置、可選設置三個方面。
總的來說，通用接觸算法的相互作用主體、接觸屬性、接觸面屬性是可以各自獨立地指定，它提供了一個更有彈性的方法去增加模型中接觸的細節。通用接觸算法允許非常自動化的接觸定義，盡管也可以采用傳統的、類似于接觸對算法的方法去交互式定義。對于傳統的接觸對算法，相對于全部包括式的自接觸（Self-contact），接觸對算法的計算效率可能更高，而且使用CAE也能比較方便地建立接觸對。因而這兩種接觸算法的選擇其實就是一個在接觸定義的便利性和計算效率性之間的平衡，它們之間的差異主要有：
一、通用接觸（General Contact）和接觸對（Contact Pair）的默認設置差異
1、接觸離散方式：通用接觸算法使用有限滑動和面對面的離散方式，而接觸對算法使用有限滑動和點對面的離散方式；
2、對殼的厚度和偏移的處理：通用接觸算法自動考慮，接觸對算法在使用點對面的離散方式時不考慮殼的厚度和偏移；
3、接觸的執行：通用接觸算法采用罰函數方法，接觸對算法在使用點對面的離散方式時采用拉格朗日乘數方法；
4、初始過盈量的處理：通用接觸算法采用無應變調整的方法消除過盈量，接觸對算法將過盈量作為穿透在第一個分析增量步處理；
5、主從面指定：通用接觸算法自動指定，接觸對算法必須由用戶指定。
當接觸對算法采用有限滑動和面對面的離散方式時，就沒有前三個差異了。
展開 