帶孔拉伸的案例
帶孔復合材料拉伸模型 ￥2
B站帶孔復合材料拉伸模型
帶孔復合材料拉伸（0厚度cohesive模型）文件 ￥3
B站帶孔復合材料拉伸（0厚度cohesive模型）文件
帶孔復合材料層合板拉伸失效VUMAT子程序 ￥45
B站子程序VUMAT，帶孔復合材料拉伸
Abaqus纖維復合材料帶孔板層合板拉伸仿真模型 ￥89
<div contenteditable="false" width="100%">
Abaqus纖維復合材料帶孔板層合板拉伸仿真模型!
</div><div contenteditable="false" width="100%">
模擬過程采用3d puck子程序
</div><div contenteditable="false" width="100%">
內附VUMAT子程序，inp文件及ODB文件
</div><p><br></p><p><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202504/attachment/83a0530cf0f442ea88ccb4c9d4c4d4b1.jpg" style="display: inline-block;" data-regular="true">
<img src="https://img.jishulink.com/202504/attachment/83a0530cf0f442ea88ccb4c9d4c4d4b1.jpg" data-mobile-src="https://img.jishulink.com/202504/attachment/83a0530cf0f442ea88ccb4c9d4c4d4b1.jpg?
展開 
各向同性硬化von Mises率無關彈塑性本構理論以及umat源代碼 ￥99
二者計算的結果對比如下：
二者等效塑性應變的演化對比圖如下：
二者塑性耗散的演化對比圖為：
4.2 帶孔板拉伸實驗
利用相同的材料性質計算一帶孔板的拉伸實驗。Abaqus計算的結果如下：
使用umat計算的結果如下：
可看到二者的結果是完全相同的。
5 參考書籍
Neto, E. A. de Souza, D. R. J. Owen, and D. Peric. , 'Computational Methods for Plasticity: Theory and Applications'
溫度依賴線彈性以及考慮熱脹冷縮的umat源代碼 ￥20
其中材料的性質與溫度相關，根據插值得到，相關代碼如下：</p><p><br></p><pre class="ql-syntax" spellcheck="false">// compute elastic properties by interpolating
double interpolate_elastic_properties(double x0, double x1, double y0, double y1, double x)
{
double fac0 = 1.0 - (x - x0) / (x1 - x0);
double fac1 = (x - x0) / (x1 - x0);
if (x < x0)
{
return y0;
}
else if (x > x1)
{
return y1;
}
else
{
return fac0 * y0 + fac1 * y1;
}
}
</pre><p><br></p><p>3 測試</p><p>3.1 帶孔板拉伸</p><p>Abaqus的計算結果為：</p><div contenteditable="false" width="100%"><figure class="figure-image" data-img
展開 沖壓件加工廠在沖壓工序設計上有什么原則?
沖壓件加工廠在沖壓件加工中，各道工序是有先后次序的，沖壓的先后次序是根據工件的形狀、尺寸精度、以及材料變形的規律來看的，在加工時候是有一定的原則的，主要有下面幾點；
1、沖壓件加工廠對于帶孔或有缺口的沖壓件，選用單工序模時，通常先落料再沖孔或缺口，選用級進模時，則落料安排為最后工序；
2、如果工件上存在位置靠近、大小不一的兩個孔，則應先沖大孔后沖小孔，以免大孔沖裁時材料變形引起小孔的形變；
3、沖壓件加工廠對于帶孔的彎曲件，在一般情況下，可以先沖孔后彎曲，當孔位位于彎曲變形區或接近變形區，以及孔與基準面有較高要求時，則應先彎曲后沖孔；
4、沖壓件加工廠對于帶孔的拉伸沖壓件，一般先拉伸后沖孔，當孔的位置在工件底部、且孔的尺寸精度要求不高時，可以先沖孔在拉深，這樣可以減少拉深次數；
5、多角彎曲件應從材料變形影響和彎曲時材料的偏移趨勢安排彎曲的順序，一般應先彎外角后彎內角；
6、對于復雜的旋轉體拉深件，一般先拉伸大尺寸的外形，后拉伸小尺寸的內形，對于復雜的非旋轉拉深尺寸的，應先拉伸小尺寸的內形，后拉伸大尺寸的外部形狀；
7、整形工序、校平工序、切邊工序等，都安排到基本成形之后；
展開 非等溫各向同性線彈性umat開發 ￥20
* y: 待求點材料參數值
function interpolate_para(npt, xp, yp, x) result(y)
integer,intent(in) :: npt
real(8),intent(in) :: xp(npt)
real(8),intent(in) :: yp(npt)
real(8),intent(in) :: x
real(8) :: y
integer :: i
if ( x < xp(1) ) then
y = yp(1)
elseif ( x > xp(npt) ) then
y = yp(npt)
else
do i = 1, npt-1
if ( x >= xp(i) .and. x <= xp(i+1) ) then
y = ( yp(i+1) - yp(i) ) / ( xp(i+1) - xp(i) ) * ( x - xp(i) ) + yp(i)
endif
enddo
endif
return
end function
5 算例
5.1 單個單元單軸拉伸
計算的應力結果對比如下（左圖為Abaqus自帶本構計算結果，右圖為子程序計算結果）：
反力曲線對比如下：
彈性應變的對比曲線為：
熱應變的對比曲線為：
5.2 帶孔平板拉伸
計算的應力結果對比如下（左圖為Abaqus自帶本構計算結果，右圖為子程序計算結果）：
6 源代碼
展開 【5月CAE培訓】Code_Aster非線性分析專題培訓通知
5月12日
上午：
（1）Code_Aster軟件介紹；
（2）Code_Aster使用指導；
（3）練習1：儲水圓柱容器壁面的靜力分析；
下午：
（4）實例練習1（續）：儲水圓柱容器壁面的靜力分析；
（5）Code_Aster命令詳解；
（6）練習2：帶孔平板拉伸的線彈性分析；
5月13日
上午：
（1）Code_Aster非線性基礎、非線性分析法中的連續法指導；
（2）練習3：帶孔平板拉伸的彈塑性分析；
（3）練習4：框架倒塌分析；
下午：
（4）練習4（續）：框架倒塌分析；
（5）Code_Aster非線性專題指導，接觸與摩擦問題；
（6）練習5：Hertz接觸問題分析；
5月14日
展開 前處理及后處理對有限元結果的影響分析
本文僅以簡單帶孔平板的拉伸分析，對比分析了網格尺寸和后處理的應力位置對結果的影響。希望對新手有所啟發和幫助。
1、首先通過力學基本理論計算了基準應力，作為有限元分析結果的標準值，計算過程如下：
2、網格尺寸對仿真結果的影響分析：
3、應力位置對結果理解的影響分析：
4、結論：
分析結構應力時，從三個層次考慮結構的受力及失效風險：
①應力分布的合理性
②最大應力的位置
③應力值的準確性
5、應用推廣：
①對于應力集中區域，應該分析單元尺寸對結果的影響；
②對于鈑金幾何邊界的應力值，建議使用單元角點應力查看；
③對于實體幾何邊界的應力值，建議使用表面單元應力查看。
文章原創，轉發請注明來源作者，@元來是你。
評論區點贊并留下郵箱，可獲得原模型和分析文章。
展開 HyperMesh與ANSYS聯合仿真（一）
下面筆者用一個簡單的
帶孔薄板拉伸（
平面應力問題）的例子來講解一下HyperMesh與ANSYS聯合仿真的關鍵步驟及注意事項。
本例仍然使用公眾號文章《ANSYS與材料力學之軸向拉伸和壓縮（六）》中使用
的模型、載荷及邊界條件。
Step1：設置求解器選項。
打開HyperWorks2020，在File中將Solver Interface設置為Ansys。
Step2：建立幾何模型。
在HyperWorks的Geometry模塊中建立帶孔薄板的平面模型如下圖所示。長為20mm，寬為10mm，孔徑為2mm。厚度設置為0.1mm（在平面單元屬性中定義）。
Step3：創建Sensor來存儲單元類型。
在Model模型樹下的空白處右擊選擇Create→Sensor，并將其命名為“PLANE”。單擊“PLANE”，將Element Type改為“PLANE182”，其他選項保持默認，從而創建出我們計算時需要的單元。
Step4：創建Property來存儲單元屬性。
在Model模型樹下的空白處右擊選擇Create→Property，并將其命名為“PLANE”。單擊“PLANE”，將Card Image改為“PLANE182p”；將Real Constants中的THK（厚度）設置為0.1。其他選項保持默認。
Step5：創建 Material來存儲材料屬性。
在Model模型樹下的空白處右擊選擇Create→Material，并將其命名為“PLANE”。單擊“PLANE”，將楊氏模量改為“200000”；將主泊松比設置為0.3。其他選項保持默認。
Step6：網格劃分。
展開 
HyperMesh與ANSYS聯合仿真（一）
下面筆者用一個簡單的
帶孔薄板拉伸（平面應力問題）的例子來講解一下HyperMesh與ANSYS聯合仿真的關鍵步驟及注意事項。
本例仍然使用公眾號文章《ANSYS與材料力學之軸向拉伸和壓縮（六）》中使用
的模型、載荷及邊界條件。
Step1：設置求解器選項。
打開HyperWorks2020，在File中將Solver Interface設置為Ansys。
Step2：建立幾何模型。
在HyperWorks的Geometry模塊中建立帶孔薄板的平面模型如下圖所示。長為20mm，寬為10mm，孔徑為2mm。厚度設置為0.1mm（在平面單元屬性中定義）。
Step3：創建Sensor來存儲單元類型。
在Model模型樹下的空白處右擊選擇Create→Sensor，并將其命名為“PLANE”。單擊“PLANE”，將Element Type改為“PLANE182”，其他選項保持默認，從而創建出我們計算時需要的單元。
Step4：創建Property來存儲單元屬性。
在Model模型樹下的空白處右擊選擇Create→Property，并將其命名為“PLANE”。單擊“PLANE”，將Card Image改為“PLANE182p”；將Real Constants中的THK（厚度）設置為0.1。其他選項保持默認。
Step5：創建 Material來存儲材料屬性。
在Model模型樹下的空白處右擊選擇Create→Material，并將其命名為“PLANE”。
展開 基于VUMAT復合材料夾層結構沖擊仿真
就比如帶孔板的拉伸和圓柱擾流，云圖真是傻傻難分。并且還真有學者找到兩者的本構共通性，用圓柱繞流來研究帶孔板。
一下子扯遠了，繼續低速沖擊。
2 模型
考慮如下模型，邊界條件為底部固支，上面板四個角點固支。面板失效基于Hashin準則判斷，芯層失效基于MISESS準則判斷。
看似簡單的問題往往暗含殺“雞”。
1） 沖頭設置為剛體，其密度的取值，不能直接賦予鋼的屬性。因為實際的沖頭結構為柱狀。建模中，處于簡化考慮，取頭部半球進行建模，為此需要根據實際沖頭質量，換算出仿真用沖頭的密度。
2） 面板和夾層之間可以綁定，如果夾層是蜂窩這類非均勻結構，用接觸屬性會比較合適，但是接觸的定義要考慮好，否則很容易穿透，或者大滑移。
3） 同樣的，沖頭和面板的接觸也要注意，網格的疏密和接觸屬性都可能造成穿透。
4） 為了防止網格過度扭曲，要對網格扭曲進行控制，也可以縮放其質量，或者對過度扭曲的單元，直接賦予高模量。
5） 載荷為速度載荷，如果已知沖擊能量，就根據沖頭質量進行速度換算，這是高中知識了。
3 VUMAT
1) 我們此次使用VUMAT最重要的目的是，實現失效區域的識別。
2) 本次VUMAT關鍵輸出，是應力的更新和損傷變量的更新。
3) 本次VUMAT的結構如下：
材料參數定義；
HASHIN準則計算；
用HASHIN準則進行失效判斷；
對失效區域的剛度進行折減；
剛度矩陣定義；
應變定義；
根據剛度矩陣和應變進行應力更新；
4 結果
應力云圖效果如下：
載荷-時間曲線如下，橫坐標單位是ms。
面板和芯層的失效各自如下：
需要指出的是，這個模型相對簡單，芯層結構比較均勻。
展開 從C3D8的uel源代碼入門Abaqus的uel編寫, 更新B-Bar修正 ￥99
</p><p>位移結果為：</p><p><img src="https://img.jishulink.com/msimage/202403/59d3be0b786fb208327fb78b42740638.png"></p><p>應變結果為：</p><p><img src="https://img.jishulink.com/msimage/202403/12da86b65a5e0bc4061254cb86ca312b.png"></p><p>應力結果為：</p><p><img src="https://img.jishulink.com/msimage/202403/388f17d6c05d463fddabdb333b62cec1.png"></p><h2>4.2 帶孔板單軸拉伸測試</h2><p>位移結果對比為（左列為Abaqus計算結果，右列為uel計算結果）：</p><p><img src="https://img.jishulink.com/msimage/202403/e8adda89596f17a0ac201d0eb821b53e.png"></p><p>應變的對比結果為（左列為Abaqus計算結果，右列為uel計算結果）：</p><p><img src="https://img.jishulink.com/msimage/202403/7521fa74e6916798eec6814b2c56095b.png"><img src="https://img.jishulink.com/msimage/202403/c7e43ad9a83375e613aa728be1fb2df7.png"></p><p>應力的對比結果為（左列為Abaqus計算結果，右列為uel計算結果）：</p><h1><img src="https://img.jishulink.com/msimage
展開 隨動硬化von Mises率無關彈塑性本構理論以及umat源代碼 ￥99
</p><p>二者應力應變滯回曲線對比如圖：</p><p><img src="https://img.jishulink.com/msimage/202402/2787199530f7487cfbf12d5376302c9c.png"></p><p>二者等效塑性應變演化曲線對比如圖：</p><p><img src="https://img.jishulink.com/msimage/202402/118f2dbd811cb1264e23d2812f932643.png"></p><p>二者塑性耗散演化曲線對比如圖：</p><p><img src="https://img.jishulink.com/msimage/202402/7f4ac8a5b98de7e140528657fc076e55.png"></p><p><br></p><p>4.2 帶孔板拉伸測試</p><p>設置的材料參數為：</p><p><img src="https://img.jishulink.com/msimage/202402/fdd4e3f61adc75058e0ad3e2a6a9967a.png"></p><p><br></p><p>使用Abaqus計算的von Mises應力以及等效塑性應變的結果如下：</p><p><img src="https://img.jishulink.com/msimage/202402/7ca4fc0978422341fabf755cb84c01ee.png"></p><p><img src="https://img.jishulink.com/msimage/202402/50dddd6a4f772a03cf5a03895601d409.png"></p><p>使用umat計算的von Mises應力以及等效塑性應變的結果如下：</p><p><img src="
展開 