abaqus木材分析的案例
利用Abaqus的UMAT子程序仿真木材蠕變現象
下面我們將以木材蠕變為例，介紹下在ABAQUS UMAT中如何實現蠕變仿真。
2. 本構理論
文獻[1]給出了木材蠕變過程中本構：
3. 算例
3.1 模型
考慮懸臂梁模型，如下圖。
3.2 邊界條件
根據蠕變的定義，模型必須現有一個穩定的載荷，因此可以分成兩個step。第一個step，完成力加載，第二個step保持載荷，實現蠕變變形的生長。
需要指出的是，蠕變通常需要在較長的時間尺度上才能有明顯的效果。比如我們要觀察100天的變形情況，那么這個時候，ABAQUS設定的總時間還是1，在UMAT里面要乘以相應的系數，給出物理時間，才能有效的實現蠕變效果。
3.3 結果
最終得到懸臂梁端部位移如下圖所示。從圖中可以看出，在最初的幾天，蠕變變形較大，隨著時間的推移，變形增加的幅度放緩，符合蠕變的特點。
參考文獻
[1] 《濕度變化和荷載作用下膠合木曲梁的工作性能研究》
[2] 《旋切板膠合木的蠕變及其對結構穩定性的影響》
最后，有abaqus二次開發相關需求歡迎大家通過我們的微信公眾號聯系我們。
微信公眾號：320科技工作室
展開 基于Hashin準則的木材受火后強度分析
木結構在火災作用下，木材的化學成分及物理特性會發生復雜的變化，導致其力學性能的改變，木構件承載能力降低。本文編寫了考慮溫度影響的木材本構vumat子程序，并對子程序的正確性進行了驗證。
木材是一種復雜的各向異性復合材料，可以分為三個方向，即縱向順紋(L)、徑向橫紋(R)、弦向橫紋(T)。其復雜的本構關系主要體現為在拉或剪力作用下發生 脆性破壞，而在壓力作用下發生塑性變形，且在橫紋壓力作用下變形較大，同時拉壓強度不相等。
屈服準則
木材是各向異性材料，且L、R、T三個方向的拉壓屈服強度不一樣，屬于拉壓非對稱材料。為了準確地預測木材的失效需要選擇合適的各向異性屈服準則，目前常用的各向異性屈服準則有：Hill準則，Hosford準則，Yamada-Sun屈服準則等。本采用Hashin準則作為木材的屈服準則
損傷演化準則
本文木材本構關系模型定義了兩種不同的損傷演化模型，受壓延性破壞采用理想彈塑性模型，受拉脆性破壞采用線彈性軟化模型.引入損傷變量D來描述木材的受損狀態。
溫度影響
木材隨著溫度的升高，發生不同程度的炭化，其強度、彈性模量、斷裂能也隨之發生變化。本文考慮了溫度對木材的模量、強度的影響，并且認為溫度對拉壓性能產生的影響不同。
根據上述相關理論編寫了abaqus vumat子程序，并通過單胞模型對子程序進行驗證。
下圖為不同溫度下單向拉壓結果
下圖為三點彎曲載荷下的破壞行為
展開 考慮溫度影響的vumat子程序在木材受火后強度分析中的應用
木結構在火災作用下，木材的化學成分及物理特性會發生復雜的變化，導致其力學性能的改變，木構件承載能力降低。本文編寫了考慮溫度影響的木材本構vumat子程序，并對子程序的正確性進行了驗證。
木材是一種復雜的各向異性復合材料，可以分為三個方向，即縱向順紋(L)、徑向橫紋(R)、弦向橫紋(T)。其復雜的本構關系主要體現為在拉或剪力作用下發生 脆性破壞，而在壓力作用下發生塑性變形，且在橫紋壓力作用下變形較大，同時拉壓強度不相等。
屈服準則
木材是各向異性材料，且L、R、T三個方向的拉壓屈服強度不一樣，屬于拉壓非對稱材料。為了準確地預測木材的失效需要選擇合適的各向異性屈服準則，目前常用的各向異性屈服準則有：Hill準則，Hosford準則，Yamada-Sun屈服準則等。本采用Hashin準則作為木材的屈服準則
損傷演化準則
本文木材本構關系模型定義了兩種不同的損傷演化模型，受壓延性破壞采用理想彈塑性模型，受拉脆性破壞采用線彈性軟化模型.引入損傷變量D來描述木材的受損狀態。
溫度影響
木材隨著溫度的升高，發生不同程度的炭化，其強度、彈性模量、斷裂能也隨之發生變化。本文考慮了溫度對木材的模量、強度的影響，并且認為溫度對拉壓性能產生的影響不同。
根據上述相關理論編寫了abaqus vumat子程序，并通過單胞模型對子程序進行驗證。
下圖為不同溫度下單向拉壓結果
下圖為三點彎曲載荷下的破壞行為
最后，歡迎通過公眾號“320科技工作室”聯系我們
展開 Abaqus 非線性屈曲分析方法 附ABAQUS分析手冊分析卷下載
當然，對于方筒這類實際上是通過顯示方法實現的，更準確的講是動力屈曲分析，所以我們還得判斷動能、塑形耗散等能量參數，才能使結果更加準確。
下載地址：ABAQUS分析手冊分析卷
Abaqus子結構與子模型分析技術 附ABAQUS結構工程分析及實例詳解文檔下載
-通過2個工程案例學習Abaqus中的子結構與子模型分析技術”
子結構與子模型技術在Abaqus中屬于模擬抽象化的范疇，所有Abaqus模型都涉及一定程度的抽象，但是除了傳統有限元的抽象方法之外，還可以通過以下幾種模擬抽象化技術來降低求解成本。
子結構
子模型
生成矩陣
對稱模型生成、結果傳遞和循環對稱模型
周期介質分析
網格劃分的梁橫截面
擴展有限元方法（XFEM）
適當地利用這些抽象化建模技術可以極大地提高Abaqus的分析效率，本期文章介紹一下子結構和子模型技術。
01
—
子結構
在有限元分析里，子結構也叫超級單元，是由多個單元組成的一個“整體單元”，它在線性分析的基礎上消除了“整體單元”中保留節點以外所有節點的自由度；子結構的系統矩陣（剛度、質量）也被縮聚成較小的矩陣，可以根據需求恢復內部求解。
很多實際工程結構都比較龐大，導致完整結構的有限元模型計算量超出計算機的硬件資源，對于具有線性響應的此類問題，可以使用子結構縮聚的方法，在一般配置的計算機上來求解完整結構的響應。
展開 基于Abaqus的建筑結構隔震分析 附ABAQUS建筑結構分析應用下載
本文通過時程分析的方法，考察隔震結構在大震作用下的性能，結果顯示，在大震作用下，結構的整體響應，無論是位移角還是結構的剪力，與小震結果都有明顯差異，隔震支座對結構性能的改善，主要體現在結構的上部，對結構的中下部則較小，且不再滿足規范中對剪力降低50%的要求。另一方面，非線性的影響會對結構的計算結果起到放大作用，使微小差異的結構方案在大震作用中表現出明顯不同的抗震性能。
下載地址 ：ABAQUS建筑結構分析應用
Abaqus超彈性材料分析 附Abaqus 分析用戶手冊材料卷下載
三、后處理
1、位移云圖
圖8 位移云圖
2、應力云圖
圖9 接觸定義
下載地址：Abaqus 分析用戶手冊材料卷
Abaqus接觸非線性在有限元計算分析中的應用 附莊茁ABAQUS非線性有限元分析與實例下載
來源：有限元在線
ABAQUS的非線性主要在有三種：幾何非線性，材料非線性以及接觸非線性。接觸非線性在ABAQUS的有限元計算分析中應用非常廣泛，特別是動態顯式的求解，只要模型中包含兩個以上相互接觸的部件，就要用到接觸非線性。
ABAQUS接觸非線性的設置主要在Interation模塊中完成，設置接觸的屬性時，可以設置摩擦系數，阻尼系數，損壞，失效準則等非線性參數，如圖1所示。
如圖2所示，在接觸定義界面，可以選擇通用接觸、面-面接觸、自接觸等各種非線性接觸方式。
在接觸編輯界面，可以選擇機械約束方式為運動學接觸算法，或是懲罰接觸方式，還可選擇滑移方式為有限滑移或小滑移，如圖3所示。
這是對模型定義非線性接觸后得到的分析結果，以供參考。
下載地址：莊茁ABAQUS非線性有限元分析與實例
展開 ABAQUS的直齒圓柱齒輪模態有限元分析 附ABAQUS有限元分析常見問題解答下載
下載地址：ABAQUS有限元分析常見問題解答
Abaqus的響應譜分析 附Abaqus頻響分析完整過程下載
在ABAQUS中，響應譜分析是分為兩步完成的，第一步需要設置一個頻率提取分析步，提取結構的前幾階固有頻率；在第二個分析步中設置響應譜分析。
值得注意的是，譜分析的激勵是在step中加載的，不需要在load中進行設置。
下載地址：Abaqus頻響分析完整過程
abaqus電池包隨機振動疲勞分析（附模型及分析流程） ￥88
本例展示基于功率譜密度曲線（PSD）的電池組疲勞分析，即針對隨機振動的疲勞壽命 分析。
1 問題設定 一塊電池組，尺寸為 70mm x 175mm x 400mm。該電池組的兩端共有 6 個端點，分別受 到垂直于電池組平面的激勵作用，且激勵的加速度功率譜密度曲線（ASD）相同。 由于在隨機振動基于線性動力學原理，因此電池，PC 材料等采用實體建模，其他鈑金 采用殼單元建模， 設定相關的 fastener 點焊單元，coupling 耦合單元和 tie 約束，建立零件 和零件之間相應的連接關系。
兩端所對應的 PSD 譜線如下圖。請注意該曲線的頻率截斷在 200Hz 處。
2 分析過程 一般來說，針對隨機振動的疲勞分析包含兩大步。第一步是在 Abaqus 中完成固有模態 和掃頻兩個計算；第二步是把這兩個計算結果與 PSD 曲線一起輸入 fe-safe，運行若干設置 后完成疲勞分析，得到相關結果。
以下內容包含完整的詳細的電池包跌落仿真分析 附件為完整教程和CAE模型文件.rar
展開 
Abaqus預應力模態分析 附Abaqus 分析用戶手冊材料卷下載
預應力模態
模態分析是一個線性攝動分析，只能進行線性求解。在動力學方程中，其載荷矩陣和阻尼矩陣為0，特征值的提取只取決于剛度矩陣和質量矩陣。而結構在外載荷的作用下剛度矩陣會發生變化，也就間接影響了結構的固有頻率。而預應力狀態下，我們不清楚剛度矩陣的變化對模態頻率的影響時，便需要進行預應力模態分析。
Abaqus預應力模態求解
分析流程如下：第一步先進行非線性靜力學求解——第二步進行模態提取
需要注意的是第一步求解時必須打開幾何非線性，即NLGEOM = YES 否則第一步求解完成后剛度矩陣不會改變，模態頻率也就不會發生變化。第二步模態求解無需設置PERTURBATION（線性攝動）或幾何非線性，軟件默認在開啟幾何非線性的后續分析步中繼續保持。
另外，第一步非線性靜力求解的材料非線性，接觸等都會對結構的剛度矩陣產生影響，進而改變模態頻率。材料如果進入塑性，相應的切向模量會降低，進而導致結構剛度矩陣變小。
靜力分析下接觸狀態的改變也會對剛度矩陣產生影響。Abaqus在進行預應力模態分析時對接觸的處理如下：第一步進行非線性接觸分析，軟件會把第一步分析結果的接觸區域作為第二步模態分析的作用區域，而第一步分析結果的接觸面分開區域不予考慮。需要注意的是，在進行第二步模態分析時，接觸區域并不是簡單的直接轉變為Tie處理，而是通過附加接觸剛度來進行求解。
Abaqus重啟動設置
重啟動分析方式是一種很便捷的模式。比如，我們需要算在預應力狀態下的模態，振動，沖擊等等一系列工況，而如果不進行重啟動分析，則每個分析工況下都需要重新計算預應力工況，對于大模型，嚴重影響計算效率；而進行重啟動設置后，預應力工況只需計算一次。
展開 abaqus有限元分析過程 附ABAQUS有限元分析常見問題解答下載
11.單元分類
1.分類方法：實體單元；殼單元；梁單元；彈簧單元；剛體單元；桁架單元；集中質量單元
也可以分為：
a.一維、二維和三維單元
b.線形、二次和三次單元
c.協調單元和非協調單元
d.傳彎單元和非傳彎單元
e.結構單元和非結構單元
下載地址：ABAQUS有限元分析常見問題解答
疲勞分析|Abaqus Goodman方法案例操作 附ABAQUS疲勞分析簡介下載
Abaqus/View結果讀取
讀取分析歷程中的最大交變應力和最小交變應力云圖
新建場變量：Alternating Stress和Mean Stress
根據公式：
在Abaqus后處理新建場變量
輸出場變量值到Excel
針對新建場，輸出單元積分點對應的交變應力和平均應力，并輸出到Excel，與Goodman圖一并繪制。
上圖，
仿真所得單元積分點落到
曲線的上方或下方，
處于上方為疲勞壽命沒達到
臨
界曲值
10
E5
次。
下載地址：ABAQUS疲勞分析簡介
abaqus電池包擠壓分析（附模型及分析流程） ￥46
1 問題設定 新能源汽車電池包擠壓分析的目的是采用 FEA 方法檢驗電池包是否可以滿足國標對電 池包擠壓性能的要求，包括電池包在擠壓過程中的結構變形、應力以及整體剛度等指標。
本 案例是利用 Abaqus2017 來建模以及求解。 電池包構件 電池包擠壓幾何模型（上下灰色的平板為剛體擠壓板）
部件的網格類型
以下內容包含完整的詳細教程，附件為完整教程文檔和CAE模型文件.rar
