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