mises
關注創建者：tony1 創建時間：2020-12-23
mises的視頻教程
使用Python批量提取多個odb的集合里面的Mises應力并寫入excel里面
本視頻詳細介紹了如何使用python代碼批量提取多個odb中的集合里面的Mises應力隨時間的變化關系，一方面可以在集合中實現提取，另一方面可以在單元編號里面實現提取，對批量處理操作有極大的幫助，可以大量節約人力物力財力。
￥99 8分鐘 40播放
查看
【ABAQUS】二次開發-----從入門到如土
起重臂多工況仿真--同時變角度和受力（for循環的嵌套使用） 混凝土泵車臂架多姿態仿真（多級旋轉，旋轉基點跟隨實例發生變化） 巖土材料非均質建模（單元list的有關操作） 第五部分：Abqus API 解讀 如何進入abaqus 的線上幫助手冊 運行abaqus 腳本的三種方式 abaqus 常用的 cmd 命令及其語法 第六部分：Abqus 后處理二次開發 獲得最大 Mises
￥200 14小時25分鐘 4584播放
查看
梁單元結構建模，optistruct求解后，hyperview查看應力沒有von mises等應力？
做一個對比： 第一個方案：1、梁單元屬性卡片選擇pbeam；2、梁單元截面類型為標準截面庫standard section library：HYPER BEAM下的截面類型。 第二個方案：1、梁單元屬性卡片選擇pbeaml；2、梁單元截面類型為標準截面庫standard section library：OPTISTRUCT下的截面類型 補加單根梁和六面體框架對比 附件是21年版本保存的
￥180 55分鐘 123播放
查看
mises的實例教程
= results.S_Mises(2858:5215);
mises_plate = results.S_Mises(7574:18166);
mises_cover2 = results.S_Mises(5216:7573);
單獨instance顯示
如果想只顯示某一個instance的結果云圖，調用一次meshplot函數即可
figure(3);
meshplot(parts(1), mises_bolt);
clim( [min(mises_bolt), max(mises_bolt)]);
t1=clim;
t1=linspace(t1(1),t1(2),13);
colorbar('ytick',t1,'Location','westoutside');
colormap("jet");
title('mises');
整個裝配體顯示
多次調用meshplot函數進行疊加繪圖即可：
figure(3);
meshplot(parts(1), mises_bolt);
meshplot(parts(2), mises_cover1);
meshplot(parts(3), mises_plate);
meshplot(parts(4), mises_cover2);
% colorbar();
clim( [min(results.S_Mises), max(results.S_Mises)]);
t1=clim;
t1=linspace(t1(1),t1(2),13);
colorbar('ytick',t1,'Location','westoutside');
colormap("jet");
title('mises');
展開 問題：
Ansys workbench進行諧響應仿真計算的后處理結果中，提供了單一頻率下的Von Mises應力查看功能和應力頻響曲線功能，但是應力頻響曲線的應力列表中沒有Von Mises應力查看項。因為Von Mises應力太常用，所以這就給我們在整個掃頻范圍內，定位Von Mises應力的最大頻率和應力值帶來一定的困難。如下所示。
需求：
希望后處理結果中可以在應力響應曲線中，有一項Von Mises應力選項。實現每個掃頻點的最大Von Mises應力和掃頻頻率的曲線圖顯示，從而一眼就看出產品在整個掃頻范圍內，哪個頻率下結構的等效應力最大。而后再通過應力云圖查看這個頻率下的Von Mises應力。
解決方法：
利用APDL命令實現。簡要流程為：首先，讀取每一個掃頻點的最大Von Mises應力值。記下應力值、頻率值和最大節點號。再統計記錄的所有掃頻點的Von Mises應力值，提取整個掃頻過程中最大應力值及其頻率。并將結果寫出到txt文件。進一步提取這個最大Von Mises應力點對應的整個掃頻范圍內的Von Mises應力曲線。
這個樣就可以在txt文檔中直接看到所有掃頻點下，結構的等效應力幅值；以及全頻段中最大Von Mises應力所在節點的等效應力掃頻曲線圖。
效果展示如下：
在結果文件夾中，會生成一個txt結果文件和一張Von Mises應力曲線圖。如此我們可以直觀注意到，在當前掃頻范圍內，結構在78.95Hz時應力最大約為17.552Mpa。
結果后處理問題示例：
Ansys workbench進可以查看某個頻率下的 Von Mises應力幅值
Ansys workbench進掃頻應力響應曲線中，應力選項卻沒有Von Mises應力選型，只能按三個方向來分別查看。
展開 熟悉abaqus的人知道，對于mises應力的問題，abaqus的解決方法是先得到精確的積分點的Mises應力的解，再通過外推的方法得到節點的解。所以當你使用xydate輸出節點MISES應力的時候，會得到多個值，但是如果輸出高斯積分點的MISES時，會得到一個值。
如果你非常喜歡用xydate，你會發現有一個功能無法實現，就是獲得模型整體的最大的Mises應力曲線。因為XYdate是基于固定節點或者單元的，而模型的最大應力單元隨著工況的變化在不斷的變換。
上面說了一些廢話，其實想表達一個意思，再abaqus中的結果文件中是只有高斯積分點的應力的。（如我理解錯誤請及時指出），所以通過python的二次開發，實現abaqus所有增量步下基于積分點的Mises的表格輸出與最大值查找成為可能。
目標1：在結構樹的XYDATE中生成一個新的date，這個date的x選項是所有的增量步（或者時間），Y選項是每個增量步（或者時間）時模型基于積分點的最大Mises應力值。
目標2：找出整個模型MISES應力最大的單元的編號，及對應的step和frame
例：一個典型的螺栓工況，共有4個step，依次是預緊、緊固、固定螺栓長度和加載
1. 打開odb文件（必須操作）
2. 使用開發的插件：MaxMisesCurve，選擇默認選項，點擊OK
3. 3.稍等幾秒。在對話框里會顯示目標2的內容
4. .接下來到結構數中，找到生成的xydate
這里說明下，X不是steptime。因為我想實現總時間，但是沒有找到相應的熟悉，所以我用了一個增量的方法。如果step1實際由3個frame，step2有2個frame。我這里的x則是從0到4.一共5個。
展開 簡單理解Mises應力分量
簡述Mises：
Mises是一種屈服準則，準則的值我們通常叫等效應力，習慣稱Mises等效應力。它遵循材料力學第四強度理論(形狀改變比能理論)。第四強度理論認為形狀改變比能是引起材料流動破壞的主要原因，鋼材等塑性材料遵循第四強度理論，結果更符合實際。
（一般材料在外力作用下產生塑性變形，以流動形式破壞時，應該采用第三或第四強度理論。壓力容器上用第三強度理論，其它多用第四強度理論。第三強度理論認為最大剪應力是引起流動破壞的主要原因,如低碳鋼拉伸時在與軸線成45度的截面上發生最大剪應力,材料沿著這個平面發生滑移,出現滑移線。這一理論比較好的解釋了塑性材料出現塑性變形的現象。形式簡單,但結果偏于安全。）
空間中的應力分量：
三維空間中任意一點應力有6個分量：
對應于S11、S22、S33、S12、S13、S23（ABAQUS一般X作為1軸，Y當成2軸，Z是3軸），那么：
S11就是X軸向的應力，正值為拉應力，負值為壓應力；
S22就是Y軸向的應力，正值為拉應力，負值為壓應力；
S33就是Z軸向的應力，正值為拉應力，負值為壓應力；
S12就是在XY平面上，沿Z向的剪力；
S13就是在XZ平面上，沿Y向的剪力；
S23就是在YZ平面上，沿向的剪力；
其中S11=S21,S13=S31,S23=S32（剪力的對稱性）。
區別于主應力：
一般情況下,通過該點的任意截面上有正應力及其剪應力作用，但有一些特殊截面,在這些截面上僅有正應力作用,而無剪應力作用。稱這些無剪應力作用的面為主截面,其上的正應力為主應力,主截面的法線叫主軸,主截面為互相正交。主應力分別以σ1，σ2，σ3表示，按數值排序（考慮正負）為：σ1≥σ2≥σ3。
展開 mises屈服條件。由于pi平面的法線垂直于八面體的一個面，此面上的八面體剪應力才會導致屈服，因此也可用八面體剪應力建立屈服條件，實質一樣，因此也稱最大八面體剪應力屈服條件。當然還有其它解釋。
mises應力是一個等效應力，注意等效的意義。在單拉時其值就等于拉伸應力。
mises的相關專題、標簽、搜索
mises的最新內容
Compute the equivalent Von Mises stress !
Calculate equivalent von Mises stress
Smises=(stress(1)-stress(2))**2+(stress(2)-stress(3))**2
1 +(stress(3)-stress(1))**2
do i=4,ntens
Smises=Smises+6.d0*stress(i)**2
end do
原始文獻：《Mechanical modelling of indentation-induced densification in amorphous silica》
該文章為了模擬非晶態二氧化硅的壓縮力學性能，把拉伸與壓縮分開處理：拉伸側采用熟悉的 von Mises 屈服，壓縮側則切換到 cap 屈服面。
）
Insert → Stress → Equivalent (von-Mises)
評估最大應力位置（注意是否出現應力奇異）
9.4 間隙變化判斷（變形 > 0.25 mm 區域）
使用Insert → Expression或User Defined Result
公式：abs(UY) > 0.25顯示為 1（需在 Mechanical 中通過插值或閾值圖實現
;">本資源包含一份 PDF 文檔和可直接編譯運行的 Fortran UMAT 代碼，具體內容為：</span></p><p class="ql-align-justify">理想彈塑性本構 + 隱式積分 + 徑向返回</p><p class="ql-align-justify">完整公式推導 + Fortran 源碼直接編譯</p><p class="ql-align-justify">von Mises
喵星人認為以下幾種應力相關不變量相對比較重要：
Mises：基于第四強度理論，用戶手冊定義如下：
Tresca：基于第三強度理論，用戶手冊定義如下：
Tresca equivalent stress, defined as the maximum difference between principal stresses.
常規性能分析
**(1) 靜態強度/剛度分析**
- 求解：應力（Von Mises）、應變、位移。
- 目標：最大應力 < 材料許用應力（鑄鐵~80–120MPa，鑄鋁~150–200MPa）。
- 識別高應力區（缸孔周邊、法蘭、軸承座）。
同樣使用約束表面自由度的方式查看導入的von mises應力，方法二 穩定很多。
Inistate，read命令使用時的地址部分需要注意的是：模塊C4：計算寫出的file.ist文件不要直接復制到D模塊的計算文件夾。
這里在反過頭來說如何獲得符合彎折預期的初始應力。
第四強度理論：我們最常用的Von mises應力（畸變能密度理論），適用絕大多數塑性金屬材料的失效評估。
公式為：
而對于各向異性的塑料材質這四種理論顯然就不在適用了，那么我們怎么判斷這類塑料材質的應力仿真結果是否滿足強度要求呢。
采用了一個率無關的 von Mises 彈塑性材料模型，其中屈服應力為 700 MPa，硬化斜率為 0.3 GPa。命名為METAL。剛性墩頭不需賦予材料。
Section截面屬性
賦予材料屬性
創建剛性墩頭2D-Wire
裝配模型
將剛性墩頭與坯料裝配到一起，剛性墩頭與坯料上邊線重合。
