各向同性
關注創建者：匿名 創建時間：2021-12-24
各向同性的視頻教程
Abaqus-umat-各向同性線彈性umat
一、視頻內容介紹 二、fortran語言基本概念 ? ? # fortran77固定格式 ? ? # IDE——vscode+Modern Fortran 三、各向同性線彈性本構 ? ? # 各向同性線彈性本構 四、umat程序編寫及應用 ? ? # umat函數頭 ? ? # umat需要實現的功能 ? ? # 線彈性umat編程實例講解 ? ? # umat與abaqus內置線彈性模型對比
免費 56分鐘 3164播放
查看
Abaqus-VUMAT-各向同性硬化彈塑性vumat
Abaqus-VUMAT-各向同性硬化彈塑性vumat，手把手教學，包含本構、離散、編程實例詳解、調用。
￥79 3小時19分鐘 2460播放
查看
各向同性的實例教程
各向同性，橫觀各向同性，正交各向異性三種線彈性umat程序
1 各向同性
各向同性線彈性材料的彈性矩陣為：
式中拉梅常數的表達式為：
因此在編寫各向同性材料的umat時，需要兩個材料參數，在這里我們使用楊氏模量E和泊松比v。
2 橫觀各向同性
橫觀各向同性線彈性材料的彈性矩陣為：
并有關系式：
可見其彈性矩陣需要5個獨立的參數，為下列5個工程常數：
下標a代表軸向，下標t代表橫向。
3 正交各向異性
正交各向異性線彈性材料的彈性矩陣為：
并有關系式：
因此對于正交各向異性材料，其彈性矩陣需要9個工程常數來確定：
4 程序
使用Fortran90編寫umat程序。由于Abaqus默認的umat子程序為Fortran77，因此為了使用f90程序，使用命令：
abaqus make library=xxx.f90
該命令可以生成相應的后綴為obj的文件，之后使用該文件即可。使用上述方法可以避免使用Fortran77進行umat的編寫。
展開 Chanboche模型是一種用于描述材料各向同性非線性隨動硬化行為的材料本構模型。該模型由Chanboche在1981年提出，其基本形式包括各向同性部分和隨動硬化本構部分。
具體而言，Chanboche模型各向同性本構部分可以用以下方程表示：
dR(p)=b(Q-R)dp
非線性隨動硬化模型可以用以下方程表示：
dx=(2/3)cdεp-rxdp
本程序已經在上一個帖子基礎上進一步完善，實現可直接輸入試驗拉伸循環曲線，計算本構參數，黑色線為計算結果，紅色為試驗循環拉伸應力應變曲線。
<p>Chanboche模型是一種用于描述材料各向同性非線性隨動硬化行為的材料本構模型。該模型由Chanboche在1981年提出，其基本形式包括各向同性部分和隨動硬化本構部分。</p><p>具體而言，Chanboche模型各向同性本構部分可以用以下方程表示：</p><p>dR(p)=b(Q-R)dp</p><p>非線性隨動硬化模型可以用以下方程表示：</p><p>dx=(2/3)cdεp-rxdp</p><p>程序基于3個背應力分量編寫，效果參見鏈接<a href="https://www.bilibili.com/video/BV1B54y1F7gS/?vd_source=9f1dda2358e63ace0b661e56fe417806" rel="noopener noreferrer" target="_blank">https://www.bilibili.com/video/BV1B54y1F7gS/?vd_source=9f1dda2358e63ace0b661e56fe417806</a>，程序為回映算法核心算法，可以修改此程序實現基于試驗數據的本構參數計算，不太會編程的可移步我的另外一個帖子，具體的<a href="https://www.yqgqt.org.cn/major/matlab" rel="noopener noreferrer" target="_blank">matlab程序</a>如下：</p>
展開 在這種情況這些值大致都是從各向同性多孔催化劑得到的。
1、問題描述
本案例使用已經畫好的體網格，導入以后的網格如下圖。
2、STAR-CCM+設置
本案例的介質是空氣，模擬狀態是穩態、湍流且不可壓縮。使用標準的K-Epsilon 模型。
（1）選擇反應類型相應的湍流模型；
（2）設置邊界條件，在Regions > Fluid > Boundaries > inlet > Physics Conditions > Turbulence Specification節點，將Method property改為Intensity + Length Scale。將進口的邊界條件改為下表：
（3）設置多孔介質。在Regions > Porous節點，將類型改為Porous Region.。在Physics Values > Porous Inertial Resistance節點，將Method property改為各向同性Principal Tensor。在Porous Inertial Resistance > Principal Tensor節點，同時選擇XX Component，YY Component和ZZ Component，將多孔慣性阻力改為25 kg/m4。在Porous Viscous Resistance > Principal Tensor節點，把XX Component，YY Component和ZZ Component改為1500 kg/m3-s計算后處理。這樣就設置好了各向同性多孔介質的多孔慣性阻力和粘性阻力。
3、計算界面的速度適量
計算域速度適量分布
多孔介質的阻力壓降
多孔介質的質量流量
本文轉自有限猿仿真博客，感謝原作者。如有侵權請立即聯系刪除。
展開 各向同性硬化von Mises率無關彈塑性本構理論以及umat源代碼
1 本構理論
1.1 率形式
對于各向同性線彈性材料，其本構方程為：
式中假設了應變張量可以分解為彈性應變和塑性應變兩部分：
因此塑性本構的關鍵在于計算塑性應變的演化。對于率無關彈塑性的本構理論，需要確定以下三個部分：
(1)：屈服條件
(2)：流動法則
(3)：硬化法則
在此采用的是 von Mises 屈服條件：
式中后繼屈服應力是等效塑性應變的函數：
流動法則為：
式中流動方向的表達式為：
硬化法則為：
1.2 Return-mapping算法
上述的本構方程均為率形式。在增量步中，給定增量應變：
首先假設該增量應變全為彈性應變，計算試驗狀態下的一些物理量：
試驗狀態下的應力
試驗狀態下的屈服函數值：
利用該試驗屈服函數值來判斷在該增量步下是否發生了塑性屈服。如果：
則說明試驗狀態即為真實狀態，即可進行更新：
反之則需要進行塑性更正，即需要計算塑性乘子的增量，利用以下非線性方程組進行計算：
可以將該非線性方程組簡化至一個非線性方程，過程如下，將該方程組中的第一式分解為球量和偏量兩部分：
因此可以計算應力為：
將上式中的第二式整理得到：
可以得到兩個張量的方向相同：
因此偏應力可以用試驗狀態的信息表示出來：
代入到最后一個一致性方程中可得：
即可利用牛頓迭代法對上述非線性方程進行求解，得到塑性乘子增量。
展開 
各向同性的相關專題、標簽、搜索
各向同性的最新內容
這種金剛石晶格微觀結構也導致各向同性的材料行為。
圖6. 金剛石晶格結構的工程常數
案例4：編織結構（布料）
11. 按照案例1的相同步驟操作。為紗線定義各向同性材料屬性（E=1000MPa, ν=0.3），為基體定義各向同性材料屬性（E=100 MPa, ν =0.3）。
12. 使用默認幾何設置定義編織結構RVE（圖7）。生成網格。
傳統有限元分析通常采用各向同性塑性模型，通過宏觀應力–應變曲線描述材料響應。但實際金屬材料并不是“均勻黑箱”：晶粒取向、滑移系激活、織構演化都會影響局部塑性變形，尤其在薄壁管壓潰這類大變形、強局部化問題中，微觀結構可能對吸能行為產生重要影響。
第四，作者還加入了一個各向同性 accommodation 項，用來描述晶界附近非晶體學協調變形的作用。這個處理非常值得重視。很多時候我們做晶體塑性，只把目光放在晶內滑移和孿晶上，卻容易忽略多晶材料中晶粒之間并不是天然完全協調的。Staroselsky 這篇文章清楚地認識到：如果不考慮這部分效應，數值計算中的應力水平會偏高，甚至難以合理匹配實驗。
摘要
分層介質組件旨在對一系列平面圖層進行嚴格而快速的分析，其中每個平面圖層后面都是均質(各向同性或各向異性)介質。這種配置在例如涂層應用中特別令人感興趣。在這個用例中，我們展示了如何在VirtualLab Fusion中定義這樣的結構，并深入探討了它的特性。
在哪里可以找到組件？
課程仔細解釋了MPPIC的構建模塊，包括平均方法、阻尼模型、各向同性修正和堆積模型。每個概念都與其在OpenFOAM中的實現相關聯，您還將檢查源代碼以了解這些模型是如何構建的。這種對代碼庫的深入研究對于想要擴展或定制求解器以用于自己的研究或工業應用的人來說特別有價值。
在整個課程中，重點放在實踐學習上。您將從零開始設置多個模擬，修改字典，運行求解器并可視化結果。
它由各向同性材料和雙折射材料交替堆疊而成。
折射率根據以下公式呈現色散：
中??1,2，??1,2和??1,2是擬合參數。
材料1和2的厚度??1和??2可由以下公式推導得出：
在這個例子中我們考慮：
和
步驟1. 對具有均勻層的反射式偏振片運行仿真并導出計算結果
在此步驟中，通過掃描入射角（θ和φ）來評估反射偏振片的反射特性。
建模任務
非球面鏡和準直物鏡
透鏡系統元件允許簡單地定義一個由光滑表面和均勻的、各向同性的介質的交替排列組成的元件。就界面和材料而言，可以從內置的目錄中選擇已有的選項，或者定制您自己的界面和材料，以獲得最大的靈活性。
同時我們可以假設纖維增強塑料是一種特殊的各向異性材料，在垂直纖維方向的平面內材料又是各向同性的。這樣Hill材料常數H、F、G、N、L、M的計算，就由、六個測試數據，變為=四個數據。
通常我們是可以查到PA基體的力學參數（拉伸屈服強度）和PA+GF20 的拉伸屈服強度。
在各向同性的情況下，使用一個非常小的周期，以確保只有0階會傳播。二氧化硅層也是根據參考文獻來定義的。
- 涂層厚度：10納米
- 涂層材料。二氧化硅
- 折射率：擴展的Cauchy模型。
?? = 1.44, ?? = 0.00422????2, ?? = 1.89?? - 05????4
- 基板材料：晶體硅
- 入射角度。
3 階 SH 是方向輻射的低通濾波器，各階對應從各向同性基礎顏色（0 階）到較高頻高光（3 階）的外觀描述。
需要注意的一個關鍵是，SH 系數無法被解碼為 IOR、粗糙度、金屬度等物理參數，它將光照與材質混合烘焙（baked），兩者無法分離。這并非工程限制，而是數學層面的根本局限——SH 的優化目標是最小化渲染誤差，而非最大化物理可解釋性。
