mises的案例
Matlab“稍作修改”Abaqus-odb結果 | 提取三維裝配體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掃頻曲線和應力幅值 ￥10
問題：
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應力選型，只能按三個方向來分別查看。
展開 使用python二次開發解決abaqus的mises曲線繪制 ￥9.9
熟悉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是一種屈服準則，準則的值我們通常叫等效應力，習慣稱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應力
mises屈服條件。由于pi平面的法線垂直于八面體的一個面，此面上的八面體剪應力才會導致屈服，因此也可用八面體剪應力建立屈服條件，實質一樣，因此也稱最大八面體剪應力屈服條件。當然還有其它解釋。
mises應力是一個等效應力，注意等效的意義。在單拉時其值就等于拉伸應力。
關于mises應力
mises屈服條件。由于pi平面的法線垂直于八面體的一個面，此面上的八面體剪應力才會導致屈服，因此也可用八面體剪應力建立屈服條件，實質一樣，因此也稱最大八面體剪應力屈服條件。當然還有其它解釋。
mises應力是一個等效應力，注意等效的意義。在單拉時其值就等于拉伸應力。
關于Von Mises ，哥們你郵箱滿了不能回復你，我發這里了希望你能看到
von_mises等效應力.doc
mises本構模型UMAT（附源代碼和詳細注釋）
STRESS
C
SMISES=(STRESS(1)-STRESS(2))*(STRESS(1)-STRESS(2)) +
1 (STRESS(2)-STRESS(3))*(STRESS(2)-STRESS(3)) +
1 (STRESS(3)-STRESS(1))*(STRESS(3)-STRESS(1))
DO 90 K1=NDI+1,NTENS
SMISES=SMISES+SIX*STRESS(K1)*STRESS(K1)
90 CONTINUE
SMISES=SQRT(SMISES/TWO)
C 計算mises應力
C HARDENING CURVE, GET YIELD STRESS
C 等效塑性應變和硬化曲線確定硬化后的屈服應力,初始等效塑性應變為0,
C 對應的屈服應力為初始屈服應力
NVALUE=NPROPS/2-1
CALL AHARD(SYIEL0,HARD,EQPLAS,PROPS(3),NVALUE)
C
C DETERMINE IF ACTIVELY YIELDING
C
IF (SMISES.GT.(1.0+TOLER)*SYIEL0) THEN
C 如果Mises應力超出了屈服應力,判斷發生了屈服,進入了塑性階段
C FLOW DIRECTION
C
SHYDRO=(STRESS(1)+STRESS(2)+STRESS(3))/THREE
C SHYDRO
展開 梁單元結構建模optistruct求解查看應力，沒有Von mises、normal stress? ￥20
本帖子是關于：整體以梁單元結構建模進行預應力模態分析，optistruct求解后查看應力結果，沒有von mises stress、normal/shear stress應力信息的原因，以及如何解決這個問題的方法。
前段時間接觸到桁架橋的結構分析，桿件橫截面主要為BOX和C型槽，C型槽的剪切中心和中性軸不重合，前處理采用梁單元cbeam建模，單元類型選擇cbar還是cbeam，可以參考：【HyperMesh寶典】之梁單元 (qq.com)。建立梁單元截面類型選擇HYPER BEAM庫下的thinwalled box和standard channel，屬性卡片選擇pbeam，求解后，hyperview查看應力結果發現只有element stress1D(s)下的CBAR/CBEAM Axial stress和long stress，沒有von mises stress、normal stress等應力。
網上搜索了一圈都沒有找到相關的問題的解決方法，也可能是我沒找全面，只能老老實實啃幫助文件，找到了關于Stress Result Written in HyperView，附上鏈接以及截圖：Stress Results Written in HyperView .h3d Format (altair.com)
展開 abaqus系列技巧17：后處理的mises SP：1是什么鬼
接下來就是我們繪制節點應力的時候，有時候會出現這樣的情況，如圖1，圖2所示
這時候我們在一個節點處發現了兩條曲線：S:MISES SP:1; ：S:MISES SP:7.這個和我們的預期有些出入，，怎么會有兩條曲線？明明我只選擇了一個點啊。
別急，讓我們來梳理下：
首先，我們用的是殼單元，注意畫重點，是殼單元：
來看下殼單元的描述，在殼單元中，有個重要的概念是積分點，即雖然是一個面，但通過積分點，可以等效成一個實體。同時一個截面可以設置不同的積分點數量，如在本個案例中，設置了7個
在abaqus中，會自動輸出最上面和最下面一個積分點應力值，所以就是SP1和SP7.
如果我講到這里，你就認為了解了，你真的不太適合學有限元，學有限元要舉一反三，力求深解。其中最主要的問題就是SP1和SP7的具體位置是殼的哪個面。
SP代表什么意思
我們在MESH模塊中用query的功能查看單元的法向，會得到下面的注解，brown的方向是正向positive。所以SP = shell positive 也可以看出SPOS。同樣的SNEG= shell negative，就是負向。
接著我們看之前發過的一個圖
殼的頂面就是SP ，那么SP1就是底板，SP5就是頂面。對于本例題，就是SP7是頂面。這可能和我們的認識有出入。那么我們實例驗證下
一個實體單元，一個殼單元，殼單元，Z正向為法向方向，可以參見上面的圖
結果查看：
取殼單元上一點，輸出應力
其中sp1＞sp7.對照實體單元的云圖，Z向負向的面應力大，結果和預期一致。
想獲得幻想飛翔最新CAE技術文章，請關注幻想飛翔公眾賬號：幻想飛翔CAE。
想獲得幻想飛翔最新CAE技術文章，請關注幻想飛翔公眾賬號：幻想飛翔CAE。
展開 各向同性硬化von Mises率無關彈塑性本構理論以及umat源代碼 ￥99
各向同性硬化von Mises率無關彈塑性本構理論以及umat源代碼
1 本構理論
1.1 率形式
對于各向同性線彈性材料，其本構方程為：
式中假設了應變張量可以分解為彈性應變和塑性應變兩部分：
因此塑性本構的關鍵在于計算塑性應變的演化。對于率無關彈塑性的本構理論，需要確定以下三個部分：
(1)：屈服條件
(2)：流動法則
(3)：硬化法則
在此采用的是 von Mises 屈服條件：
式中后繼屈服應力是等效塑性應變的函數：
流動法則為：
式中流動方向的表達式為：
硬化法則為：
1.2 Return-mapping算法
上述的本構方程均為率形式。在增量步中，給定增量應變：
首先假設該增量應變全為彈性應變，計算試驗狀態下的一些物理量：
試驗狀態下的應力
試驗狀態下的屈服函數值：
利用該試驗屈服函數值來判斷在該增量步下是否發生了塑性屈服。如果：
則說明試驗狀態即為真實狀態，即可進行更新：
反之則需要進行塑性更正，即需要計算塑性乘子的增量，利用以下非線性方程組進行計算：
可以將該非線性方程組簡化至一個非線性方程，過程如下，將該方程組中的第一式分解為球量和偏量兩部分：
因此可以計算應力為：
將上式中的第二式整理得到：
可以得到兩個張量的方向相同：
因此偏應力可以用試驗狀態的信息表示出來：
代入到最后一個一致性方程中可得：
即可利用牛頓迭代法對上述非線性方程進行求解，得到塑性乘子增量。
展開 
馮-米塞斯(von Mises)塑性準則小白案例 ￥10
這個案例是一個驗證馮-米塞斯(von Mises)塑性準則小白案例,結構是個立方體1*1*1，主要是對本構模型的一個回顧。因為7、8年沒怎么摸過了（之后想用abaqus編一些自己的本構模型），所以先做了這么一個驗證模型。付費內容里面有完整的cae操作流程，并給了流程的inp文件，可以下載。
1. Part 創建 （3D）
Next:
點擊Rectangle 分別輸入坐標（0，0）和（1，1），
正方形建好后重新點，底部點done
Next：Depth 改為1，然后點ok，出現1*1*1的立方體
拓撲優化結構MISES應力分布3D顯示MATLAB代碼 ￥100
拓撲優化結構MISES應力分布3D顯示MATLAB代碼
用python保存odb文件中每一幀的mises應變
直接在abaqus里面調用就可以了
from abaqus import *
from abaqusConstants import *
from odbAccess import *
import numpy as np
import os
SavePath = 'G:/Abaqus_jobFile/6/'
odb = openOdb(path='G:/Abaqus_jobFile/6.odb ')
for i in range(len(odb.steps['static'].frames)):
frame = odb.steps['static'].frames[i]
#print frame.frameId
fileName = SavePath + str(frame.frameId)+'.txt'
with open(fileName,'a') as f:
s_field = frame.fieldOutputs['S'].values
for j in s_field:
data = str(j.elementLabel)+'\t' +str(j.mises) + '\n'
f.write(data)
展開 ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值
ABAQUS中mises應力云圖顯示的最大值還不到屈服應力值為啥還有PEEQ值，PEEQ云圖有變形值
