abaqus接觸剛度
關注創建者：王靖雯 創建時間：2023-02-27
abaqus接觸剛度的視頻教程
詳解abaqus接觸分析之二：abaqus中基于面的接觸詳
詳細講解了abaqus中基于面的接觸行為： 接觸模型 接觸離散化（N-S S-S的接觸原理，優劣勢，選用建議） 接觸施加方法（罰剛度、拉格朗日法的使用方法和注意事項，已經如何通過他們來控制穿透和收斂） 體間相互滑移（有限滑移和小滑移的原理和注意事項） 等。
￥25 1小時10分鐘 2143播放
查看
abaqus接觸剛度的實例教程
實體接觸（Contact Entity）
這一接觸類型用于可變形（物體的）節點和剛性幾何面之間，在計算過程中用到了定義幾何面的解析方程。與通常的需要劃分成網格的面段相比，這是一種高級的算法。用于抵抗穿透的接觸力通過罰函數法來獲得。這一方法一般用來在LS-DYNA中模擬那些通常由規則曲面近似而來的假人模型。為了清楚地可視化計算結果，程序會自動在剛性曲面上劃分網格，但是這些網格并不參與計算。解析剛性曲面可以是如下幾種：
無限或有限的平面；
球面；
圓柱面；
超橢球面；
圓環面；
加載曲線中定義的線；
CAL3D/MADYMO 平面；
CAL3D/MADYMO 橢球面；
來自VDA文件的曲面；
來自IGES文件的曲面。
接觸剛度的計算
從節點和主面段之間的接觸關系由兩者之間的線性彈簧表示，彈簧的剛度決定了將施加給從節點和主面段的接觸力。下面簡述兩種計算接觸剛度的方法：
1.罰函數法(SOFT=0 in Optional Card A in '*CONTACT_')
接觸面剛度計算公式如下：
這是默認的計算方法，其中利用接觸面的尺寸和材料的力學性能計算剛度，因此當接觸面之間的材料剛度參數為同一數量級時計算效率很高。當參與接觸的兩個部件的材料不同時，接觸剛度大概是兩個材料剛度中的較小值，如果接觸剛度太小的話則會計算出現錯誤。這種錯誤一般出現在很軟的泡沫材料和金屬材料之間的接觸上。因此我們不建議使用SOFT = 0選項，除非用戶之前用這個選項計算時并沒有出現錯誤。
展開 實體接觸
接觸剛度的計算
從節點和主面段之間的接觸關系由兩者之間的線性彈簧表示，彈簧的剛度決定了將施加給從節點和主面段的接觸力。下面簡述兩種計算接觸剛度的方法：
1.罰函數法(SOFT=0 in Optional Card A in ‘*CONTACT_’)
接觸面剛度計算公式如下：
接觸面剛度計算
這是默認的計算方法，其中利用接觸面的尺寸和材料的力學性能計算剛度，因此當接觸面之間的材料剛度參數為同一數量級時計算效率很高。當參與接觸的兩個部件的材料不同時，接觸剛度大概是兩個材料剛度中的較小值，如果接觸剛度太小的話則會計算出現錯誤。這種錯誤一般出現在很軟的泡沫材料和金屬材料之間的接觸上。因此我們不建議使用SOFT = 0選項，除非用戶之前用這個選項計算時并沒有出現錯誤。
2.軟約束法 (SOFT=1 & 2 on Optional Card A in *CONTACT_ )
這種方法不是默認的，通過參與接觸的節點質量和全局時間步長度來計算接觸剛度。這樣計算出來的剛度與材料參數無關，因此適用于不同材料之間的接觸。如下所示，剛度由節點質量除以時間步的平方并乘以一個確保計算穩定的縮放系數而得。
剛度的計算
通常情況下，對于金屬和金屬之間的接觸，不論是SOFT = 0還是SOFT = 1，計算的結果都是相近的。對于軟泡沫和金屬之間的接觸，SOFT = 1選項會給出高出1-2個數量級的接觸剛度。我們推薦在不同材料的接觸之間使用SOFT = 1選項。
和SOFT = 1類似，SOFT = 2利用了基于節點質量和時間步的罰函數剛度。SOFT = 2調用了基于面段的接觸算法，這種算法起源于Belytschko和他的同事提出的彈球接觸（Pinball contact）。在這種算法下，程序會計算面段之間的接觸問題，而不是傳統的節點-面段接觸問題。
展開 劃分網格
定義邊界條件
5.2 預應力模態
預應力模態，按我的理解就是，假設t=0時，結構的剛度矩陣、質量矩陣為M0,K0；t=t1時，結構（分析對象）收到外部激勵的作用，使得結構的剛度矩陣，質量矩陣發生改變，結構的剛度矩陣、質量矩陣為M1,K1；t=t2時，開始進行特征值提取，此時求解的是t1狀態的結果。
從上面這個理解出發，ABAQUS預應力模態只要在frequency分析步之前進行General,Static分析步,打開NLGeom選項(分析過程中剛度矩陣會不斷變化)。
提取單元剛度矩陣：
【ABAQUS 二次開發筆記】輸出單元剛度矩陣 - hayden_william - 博客園
以上均為我的一點理解，不一定完全正確，本文僅作為個人學習記錄之用，其他概不負責。
展開 ABAQUS 剛度矩陣 ￥1000
我有個abaqus的問題，你們幫幫我出出主意。 是這樣的：
1，TestElement.py 是編寫的測試代碼，可輸出8結點線性單元的剛度矩陣。
2，abaqus文件Job-testing.inp ，運行可以輸出單元剛度矩陣。
問題是：他們的結點坐標/排序，材料參數都一致， 但得到的剛度矩陣就不一樣。
需要：代碼輸出的剛度矩陣與abaqus得到的剛度矩陣一致或者基本一致。
本期給大家推薦一款由木木自研的小工具：TransferMatrix，主要用于導出 Abaqus 中各種剛度矩陣，方便和自研的程序進行實時對比，基于 Pyside6 搭建的軟件界面，用戶僅需導入 inp 文件，就可以自動調用電腦內的 Abaqus 進行計算分析，不受限于 Abaqus 版本，可以導出：
單元剛度矩陣
單元質量矩陣
單元分布節點荷載列陣
整體剛度矩陣
整體質量矩陣
可選是否導出后打開 Matlab ，并保存為 mat 文件
默認保存的格式為 csv，用戶可以直接在TransferMatrix 中查看
可選是否保存整體剛度/質量矩陣為稀疏矩陣/全矩陣形式
軟件界面見下圖：
若勾選導出后打開matlab，可直接查看：
軟件實現原理
使用注意事項
軟件啟動速度較慢，需要有點耐心，鄙人對于 PySide 的打包技術認識有限
若 inp 文件內沒有密度參數，是不能導出質量矩陣的，建議在動力類型的分析步中進行導出質量矩陣
非協調單元和雜交單元只能輸出質量矩陣
在輸出載荷列陣時，Abaqus 產生的是單元節點的分布載荷，如果沒有這中載荷類型，也是無法導出的
建議再導入 inp 文件之前，確保這個文件能在 Abaqus 上跑通，本軟件只是負責轉換剛度矩陣
不支持在TransferMatrix 中查看整體剛度/質量矩陣，因為一般情況下全局矩陣尺寸很大，可選擇導出至 Matlab 中進行查看
如果出現風險提示，請選擇信任：
軟件同級目錄放置的config.json文件，里面配置的是本地的 abaqus.bat 和 matlab.exe 目錄，建議用 everything 工具搜索出這兩個的安裝位置，復制進去即可。
展開 
abaqus接觸剛度的相關專題、標簽、搜索
abaqus接觸剛度的最新內容
[圖片]
Job-3.inp
模型一直在提示這個錯，哪位大佬可以幫忙解決一下，有償，主要問題就在旋輪和坯料接觸的地方，接觸算法修改過，也不行，一上來就有錯，但是可以算，算一段時間后就因為錯誤而中斷了。
目錄
初始過盈產生的原因
解決初始過盈的方法
在 Abaqus 中指定不同選項
壓力工況下的性能驗證
總結
在本文中，我們將以軸對稱 O 型圈為例，闡述并展示 Abaqus 處理接觸表面初始重疊（即初始過盈）的多種方法。
<p>Cohesive作為ABAQUS中常用的粘結技術，無論在模擬粘結界面（例如新舊混凝土疊合面、復合材料粘結界面）或是全局粘結單元（例如模擬細觀混凝土開裂）具有較廣泛的應用。今天喵星人從官方的用戶手冊中選取了幾個Cohesive基礎而又關鍵的知識點，幫助大家修煉Cohesive內功。</p><p class="ql-align-center"><strong>01</strong>內聚力單元/接觸區別
<p><span style="color: rgb(0, 0, 0);">接觸問題是ABAQUS數值模擬中常遇到的典型分析問題，在結構工程中裝配式結構（包含裝配式混凝土結構、裝配式鋼結構、裝配式組合結構等）均會涉及到大量的接觸面，且接觸面的屬性設置可能存在不同。當大家在操作這些錯綜復雜的接觸面時難免無從下手，今天就看喵星人如何輕松拿捏錯綜復雜的接觸面吧~</span></p><p class="ql-align-center
[圖片]
本期給大家推薦一款由木木自研的小工具：TransferMatrix，主要用于導出 Abaqus 中各種剛度矩陣，方便和自研的程序進行實時對比，基于 Pyside6 搭建的軟件界面，用戶僅需導入 inp 文件，就可以自動調用電腦內的 Abaqus 進行計算分析，不受限于 Abaqus 版本，可以導出：
單元剛度矩陣
單元質量矩陣
單元分布節點荷載列陣
整體剛度矩陣
ABAQUS 剛度矩陣11個月前
我有個abaqus的問題，你們幫幫我出出主意。 是這樣的：
1，TestElement.py 是編寫的測試代碼，可輸出8結點線性單元的剛度矩陣。
2，abaqus文件Job-testing.inp ，運行可以輸出單元剛度矩陣。
問題是：他們的結點坐標/排序，材料參數都一致， 但得到的剛度矩陣就不一樣。
需要：代碼輸出的剛度矩陣與abaqus得到的剛度矩陣一致或者基本一致
