abaqus沙漏剛度的案例
ABAQUS疑難雜癥之沙漏剛度
源文件可在公眾號:‘易木木響叮當’ 內回復“沙漏剛度”自動獲取。
自研有限元程序的減縮積分單元如何添加沙漏剛度（理論解釋+數值實現）
公式排版、代碼排版效果不佳，所以上傳的圖片，見諒
abaqus DEM做的沙漏
*particle generator flow，幫助文檔上有說明，
在本平臺中也看到玉蝴蝶做的沙漏，就想自己也做做，生成粒子沒有問題，卡在旋轉上，真的卡了一段時間，最近看到abaqus官網上有個例子，看了看，找到問題關鍵，終于得以解決，貼出來大家一起欣賞吧
，
ABAQUS模擬沙漏
ABAQUS模擬沙漏
轉載，Abaqus分析中的沙漏效應
沙漏的定義
沙漏（hourglassing）的產生是一種數值問題，是單元自身存在的一種數值問題。一般出現在采用縮減積分單元的情況下：
比如一階四邊形縮減積分單元，該單元有四個節點 “o”，但只有一個積分點“*”。而且該積分點位于單元中心位置，此時如果單元受彎曲或者受剪切作用，則必然會發生變形，如下圖所示。
對于單積分點線性單元，單元雖然受力后產生變形，但并沒有產生應變能--也叫零能量模式。在這種情況下，單元沒有剛度，所以不能抵抗變形，顯然這樣的結論是不合理的，所以必須避免這種情況的出現，需要加以控制，既然沒有剛度，就要施加虛擬的剛度以限制沙漏模式的擴展---人為加的沙漏剛度就是這么來的。
2.沙漏的控制方法
目前常用的沙漏控制算法大致分為兩類：粘性阻尼算法和彈性剛度算法。這兩種算法分別通過引入沙漏變形方向上的阻尼約束力和剛度約束力來控制沙漏變形。
由于引入了沙漏控制力，同時就會產生沙漏能量損失，對于系統的能量平衡產生影響。在某些工程問題中，采用沙漏控制方法并不能完全解決沙漏問題，對于這類問題，可采用多點積分的單元來解決，當然計算成本也會大大增加。
3.ABAQUS中沙漏的設置
在ABAQUS/CAE中，可以方便地在Element Type界面下進行沙漏的設置。
1、Distortion control：只用于explicit分析。
當選擇 YES時，激活防止負體積單元出現或其他可壓縮材料的過度變形，這對超彈材料是默認的。Distortion control參數對線性動力學不可用，并且不能防止單元由于時間不穩定、沙漏不穩定或不切實際的物理變形造成的扭曲。
展開 Abaqus分析中的沙漏效應
沙漏的定義
沙漏（hourglassing）的產生是一種數值問題，是單元自身存在的一種數值問題。一般出現在采用縮減積分單元的情況下：
比如一階四邊形縮減積分單元，該單元有四個節點 “o”，但只有一個積分點“*”。而且該積分點位于單元中心位置，此時如果單元受彎曲或者受剪切作用，則必然會發生變形，如下圖所示。
對于單積分點線性單元，單元雖然受力后產生變形，但并沒有產生應變能--也叫零能量模式。在這種情況下，單元沒有剛度，所以不能抵抗變形，顯然這樣的結論是不合理的，所以必須避免這種情況的出現，需要加以控制，既然沒有剛度，就要施加虛擬的剛度以限制沙漏模式的擴展---人為加的沙漏剛度就是這么來的。
2. 沙漏的控制方法
目前常用的沙漏控制算法大致分為兩類：粘性阻尼算法和彈性剛度算法。這兩種算法分別通過引入沙漏變形方向上的阻尼約束力和剛度約束力來控制沙漏變形。
由于引入了沙漏控制力，同時就會產生沙漏能量損失，對于系統的能量平衡產生影響。在某些工程問題中，采用沙漏控制方法并不能完全解決沙漏問題，對于這類問題，可采用多點積分的單元來解決，當然計算成本也會大大增加。
3.ABAQUS中沙漏的設置
在ABAQUS/CAE中，可以方便地在Element Type界面下進行沙漏的設置。
1、Distortion control：只用于explicit分析。
當選擇 YES時，激活防止負體積單元出現或其他可壓縮材料的過度變形，這對超彈材料是默認的。Distortion control參數對線性動力學不可用，并且不能防止單元由于時間不穩定、沙漏不穩定或不切實際的物理變形造成的扭曲。
展開 abaqus離散元沙漏模型
*Particle Generator, name=dem1, type=PD3D,
Maximum Number of Particles=2000
*Particle Generator Inlet, surface=inlet1
*Particle Generator Mixture
dem1,
*Discrete Section, elset=dem1, density=0.25, alpha=7.0
PDF1,
*Discrete Elasticity
2.0E+10, 0.25
*Probability Density Function, name=PDF1, TYPE=DISCRETE
0.5, 1.0
______
*Particle Generator Flow, generator=Particle1-1.dem1
FlowSpeed, MassFlow
**
*Contact
*Contact Controls Assignment, rotational terms=STRUCTURAL
*Contact Inclusions
dem1, loudou-1.in
dem1, dem1
*Contact Property Assignment
dem1, loudou-1.in, P1F
dem1, dem1, P11
*dload
particle1-1.dem1, GRAV, 9800., 0., -1., 0.
展開 ABAQUS動態分析中的能量平衡、沙漏及結果評估
ABAQUS動態分析中的能量平衡、沙漏及結果評估
Abaqus使用umat子程序的沙漏問題
在使用umat子程序時，采用縮減積分單元后，沙漏控制剛度是通過材料屬性中的彈性性質定義的，這些剛度基于材料初始剪切模量的值。但是在使用umat時，Abaqus對程序輸入文件進行預處理時得不到剪切模量的數值，所以這時候必須通過hourglass stiffness定義具有沙漏模式的單元的沙漏控制剛度。
在Abaqus/CAE 的單元模塊即可定義沙漏控制模式，如圖所示。
ABAQUS中的單元選擇-理解剪切自鎖和沙漏
未能對剛度矩陣進行精確積分的誤差可以通過加密網格來彌補，也即加密縮減積分網格的數量可以有效改善計算結果。
但是同時注意到，單元中心點（積分點）x和y方向上的長度也沒有變化（圖4），基于該點所計算的整個單元的應變能將等于0，單元表現為零剛度，單元過軟，稱為沙漏hourglass現象。在網格較粗時，零剛度現象將通過網格進一步擴散，從而產生無意義的結果。為了限制沙漏現象的擴展，ABAQUS引進了“防沙漏剛度”Hourglass stiffness，一般情況下采用默認值即可，如果確有需要可在圖1中的Hourglass control選項中設置。
圖4
四、小結
如果模型中有比較明顯的彎曲現象，為避免出現剪切自鎖現象，優先選擇二階單元，或者采用縮減積分方案（網格需要更細，通常厚度方向4層以上）。
來源： ABAQUS在巖土工程中的應用
展開 Abaqus有限元仿真分析中的沙漏控制方法與設置
abaqus中沙漏的產生是一種數值問題，單元自身存在的一種數值問題，舉個例子，對于單積分點線性單元，單元受力變形沒有產生應變能，也叫0能量模式，在這種情況下，單元沒有剛度，所以不能抵抗變形，不合理，所以必須避免這種情況的出現，需要加以控制，既然沒有剛度，就要施加虛擬的剛度以限制沙漏模式的擴展，人為加的沙漏剛度就是這么來的。下面，小編就給大家分享一下"Abaqus有限元仿真分析中的沙漏控制方法與設置"。
沙漏的定義
沙漏（hourglassing）的產生是一種數值問題，是單元自身存在的一種數值問題。一般出現在采用縮減積分單元的情況下：比如一階四邊形縮減積分單元，該單元有四個節點“o”，但只有一個積分點"*"。而且該積分點位于單元中心位置，此時如果單元受彎曲或者受剪切作用，則必然會發生變形，如下圖所示。
單元原始狀態、單元受剪切作用變形、單元受彎曲作用變形
對于單積分點線性單元，單元雖然受力后產生變形，但并沒有產生應變能--也叫零能量模式。在這種情況下，單元沒有剛度，所以不能抵抗變形，顯然這樣的結論是不合理的，所以必須避免這種情況的出現，需要加以控制，既然沒有剛度，就要施加虛擬的剛度以限制沙漏模式的擴展，人為加的沙漏剛度就是這么來的。
沙漏的控制方法
目前常用的沙漏控制算法大致分為兩類：粘性阻尼算法和彈性剛度算法。這兩種算法分別通過引入沙漏變形方向上的阻尼約束力和剛度約束力來控制沙漏變形。
由于引入了沙漏控制力，同時就會產生沙漏能量損失，對于系統的能量平衡產生影響。在某些工程問題中，采用沙漏控制方法并不能完全解決沙漏問題，對于這類問題，可采用多點積分的單元來解決，當然計算成本也會大大增加。
ABAQUS中沙漏的設置
在ABAQUS/CAE中，可以方便地在ElementType界面下進行沙漏的設置。
展開 
【ABAQUS模態動力學】Composite&abaqus 預應力模態分析&輸出單元剛度矩陣
劃分網格
定義邊界條件
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得到的剛度矩陣一致或者基本一致。
有限元理論基礎及Abaqus內部實現方式研究系列3：S4殼單元剪切自鎖和沙漏控制
圖1：剪切自鎖
圖2：沙漏
===S4殼單元剪切自鎖和沙漏控制研究總結===
完全積分單元才有剪切自鎖，雖然Abaqus的S4單元是完全積分，但內部已經做了修正完全消除了剪切自鎖，所以不需要用戶做任何設置。
減縮積分單元才有沙漏現象，Abaqus的S4R默認增加一個人工的沙漏剛度來控制沙漏現象，如果發現結果還是不理想，那么需要采用其它建模方法才能控制沙漏了。
Abaqus針對剪切自鎖和沙漏控制做的修正如下：
單元類型
元素
修正情況
剪切自鎖
S4
薄膜剛度
假設應變方法修正
面外彎曲剛度
無
面外橫向剪切剛度
假設應變方法修正
S4R
所有項
無
沙漏控制
S4
薄膜剛度
無
面外彎曲剛度
無
面外橫向剪切剛度
沙漏控制，和Belytchko公式不一致，暫時沒研究
S4R
薄膜剛度
和Belytchko公式一致，因子取0.005G
面外彎曲剛度
和Belytchko公式一致，因子取0.00375G*4
面外橫向剪切剛度
沙漏控制，和Belytchko公式不一致，暫時沒研究
詳細研究方法，見附件：
有限元理論基礎及Abaqus內部實現方式研究系列3：S4殼單元剪切自鎖和沙漏控制(SnowWave02 20171018).pdf
展開 TransferMatrix：導出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 工具搜索出這兩個的安裝位置，復制進去即可。
展開 