單元沙漏 abaqus的案例
有限元理論基礎及Abaqus內部實現方式研究系列3：S4殼單元剪切自鎖和沙漏控制
這兩者具有相似性，所以我們在本文中一起研究Abaqus中線性殼單元S4針對這兩種情況下的內部實現方式。剪切自鎖和沙漏現象影響的是剛度矩陣和應力，我們研究方式是在自編程序iSolver中根據成熟的消除剪切自鎖和沙漏控制的理論實現剛度矩陣的修正，通過比較同一模型的Abaqus的剛度矩陣結果，結合幫助文檔猜測Abaqus軟件單元消除剪切自鎖和控制沙漏的內部實現方法。
圖1：剪切自鎖
圖2：沙漏
===S4殼單元剪切自鎖和沙漏控制研究總結===
完全積分單元才有剪切自鎖，雖然Abaqus的S4單元是完全積分，但內部已經做了修正完全消除了剪切自鎖，所以不需要用戶做任何設置。
減縮積分單元才有沙漏現象，Abaqus的S4R默認增加一個人工的沙漏剛度來控制沙漏現象，如果發現結果還是不理想，那么需要采用其它建模方法才能控制沙漏了。
展開 ABAQUS中的單元選擇-理解剪切自鎖和沙漏
但是同時注意到，單元中心點（積分點）x和y方向上的長度也沒有變化（圖4），基于該點所計算的整個單元的應變能將等于0，單元表現為零剛度，單元過軟，稱為沙漏hourglass現象。在網格較粗時，零剛度現象將通過網格進一步擴散，從而產生無意義的結果。為了限制沙漏現象的擴展，ABAQUS引進了“防沙漏剛度”Hourglass stiffness，一般情況下采用默認值即可，如果確有需要可在圖1中的Hourglass control選項中設置。
圖4
四、小結
如果模型中有比較明顯的彎曲現象，為避免出現剪切自鎖現象，優先選擇二階單元，或者采用縮減積分方案（網格需要更細，通常厚度方向4層以上）。
來源： ABAQUS在巖土工程中的應用
展開 有限元筆記#2：什么是沙漏現象？為什么減縮積分線性單元會存在沙漏問題？
如果使用了減縮積分的線性單元，即使不是在純彎曲加載模式下，其得到的應力應變值相比理論預示值應該要小（我推測的^_^，沒空詳細證實），所以用這樣的數據構造的切線剛度矩陣相比其他單元構造的切線剛度矩陣要小，這也許就是通常所說的出現沙漏問題的單元“太軟”的緣故。
結語：本文算不得什么，只是從公式上加深了商業軟件使用者對沙漏這一現象的了解，稍微知其所以然罷了。如果要進一步探究如何防止沙漏問題，要構造怎樣的位移模式，需要更多功夫，可見如下博文：
FEMer，公眾號：易木木響叮當
減縮積分單元、沙漏控制與自定義單元：與Abaqus C3D8R單元的精度對比之旅
參考資料：
《有限元分析及應用》曾攀，清華大學出版社，2004.
《基于ABAQUS的有限元分析和應用》莊茁等，清華大學出版社2008.
展開 自研有限元程序的減縮積分單元如何添加沙漏剛度（理論解釋+數值實現）
公式排版、代碼排版效果不佳，所以上傳的圖片，見諒
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疑難雜癥之沙漏剛度
沙漏模式是怎樣產生的？
沙漏模式指的是非物理的零能變形模式，簡而言之有變形沒有應力或應變。如果采用單點積分（積分點在等參元中心），在某些情況下節點位移不為零，即單元有形變，但插值得到的應變卻為零。主要出現在CPS4R、CAX4R、C3D8R等線性減縮積分單元中，因為線性減縮積分在每個方向上只有一個積分點，應變插值后有可能為零。打個比方說，一個正方體單元變形為等腰梯形，節點位移相等卻方向相反，各點的形函數為零，所以插值結果為零，這樣內能計算結果也為零，即計算認識單元沒有變形，與事實相反。
如何判斷出現沙漏模式？
方法一：
在Visualization功能模塊中選擇菜單Result、History Output，繪制偽應變能ALLAE(artificial strain energy)和內能ALLIE(internal energy)的曲線，當偽應變能ALLAE約占內能ALLIE的1%時，表明沙漏模式對計算結果影響不大；當偽應變能超過內能10%時，分析就是無效的，必須采取措施加以解決。
方法二：
查看單元變形過程：如果有單元變形明顯異常，或有單元變成交替出現的梯形形狀，一般是出現沙漏模式。
應該怎樣避免？
細化網格。
展開 ABAQUS動態分析中的能量平衡、沙漏及結果評估
ABAQUS動態分析中的能量平衡、沙漏及結果評估
Abaqus使用umat子程序的沙漏問題
在使用umat子程序時，采用縮減積分單元后，沙漏控制剛度是通過材料屬性中的彈性性質定義的，這些剛度基于材料初始剪切模量的值。但是在使用umat時，Abaqus對程序輸入文件進行預處理時得不到剪切模量的數值，所以這時候必須通過hourglass stiffness定義具有沙漏模式的單元的沙漏控制剛度。
在Abaqus/CAE 的單元模塊即可定義沙漏控制模式，如圖所示。
Abaqus有限元仿真分析中的沙漏控制方法與設置
abaqus中沙漏的產生是一種數值問題，單元自身存在的一種數值問題，舉個例子，對于單積分點線性單元，單元受力變形沒有產生應變能，也叫0能量模式，在這種情況下，單元沒有剛度，所以不能抵抗變形，不合理，所以必須避免這種情況的出現，需要加以控制，既然沒有剛度，就要施加虛擬的剛度以限制沙漏模式的擴展，人為加的沙漏剛度就是這么來的。下面，小編就給大家分享一下"Abaqus有限元仿真分析中的沙漏控制方法與設置"。
沙漏的定義
沙漏（hourglassing）的產生是一種數值問題，是單元自身存在的一種數值問題。一般出現在采用縮減積分單元的情況下：比如一階四邊形縮減積分單元，該單元有四個節點“o”，但只有一個積分點"*"。而且該積分點位于單元中心位置，此時如果單元受彎曲或者受剪切作用，則必然會發生變形，如下圖所示。
單元原始狀態、單元受剪切作用變形、單元受彎曲作用變形
對于單積分點線性單元，單元雖然受力后產生變形，但并沒有產生應變能--也叫零能量模式。在這種情況下，單元沒有剛度，所以不能抵抗變形，顯然這樣的結論是不合理的，所以必須避免這種情況的出現，需要加以控制，既然沒有剛度，就要施加虛擬的剛度以限制沙漏模式的擴展，人為加的沙漏剛度就是這么來的。
沙漏的控制方法
目前常用的沙漏控制算法大致分為兩類：粘性阻尼算法和彈性剛度算法。這兩種算法分別通過引入沙漏變形方向上的阻尼約束力和剛度約束力來控制沙漏變形。
由于引入了沙漏控制力，同時就會產生沙漏能量損失，對于系統的能量平衡產生影響。在某些工程問題中，采用沙漏控制方法并不能完全解決沙漏問題，對于這類問題，可采用多點積分的單元來解決，當然計算成本也會大大增加。
ABAQUS中沙漏的設置
在ABAQUS/CAE中，可以方便地在ElementType界面下進行沙漏的設置。
展開 abaqus實體-梁單元，實體-實體單元，梁-梁單元鉸接設置
使用多點約束MPC，實現實體-梁單元，實體-實體單元，梁-梁單元鉸接如何設置，實體單元梁彎矩曲線怎么提取？可下載附件，也可觀看視頻。
https://www.yqgqt.org.cn/college/video/c15810?nagivator=course
abaqus實體-梁單元，實體-實體單元，梁-梁單元鉸接設置.rar
BCC點陣結構梁單元Abaqus壓縮仿真模擬-顯示動力學質量縮放 ￥19.89
本文通過abaqus顯示動力學的方法對BCC結構進行壓縮仿真模擬，同時為減小計算量，采用梁單元模擬點陣結構，壓頭設置為剛性面，添加質量縮放，加快運算速度，為點陣結構壓縮模擬提供一種便捷方法。
1. 建立BCC點陣模型,以單胞尺寸5X5X5為例。
a.首先建立立方體實體，然后對實體進行處理，得到點陣單胞點陣結構。
b.建立單胞BCC梁單元點陣模型,然后進行刪除面的操作，得到單胞BCC點陣結構，接下來進行陣列操作，得到最大外形尺寸為25x25x25的點陣壓縮模擬試件。
C.建立剛性壓板，設置參考點，模擬萬能試驗機壓頭，剛性單元不參與計算，不影響計算結果，加快運算速度。
2. 裝配，按壓縮試驗進行裝配，從上到下依次為壓板-點陣-壓板。
3.設置材料屬性，本文為鈦合金TC4，密度4.43e-9Tone/mm3,彈性模量為118000MPa，泊松比0.3，應力應變值見下表所示。
設置截面屬性Beam，定義截面半徑0.5mm
指派截面，定義方向[1,2,3]，完成材料屬性設置。
4.設置分析步Dynamic,Explicit，時間設置為5s，以每秒1mm的速度進行壓縮模擬，開啟質量縮放為1e-5，歷程輸出勾選位移和力，以便輸出力-位移曲線，然后計算相應的應力-應變曲線。
5.設置相互作用-切向行為和法向行為，摩擦系數為0.3，設置通用接觸。
以下部分為付費部分
展開 