abaqus 分析結(jié)果的案例
基于ABAQUS分析結(jié)果的Isight試驗(yàn)數(shù)據(jù)擬合
圖19 Optimization History選項(xiàng)卡
30.優(yōu)化任務(wù)運(yùn)算完成后會(huì)給出所有優(yōu)化方案中的最優(yōu)結(jié)果，可以看出從第33次迭代開始衡量指標(biāo)趨于平穩(wěn)，但在第36次迭代中達(dá)到了最小值（自動(dòng)綠色底紋標(biāo)注）。最優(yōu)迭代結(jié)果也將會(huì)以最終迭代結(jié)果最為最后輸出。
圖20 History選項(xiàng)卡
31.初始模擬、最終模擬和試驗(yàn)數(shù)據(jù)對(duì)比見圖21.可以看出
圖21 初始、最終模擬結(jié)果與試驗(yàn)數(shù)據(jù)對(duì)比
四、結(jié)論
從擬合結(jié)果可以看出，Isight的數(shù)據(jù)匹配可以有效地?cái)M合材料參數(shù)，主要材料參數(shù)與試驗(yàn)誤差明顯減少，見表2，如彈性模量誤差減小3.75%，軸心受壓強(qiáng)度誤差減小16.67%，這兩個(gè)參數(shù)對(duì)于模態(tài)分析和結(jié)構(gòu)承載力分析影響顯著?？焖儆行У仄ヅ涞秸鎸?shí)材料參數(shù)可以提高建模效率和分析準(zhǔn)確性。
基于ABAQUS分析結(jié)果的Isight試驗(yàn)數(shù)據(jù)擬.ppt
模型文件.rar
五、設(shè)備情況及計(jì)算耗時(shí)
CPU：AMD Athlon(tm) II X4 640 Processor 3.0 Ghz
內(nèi)存：4GB
計(jì)算耗時(shí)：406s
展開 Matlab“稍作修改”Abaqus-odb結(jié)果！！！
Abaqus是一款強(qiáng)大的有限元分析軟件，隨著版本的更迭，產(chǎn)品逐漸包含了前處理模塊、求解器、后處理模塊，用戶上手難度越來越小。
Matlab同樣也是一款強(qiáng)大的商業(yè)數(shù)值計(jì)算軟件，其可視化能力強(qiáng)的一批，作為興趣研究，可多玩玩這些工具聯(lián)合使用的效果，這次將這兩款工具聯(lián)合使用，介紹Matlab如何讀取Abaqus的.inp、.odb文件？
故事的開頭
木木平時(shí)喜歡玩一些有限元的東西，有一次在模型分析完畢后，我突發(fā)奇想：能否自定義obd文件，使之顯示成我想要的樣子？
我懷揣著這個(gè)想法，嘗試修改.odb文件，可是當(dāng)我打開文件后，傻眼了...
某.odb文件
這看不懂啊，好像Abaqus在說：少年，別動(dòng)我的數(shù)據(jù)！
我偏要對(duì)其作出修改！
進(jìn)入正題
主要編制了兩個(gè)函數(shù)文件：loadinp用于讀取.inp文件，meshplot用于可視化繪圖。該程序是國(guó)外一個(gè)大佬編制的小工具，結(jié)合具體案例可靈活使用。
Abaqus案例模型介紹
邊界條件及單元類型可見下圖，藍(lán)色區(qū)域使用CPS4單元，黃色區(qū)域使用CPS3單元
邊界條件
Abaqus分析結(jié)果如下，我將兩個(gè)Part分開展示，因?yàn)楸敬翁峁┑腗atlab程序不能完成裝配功能，為了與Matlab結(jié)果對(duì)比分析，只能分Part顯示。
Abaqus分析結(jié)果
Matlab讀入結(jié)果文件
導(dǎo)入inp文件，對(duì)文件中的單元節(jié)點(diǎn)信息進(jìn)行繪制部件模型。
展開 python基于分析結(jié)果的分析運(yùn)算
數(shù)學(xué)運(yùn)算允許對(duì)Abaqus分析結(jié)果進(jìn)行線性疊加或其他推導(dǎo)運(yùn)算。規(guī)則：
1） 允許對(duì)張量分量或向量進(jìn)行數(shù)學(xué)運(yùn)算
2） 可以根據(jù)分量來計(jì)算不變量
3） 不支持FieldOutput對(duì)象與HistoryOutput對(duì)象、FieldValue對(duì)象和HistoryOutput對(duì)象之間進(jìn)行數(shù)學(xué)運(yùn)算
4） 不支持兩個(gè)向量對(duì)象、兩個(gè)張量對(duì)象之間的乘除法運(yùn)算
5） 數(shù)學(xué)運(yùn)算表達(dá)式中的數(shù)據(jù)必須相容：向量和張量不支持求和；三維表面張量surface tensor和三維平面張量planar tensor不支持求和；積分點(diǎn)的結(jié)果數(shù)據(jù)不能夠與單元節(jié)點(diǎn)結(jié)果數(shù)據(jù)求和。
6） 如果通過getSubset方法獲取數(shù)學(xué)運(yùn)算表達(dá)式中的場(chǎng)變量，則必須使用getSubset方法按照相同順序獲取每個(gè)場(chǎng)變量。
7） 三角函數(shù)中各個(gè)參數(shù)的單位是弧度
8） 如果建立了局部坐標(biāo)系，張量的數(shù)學(xué)運(yùn)算基于局部坐標(biāo)系，否則，則給予整體坐標(biāo)系。
9） 如果FieldValue對(duì)象中的數(shù)據(jù)類型均相同，則允許對(duì)模型中不同位置的FieldValue對(duì)象進(jìn)行數(shù)學(xué)運(yùn)算；如果FieldValue對(duì)象的位置不同，則計(jì)算得到的FieldValue對(duì)象值與位置無關(guān)。如果參加運(yùn)算的FieldValue對(duì)象采用不同的局部坐標(biāo)系，計(jì)算時(shí)就不再考慮局部坐標(biāo)系的影響，因此，計(jì)算得到的fieldValue對(duì)象沒有基于任何局部坐標(biāo)系。
10） 不支持共軛數(shù)據(jù)（例如復(fù)數(shù)分析結(jié)果的虛部）的數(shù)學(xué)運(yùn)算。
※一般情況下，場(chǎng)變量FieldOutput對(duì)象比場(chǎng)變量值FieldValue object對(duì)象進(jìn)行數(shù)學(xué)運(yùn)算的效率高。得到FieldOutput對(duì)象的步驟：
① 在輸出數(shù)據(jù)庫(kù)中創(chuàng)建FieldOutput對(duì)象
② 調(diào)用addData方法向FieldOutput對(duì)象中添加計(jì)算得到的場(chǎng)變量對(duì)象。
展開 利用odb結(jié)果 使用fesafe進(jìn)行疲勞分析的基本操作流程 ￥29
基于abaqus分析結(jié)果odb文件，采用截圖+文字說明的方式介紹fe-safe進(jìn)行疲勞計(jì)算的操作流程。
ABAQUS分析參數(shù)設(shè)置對(duì)混凝土結(jié)構(gòu)分析結(jié)果的影響
采用一維梁?jiǎn)卧M混凝土梁壓彎試驗(yàn)，因?yàn)槭芾セ炷敛牧鲜掷_裂容易產(chǎn)生模型節(jié)點(diǎn)分析受力不平衡導(dǎo)致分析終止。針對(duì)混凝土材料特性，如果在默認(rèn)設(shè)置的基礎(chǔ)上略微調(diào)整分析設(shè)置參數(shù)，可能會(huì)明顯提高計(jì)算收斂性，達(dá)到預(yù)期目標(biāo)。下圖是混凝土梁分析參數(shù)調(diào)整的情況，和對(duì)分析結(jié)果的影響，希望能拋磚引玉。
調(diào)整前
*Boundary,amplitude=Cyclic
Load, 1, 1, -2000.
*Output, field
*Node Output
U,
調(diào)整后
*Boundary,amplitude=Cyclic
Load, 1, 1, -2000.
*Controls, ANALYSIS=DISCONTINUOUS
*Controls, reset
*Controls, parameters=line search
5, , , , 0.15
*Controls, parameters=field, field=displacement
0.05, 0.05
*Output, field
*Node Output
U,
調(diào)整參數(shù)含義：
①最大線性搜索步數(shù)設(shè)為5（即使用擬牛頓法）；
②線性搜索修正系數(shù)設(shè)為0.15；
③不平衡力與當(dāng)前平衡力范數(shù)容許比調(diào)整為0.05；
④最大修正值與對(duì)應(yīng)的增量值的容許比值調(diào)整為0.05；
結(jié)論：調(diào)整分析設(shè)置參數(shù)后模擬的終止加載力明顯增大。調(diào)整之前提前結(jié)束加載是因?yàn)楣?jié)點(diǎn)不平衡力超出容許值引起的。但對(duì)于混凝土材料來講，應(yīng)變積累導(dǎo)致突然的允許應(yīng)力下降，極易引起節(jié)點(diǎn)不平衡力增加，導(dǎo)致分析進(jìn)程結(jié)束。如果在允許范圍內(nèi)提高節(jié)點(diǎn)不平衡力容許范圍，則可以明顯增加加載幅度，達(dá)到預(yù)定分析目標(biāo)。
展開 Abaqus沖壓-回彈過程仿真詳細(xì)教程，顯式分析到隱式分析的結(jié)果傳遞方法 ￥99.9
沖壓回彈分析會(huì)涉及顯式求解器到隱式求解器之間的結(jié)果傳遞設(shè)置，這樣能夠?qū)F(xiàn)實(shí)中的力學(xué)過程進(jìn)行拆解，利用適當(dāng)?shù)那蠼馄?em>分析計(jì)算其對(duì)應(yīng)擅長(zhǎng)處理的的過程（動(dòng)態(tài)過程、穩(wěn)定過程），從而使整個(gè)分析效率極大地提高。
圖1 沖壓示意圖（1/4模型）
如圖1所示，毛坯（藍(lán)色）位于夾具（綠色）和模具（黃色）之間，沖頭（紅色）以一定的速度沖擊毛坯，毛坯在壓力和模具約束作用下發(fā)生一定的變形（沖壓過程）；隨后沖頭與夾具向上運(yùn)動(dòng)，卸載后的毛坯回彈并保留一定的永久變形（回彈過程），產(chǎn)品沖壓成型過程結(jié)束。
圖2 材料加、卸載的力學(xué)過程
材料加、卸載的過程中產(chǎn)生了彈性變形和塑性變形，分析時(shí)，通過Abaqus/Explicit分析其沖壓過程，再將分析結(jié)果作為初始狀態(tài)繼承給Abaqus/Standard進(jìn)行回彈分析。由于對(duì)稱性，使用一個(gè)1/4模型解決這個(gè)問題，全部采用殼單元。
展開 abaqus的三維和軸對(duì)稱模型分析的結(jié)果差異
一直懷疑abaqus在用三維模型和軸對(duì)稱模型分析同樣的東西的結(jié)果，
在動(dòng)力分析時(shí)軸對(duì)稱結(jié)果非常不可靠，與現(xiàn)場(chǎng)實(shí)測(cè)相差10倍，
而三維比較接近現(xiàn)場(chǎng)實(shí)測(cè)結(jié)果。
為此建了個(gè)簡(jiǎn)單的模型，用abaqus6.12做的，inp也附上，
大家一起探討一下。
直徑2m、高0.5m的圓柱體，彈性材料，彈性模量35E6Pa，泊松比0.35，柱頂面作用一個(gè)圓形荷載，1E6Pa，計(jì)算柱頂面中心點(diǎn)的最大位移。
分別用三維模型和軸對(duì)稱模型來模擬，結(jié)果見下面兩個(gè)圖，三維的頂面中性點(diǎn)位移1.026E-2，軸對(duì)稱1.151E-2。
inp.zip
---------------------------------------------------------------------------------------------------------------------------------------------
技術(shù)鄰?fù)扑]：
Hypermesh精典問答 （經(jīng)典加精品）
技術(shù)鄰學(xué)院:abaqus復(fù)合材料與cohesive教學(xué)視頻發(fā)布
Abaqus 二維hashin失效模型案例（附inp）
展開 Abaqus諧響應(yīng)分析結(jié)果的輸出
對(duì)于諧響應(yīng)分析，分析結(jié)果需要查看結(jié)構(gòu)在整個(gè)掃頻區(qū)間內(nèi)的響應(yīng)情況，例如：位移、應(yīng)力、應(yīng)變等。因此需要繪制整個(gè)頻率區(qū)間結(jié)構(gòu)的響應(yīng)曲線，方法如下。
進(jìn)入后處理模塊，選擇Creat XY date 如下圖所示：
如上圖2-4所示在Creat XY Date可以指定曲線縱坐標(biāo)的變量，可以是加位移、應(yīng)力、加速度、速度等，還可以定義繪制指定的分析步或頻率的結(jié)果，最后選擇節(jié)點(diǎn)，點(diǎn)擊Polt即可繪制該點(diǎn)相應(yīng)變量的曲線，如圖2-5所示。在利用Abaqus進(jìn)行諧響應(yīng)分析時(shí)，需要注意以下幾點(diǎn)：
①對(duì)于需要利用模態(tài)進(jìn)行多次諧響應(yīng)分析的問題，可以使用重啟動(dòng)，只需進(jìn)行一次模態(tài)分析，其它分析可在模態(tài)分析的基礎(chǔ)上進(jìn)行，可大大的減少計(jì)算成本。
②在利用模態(tài)疊加法進(jìn)行諧響應(yīng)分析時(shí)，模態(tài)頻率一定要將諧載荷的頻率段包含在內(nèi)，即，要包含掃頻區(qū)間。與此同時(shí)，模態(tài)三個(gè)方向的有效質(zhì)量要達(dá)到90%。
③在利用模態(tài)疊加法進(jìn)行諧響應(yīng)分析時(shí)，若結(jié)構(gòu)的模態(tài)頻率間隔較小，在諧響應(yīng)分析步中，在滿足計(jì)算要求的情況下，盡可能較小的定義特征頻率間的頻率點(diǎn)數(shù)與偏置。
展開 Abaqus動(dòng)態(tài)分析中，如何快速查看整個(gè)響應(yīng)過程中場(chǎng)輸出結(jié)果的最值 ￥9.9
<p>需求：動(dòng)態(tài)分析（基于模態(tài)的瞬態(tài)動(dòng)態(tài)響應(yīng)分析、顯示動(dòng)態(tài)分析等）中結(jié)果的響應(yīng)也是一個(gè)動(dòng)態(tài)的過程，不確定哪個(gè)時(shí)刻的結(jié)果是最大值或者最小值，或者說想知道整個(gè)響應(yīng)過程中的最大值、最小值是多少。結(jié)果輸出中是不會(huì)直接輸出的，只能看到每幀場(chǎng)輸出中的最值，又不可能自己逐幀場(chǎng)輸出結(jié)果里去看，然后找到所有幀中的最值，那么Abaqus軟件內(nèi)如何實(shí)現(xiàn)呢？</p><p><br></p><p><span style="background-color: rgb(255, 255, 255); color: rgb(25, 27, 31);">原創(chuàng)聲明：未經(jīng)本人同意，禁止抄襲、二次創(chuàng)作及轉(zhuǎn)載！</span></p>
展開 基于Abaqus的螺栓預(yù)緊力仿真分析
分析背景
Case1：鑄鋁卷邊部分分析
Case1邊界條件：
1、20KN針對(duì)兩個(gè)螺栓的預(yù)緊力（下圖藍(lán)色箭頭）
2、100N殼體內(nèi)部聚氨酯內(nèi)置件對(duì)鋁蓋的壓力（下圖紅色箭頭）
3、14KN來自減震器活塞桿對(duì)蓋子的載荷（下圖深紅色箭頭）
Case2：鋁殼底部受力分析
Case2邊界條件：
1、13.5KN來自減震器活塞桿對(duì)鋁殼底部的載荷（下圖黃色箭頭）
2、56.2KN來自jounce bumper針對(duì)鋁殼底部的載荷（下圖紅色箭頭）
幾何模型和FEA模型
分析結(jié)果
Case1：鑄鋁卷邊部分分析
鋁殼卷邊位置最大應(yīng)力143.3MPa，超過了材料的屈服強(qiáng)度（130MPa），在上圖紅色區(qū)域會(huì)有少許塑性變形。
鋁蓋上最大應(yīng)力199.6MPa，也超過了材料的屈服強(qiáng)度，在局部位置會(huì)有塑性變形。
Abaqus的分析結(jié)果
鋁殼卷邊位置，由于abaqus采用C3D10M單元，部分位置有應(yīng)力奇異現(xiàn)象，應(yīng)力值偏大，但應(yīng)力的整體分布與前者基本一致，從應(yīng)變圖可以看出來。
Case2：鋁殼底部受力分析
鋁殼應(yīng)力云圖
鋁殼根部有較大應(yīng)力，同時(shí)會(huì)產(chǎn)生較小的塑性變形。由于ALSi12本身延展性較差，局部位置會(huì)出現(xiàn)失效。
鋁殼應(yīng)力云圖（abaqus結(jié)果）
應(yīng)力分布基本一致，產(chǎn)生塑性應(yīng)變的區(qū)域也差不多。
結(jié)論
工況一條件下，鋁殼卷邊位置會(huì)有很小的區(qū)域有塑性變形，不影響整體使用功能。
工況二條件下，鋁殼大部分位置應(yīng)力都超過了屈服強(qiáng)度，而且局部位置應(yīng)力超過抗拉強(qiáng)度，有失效的風(fēng)險(xiǎn)。
在實(shí)際工程應(yīng)用中，常常使用力矩扳手或超聲波測(cè)力法進(jìn)行測(cè)量螺栓預(yù)緊力，但是這兩種方法局限性很大，測(cè)量結(jié)果并不太可靠。
展開 Matlab“稍作修改”Abaqus-odb結(jié)果 | 提取三維裝配體Mises應(yīng)力
Abaqus分析結(jié)果如下（取step-1-increment28）：
Matlab讀入網(wǎng)格文件
導(dǎo)入inp文件，對(duì)文件中的單元節(jié)點(diǎn)信息進(jìn)行繪制部件模型。
Abaqus與Dyna電池包沖擊分析結(jié)果對(duì)比 ￥20
對(duì)比二：不同求解器的結(jié)果對(duì)比
位移結(jié)果動(dòng)畫
計(jì)算結(jié)果
能量曲線對(duì)比
Dyna能量曲線
Abaqus能量曲線
為了更好的解讀計(jì)算結(jié)果，我們先查看鈑金件的材料曲線
塑性段材料曲線
由計(jì)算結(jié)果可知：
1. Dyna中減縮積分單元類型2和全積分單元16結(jié)果相差較大。16號(hào)單元加上沙漏類型8可以很好的處理單元扭曲變形的問題；2號(hào)單元由于采用面內(nèi)單點(diǎn)積分，計(jì)算速度很快，對(duì)于大變形問題的計(jì)算最穩(wěn)定，但對(duì)于單元扭曲變形非常敏感，計(jì)算精度相對(duì)較差。所以推薦使用單元類型16，這也是各主機(jī)廠通常使用的單元類型。
2. Abaqus減縮積分S4R和全積分S4位移結(jié)果基本完全一致，單元應(yīng)力應(yīng)變結(jié)果相近（減縮積分應(yīng)力應(yīng)變偏?。e分S4外推到節(jié)點(diǎn)的應(yīng)力應(yīng)變值差異很大，而且外推后應(yīng)力應(yīng)變不再滿足材料本構(gòu)關(guān)系（應(yīng)力340MPa對(duì)應(yīng)的塑性應(yīng)變?yōu)?%）。所以推薦使用整體S4R+局部關(guān)心區(qū)域S4的方式進(jìn)行建模，這樣可以兼顧計(jì)算速度和計(jì)算精度，當(dāng)然結(jié)果處理不要用外推到節(jié)點(diǎn)的應(yīng)力應(yīng)變。
3. Dyna中16號(hào)單元與Abaqus中S4的應(yīng)力應(yīng)變結(jié)果相近，誤差最大的反而是位移。能量曲線趨勢(shì)一致，但數(shù)值上相差較大，可能是由于各自的接觸算法和單元理論假定的不同導(dǎo)致。
歡迎關(guān)注微信公眾號(hào)：共享CAE
付費(fèi)內(nèi)容為相應(yīng)計(jì)算文件，包含完整的電池包模型，Dyna和abaqus帶損傷的材料韌性斷裂失效設(shè)置，abaqus explicit通用接觸防止接觸厚度自動(dòng)減小的設(shè)置，及其他常規(guī)的電池包沖擊分析設(shè)置，供大家學(xué)習(xí)。
展開 Abaqus熱分析關(guān)于結(jié)果中溫度的后處理二次開發(fā) ￥5
在abaqus中，進(jìn)行熱分析的時(shí)候，經(jīng)常要對(duì)得到的溫度結(jié)果進(jìn)行處理。如下圖所示：
我們要對(duì)這個(gè)面上的溫度數(shù)據(jù)，進(jìn)行處理，比如：
（1）我要找出這個(gè)面上溫度符合某一條件的區(qū)域。
（2）我要找出某一時(shí)刻，這一面上所有節(jié)點(diǎn)溫度的平均值。
要解決以上兩個(gè)問題，最好的辦法就是用python對(duì)ABAQUS的ODB文件進(jìn)行訪問處理。用python來處理ODB文件中的歷史輸出數(shù)據(jù)。我寫了一個(gè)函數(shù)，用以解決以上兩個(gè)方面的問題。大家可以從代碼得到啟發(fā)，思考如何對(duì)abaqus的歷史輸出數(shù)據(jù)進(jìn)行訪問。有相關(guān)ABAQUS二次開發(fā)方面的問題，我們可以探討探討。QQ:1917665644
具體代碼如下：
展開 ABAQUS動(dòng)態(tài)分析中的能量平衡、沙漏及結(jié)果評(píng)估
ABAQUS動(dòng)態(tài)分析中的能量平衡、沙漏及結(jié)果評(píng)估
考慮壩體-庫(kù)水相互作用的重力壩模態(tài)分析--對(duì)比分析ANSYS和ABAQUS重力壩流固耦合模態(tài)結(jié)果
ANSYS壩體-庫(kù)水有限元模型
ABAQUS壩體-庫(kù)水幾何模型及約束條件
（1）首先對(duì)比分析，ANSYS與ABAQUS混凝土重力壩空庫(kù)模態(tài)分析結(jié)果
（2）壩體-庫(kù)水流固耦合模型，模態(tài)分析結(jié)果對(duì)比
ANSYS與ABAQUS計(jì)算前12階模態(tài)對(duì)比分析結(jié)果
階數(shù)
ABA滿庫(kù)自振頻率/Hz
ANS滿庫(kù)自振頻率/Hz
ANSYS-ABAQUS誤差百分比
1
5.5487
5.5717
0.004145115
2
6.7567
6.7702
0.001998017
3
9.5858
9.6043
0.001929938
4
13.922
14.013
0.006536417
5
16.358
16.42
0.003790194
6
17.76
17.786
0.001463964
7
19.648
19.741
0.004733306
8
展開 