
發布
注冊
/
登錄玻璃在abaqus中參數的案例
Abaqus中材料參數隨機場實現
3、隨機參數的賦值
在abaqus中,前處理都是利用python進行的,且在CAE中每一個操作,都會在abaqus的工作目錄中的abaqus.rpy中生成相對應的代碼。不熟悉abaqus某一步操作代碼時候,都可以在cae中進行一步操作,然后去看abaqus.rpy中生成的最新的一句代碼。另外,還可以使用file->macro manager錄制相應cae中操作的宏代碼,錄制完成后會在工作目錄下面生成abaqusMacros.py文件,里面是相應操作的python代碼。Abaqus中運行python代碼有兩種方式,第一種是在abaqus下面的輸入框中輸入代碼,如下圖;還有一種方法是在abaqusCAE中,file->run script。
展開 abaqus中參數化BCC點陣 ¥5
點陣機構參數(例如桿件直徑、單胞大小等)較多,當修改參數時建模工作量較大,因此作者基于Abaqus開發了BCC點陣結構的快速建模插件。
點陣結構有兩種模型形式:一種是用線表示,該方法單元數量少,便于計算;另一種是實體模型,網格數量大,計算費時,主要用于輸出3D打印的幾何模型。圖1和圖2是改變兩種參數后得到的點陣結構
圖1 553實體點陣結構
圖2 333 實體點陣結構
圖3 梁單元點陣結構
操作過程如圖4所示,在abaqus標題欄選擇 File > Script讀取后綴為.py的腳本文件。在彈出的窗口中分別輸入參數。
圖4 腳本運行過程
——————————————————————————————————————————————————————————————————————————————————
注:附件中為txt文本,下載后將后綴改為.py即可。
附件中BCCBeam.txt為梁單元模型
附件中的solidBCC.txt為實體單元模型,如圖4
展開 利用Isight/DOE辨識Abaqus分析模型中的關鍵參數
Isight中的試驗設計方法DOE模塊用于辨識關鍵參數、構件經驗公式以獲得最佳設計。DOE運用包括三個步驟:試驗計劃、執行試驗和結果分析。
Isight/DOE與Abaqus結合,可以將Abaqus分析模型中的參數進行試驗設計,辨識出關鍵參數,從后處理Pareto圖等中分析出參數對輸出響應的影響大小。以沖擊力作用的鋼架結構為例進行說明,其中以兩種材料Al和STEEL的彈性模量和泊松比為設計參數,試驗測得鋼架結構在沖擊力作用下的位移響應曲線與仿真結果曲線之差為輸出響應。下圖為Abaqus分析模型。
在Isight中建立Abaqus工作流,如下圖所示。
Abaqus工作流
Isight中,有9種DOE方法,如參數試驗、全因子設計、部分因子設計、正交數組、中心組合設計、Box-Behnken設計、拉丁超立方設計、優化拉丁超立方設計及自定義數據文件,用戶也可以進行二次開發,自編DOE方法。用戶可以根據自己的需求選擇合適的試驗設計方法。
Isight提供后處理工具,幫助進行試驗設計的后處理工作,有試驗數據表格、散點圖、ANOVA分析表、Pareto圖、主效應圖、交互效應圖及相關性圖等。下圖為本鋼件結構的后處理。
利用IsightDOE辨識Abaqus分析模型中的關鍵參數.pdf
展開 Abaqus定義幅值曲線(Amplitude)中的參數設置
而在Abaqus/Explicit中則沒有默認進行光滑處理,除非是繼承自相關的有限時間增量步。這個參數僅僅當需要時間微分(對位移或者速度邊界條件在直接積分動力分析中)且忽略選項中其它的使用才應用。
smoothing所指定數值的意義解釋見下文。
通常通過表格定義的幅值函數為分段線性的,當指定了Smoothing的值時,就確定了在函數尖點兩端進行光滑處理的鄰域區間長度。可見下圖:
圖3 分段線性函數的光滑化處理比較(位移、速度和加速度)
即,在尖點的兩側各指定了一段區間長度:
t=用戶指定的smoothing值*尖點兩側線性區間長度t1、t2的較小值
在這個長度為2t的區間內進行二次插值。圖3中顯示了對位移幅值分別求一階導數(速度)、二階導數(加速度)的不同效果。
取值規定:
Smoothing值的范圍是[0,0.5]。
Use solver default選項為,對于Abaqus/Standard取0.25,對于Abaqus/Explicit取0.0;
Specify則需用戶自行指定數值(在包含較大時間間隔(time interval)的幅值曲線中,建議取0.05,以避免嚴重偏離用戶的幅值定義)。
轉自公眾號——ABAQUS大世界
旨在分享,若侵即刪.
展開 
[原創]Abaqus中材料參數不為常數的Python處理方法
[原創]Abaqus中材料參數不為常數的Python處理方法
在有限元計算中,可能會碰到材料參數(這里以線彈性材料彈性模量E為例)在模型中并不為常值的情況,例如,E與某一坐標相關。大家知道,abaqus材料中材料是在*section對截面屬性定義的時候對單元組進行賦值的。并且GUI中也沒有提供變參數的接口(這里順路提一下,載荷是可以通過函數給出變載荷的,GUI中有接口)。這就需要我們采取一些特殊的措施。在CAE-lab前文《Abaqus中計材料參數不為常數的兩種處理方法》中, 給大家介紹了‘編輯inp文件’和‘使用Umat’兩種方法,這里就給大家介紹一下用python腳本實現的方法。
方法
要求
難度
編輯inp
會一門編程語言,能夠完成對文本文件的讀取和寫入
4星
使用Umat
會使用abaqus的Umat用戶子程序,計算文件中不需要使用其它的用戶子程序
3星(如果對abaqus用戶子程序不熟悉,則為5星)
Python腳本
會python基礎,對python編寫abaqus腳本了解
3星(如果對python腳本不熟悉,則為5星)
Python大法在abaqus前后處理中具有舉足輕重地位和天然的優勢,通過python給單元賦材料屬性非常簡單,本例中,假設每個單元的彈性模量不一樣,并且是一個函數,我們的思路就是遍歷所有的單元,給每一個單元定義一種材料。具體實現如下:
1.模型
以簡單的Cube模型為例,如圖1,注意由于我們是遍歷單元賦材料,Abaqus中模型必須在劃分好網格后在進行賦材料屬性。
展開 ABAQUS中混凝土CDP模型的參數計算(下附python程序代碼) ¥19
<p>因為要仿真混凝土破壞實驗,考慮用abaqus里面的CDP模型,查閱了相關資料進行了理論總結,并根據理論編寫計算程序。</p><p>ABAQUS中CDP 模型中采用的是混凝土在單軸受力狀態下的應力和非彈性應變,非彈性應變根據混凝土的單軸應力-應變曲線換算。
CDP模型參數在ABAQUS中應用-小白案例(附inp和案例cae文件) ¥9
這個帖子的重點放在cdp模型參數的測試上,所以在abaqus中建立一個單位立方體進行計算,得到壓應力應變如下:
立方體大小是1*1*1。
如何在abaqus建立方體在前面一個帖子中寫過,在此不再重復。Cdp模型參數如何計算在上一篇帖子中詳細說明,在此直接拿過來用。
1、 材料設置,
1.首先設置彈性參數:
2.再設置塑性參數,菜單欄里找到Mechanical->Plasticity->Concrete Damaged Plasticity,設置如下參數,可微調:
展開 ABAQUS UEL中UMAT材料參數設置問題,等效迭代數(Equll Iter)的含義
在UMAT設置中,一般來說是設置楊氏模量為1E-11,泊松比為0.3。請問一下為什么要這樣設置呢?既然只是借用實體單元的積分點來傳遞 UEL計算得到的SDV,通過一個common來傳遞,那這個參數為什么要設置這么小?
我最近做了一些工作,發現如果把楊氏模量設置為實際的材料參數,得到的結果會和實際情況相比偏小,設置為1E-11的話,結果就基本一致。
問題二:implicit算法是需要設置時間步長,一般有初始值,最大值,最小值;也可以設置為固定步長。我最近的工作設置為了固定步長,在作業監控器界面出現了下圖的情況,沒有不連續的迭代,但是等效迭代欄(Equll Iter)出現了大于1的情況,請問一下大家,這個數代表了什么?我的理解是第一次嘗試計算不收斂,然后嘗試第二次,但是我設置了固定步長,那么第二次為什么會出現收斂?
展開 ABAQUS中材料參數不為常數的兩種處理辦法——Part1
在有限元計算中,有時候會碰到材料參數在離散區域內并不是常數的情況,而有限元軟件中一般也沒有給出相應的接口。本文通過兩種方法,即編輯inp文件和使用Umat的方法來實現材料參數與坐標相關的計算。其中,Part1 介紹了兩種方法的需求和基本思路,并就編輯inp方法做了介紹。Part2 介紹Umat方法和算例。由于時間關系,這里先更新Part1 下次更新Part2。
Part 1
1目的
在有限元計算中,可能會碰到材料參數(這里以線彈性材料彈性模量E為例)在模型中并不為常值的情況,例如,E與某一坐標相關。大家知道,abaqus材料中材料是在*section對截面屬性定義的時候對單元組進行賦值的。并且GUI中也沒有提供變參數的接口(這里順路提一下,載荷是可以通過函數給出變載荷的,GUI中有接口)。這就需要我們采取一些特殊的措施。本教程給出2種行之有效的方法,它們分別是:(1)編輯inp文件和(2)使用Umat。
2方法
方法
要求
難度
編輯inp
會一門編程語言,能夠完成對文本文件的讀取和寫入
4星
使用Umat
會使用abaqus的Umat用戶子程序,計算文件中不需要使用其它的用戶子程序
3星(如果對abaqus用戶子程序不熟悉,則為5星)
上表中給出兩種方法的一些基本要求,編輯inp文件的方法要求更低,它通過直接編輯inp文件的方法來解決問題。而使用Umat的方法要求更高,需要用戶對子程序有一定的了解,并且也需要會一門編程語言。
展開 【JC本構插件】abaqus中如何確定Johson-Cook本構A、B和n等參數 ¥19.89
當我們用abaqus模擬沖擊動力學問題時,經常會考慮使用Johson-Cook本構,而正確輸入材料本構的各參數,對我們的仿真結果意義重大,今天我們就來介紹下abaqus中JC本構的各參數識別問題。
Johnson-Cook塑性模型是一種具有硬化規律和速率依賴的解析形式的米塞斯塑性模型,主要適用于許多材料的高應變率變形模擬,包括大多數金屬。
通常用于絕熱瞬態動態模擬;與Abaqus/Explicit中的Johnson-Cook動態失效模型結合使用;Abaqus/Explicit中,可以結合拉伸破壞模型來模擬拉伸剝落或壓力斷口;可與漸進損傷和失效模型(漸進損傷和失效)結合使用,以指定不同的損傷起始準則和損傷演化規律,同時允許材料剛度的漸進退化和網格單元的移除;必須與線彈性材料模型(線性彈性行為)或狀態方程材料模型(狀態方程)結合使用。
下面是JC本構的一般表達式,該模型中主要確定A、B、n、C和m等參數。可以看到J-C本構的主體由三部分構成,分別表征了材料的應變硬化、應變速率硬化(強化)以及溫度軟化,可以概括為“兩硬一軟”。
A-參考應變率和參考溫度下的初始屈服應力,B和n-材料應變硬化模量和硬化指數,C-材料應變率強化參數,m-材料熱軟化指數。
查幫助文檔可以知道各參數含義如下:
當我們不考慮應變速率和溫度影響時,該表達式就簡化為下面的表達式:
如果我們確定了參數A、B和n,那么我們在abaqus中就能輸入相應的JC參數,重點來了!
展開