abaqus自由變形的案例
批量提取Abaqus的節(jié)點(diǎn)坐標(biāo)(初始坐標(biāo)、指定Step下的變形量、變形后節(jié)點(diǎn)坐標(biāo)) ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發(fā)來批量提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中指定Step下的Set節(jié)點(diǎn)集變形量。通過詳細(xì)的步驟說明、代碼示例和圖片展示，您將學(xué)會(huì)如何使用該腳本，自動(dòng)化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，U1，U2）。</p><p>如果還需要按Increment提取每個(gè)增量下的變形后的節(jié)點(diǎn)坐標(biāo)的話，在提取變形量的基礎(chǔ)上，與初始坐標(biāo)進(jìn)行簡單的計(jì)算就可以求得坐標(biāo)。 （備注：該代碼只提取了x，y方向的變形量）</p><h2>1. 問題描述</h2><p>在工程仿真和分析領(lǐng)域，提取ABAQUS輸出數(shù)據(jù)庫（ODB）文件中的節(jié)點(diǎn)集變形量是一項(xiàng)常見任務(wù)。然而，手動(dòng)提取這些數(shù)據(jù)是一項(xiàng)繁瑣且容易出錯(cuò)的工作。因此，需要一種自動(dòng)化的方法來批量提取指定步驟下按節(jié)點(diǎn)集組織的變形量數(shù)據(jù)。</p><h2>2. 實(shí)例展示</h2><p>假設(shè)我們有一個(gè)名為`example.odb`的ODB文件，其中包含名為`Step-x`的步驟和名為`Set-x`的節(jié)點(diǎn)集。運(yùn)行以上代碼后，腳本會(huì)自動(dòng)將該步驟下節(jié)點(diǎn)集的變形量提取出來，并保存為`NodalDisplacement.csv`文件。
展開 abaqus的Fortran子程序中使用自由格式
在上一篇文章Fortran語言的自由格式與固定格式中，曾經(jīng)討論過Abaqus的Fortran子程序在默認(rèn)狀態(tài)下只能使用固定格式，而實(shí)際上自由格式更為靈活好用，其一行不受72個(gè)字符的限制，并且可以將多個(gè)語句寫在同一行，同時(shí)續(xù)行符相對(duì)于固定格式也更為好用。本文主要描述下如何實(shí)現(xiàn)在abaqus中采用Fortran自由格式編程。
以幫助文檔中的經(jīng)典塑性UMAT子程序?yàn)槔? （1）打開Abaqus documentation，點(diǎn)擊Abaqus Verification Guide，并進(jìn)入U(xiǎn)MAT AND UHYPER下的input files；
（2）下載inp文件及f子程序，用文本編輯器打開umatmst3.f，可以看出其采用的是固定格式編寫；
（3）修改abaqus的運(yùn)行環(huán)境，加入支持自由格式的設(shè)置：找到abaqus環(huán)境的配置文件，對(duì)于abaqus2020和abaqus2021版本，該文件為以下路徑的win86_64.env文件：
對(duì)于比較早的abaqus版本如abaqus6.14，該文件可能叫abaqus_v6.env?？傊ㄟ^搜索等方式找到該文件并打開即可。
（4）在.env文件中加入使用自由格式的命令：
（5）將子程序修改為自由格式，需要修改的內(nèi)容主要包括 續(xù)航符，注釋等。
展開 彈簧下樓梯Abaqus仿真，奇怪的自由落體行為之謎
可能大家都是小時(shí)候玩的了，今天咱們就用Abaqus來玩玩slinky下樓梯和“詭異”的自由落體，回味一下童年的樂趣。
案例1. slinky下樓梯仿真
它的軸向剛度 =0.002829095N/mm
一般的彈簧軸向剛度是多大呢？1~200N/mm !!! 相比之下，slinky幾乎沒有軸向剛度可言！
正因?yàn)檫@樣，slinky第一階梯走完之后，才有可能由收尾慣性段抵抗很小的軸向剛度，甩向下一階梯，并重復(fù)上一狀態(tài)。這是一個(gè)常規(guī)的動(dòng)力學(xué)分析，材料屬性是普通的塑料線彈性本構(gòu)，模型的重點(diǎn)還是在邊界條件的處理上，建模時(shí)以下兩點(diǎn)需要注意：
殼單元可以節(jié)省計(jì)算量，最開始我嘗試的思路是采用更簡單的beam，但是接觸處理起來不太方便，通過嘗試發(fā)現(xiàn)使用殼時(shí)彈簧會(huì)站的比較穩(wěn)；
設(shè)置合適的邊界條件，保證走下第一個(gè)臺(tái)階時(shí)橫向位移不要超過臺(tái)階寬度，不然slinky很容易跑偏或整體滑下樓梯。
案例2. “詭異”的自由落體仿真
像上圖那樣丟過slinky的人可能會(huì)發(fā)現(xiàn)，它的底端，會(huì)神奇地表現(xiàn)出“反重力”的特點(diǎn)，是不是和你頭腦中的自由落體的感覺不太一樣？
為什么會(huì)這樣呢？
讓我們用Abaqus來一探究竟！
首先，要明白一個(gè)概念：初始應(yīng)力狀態(tài)。一般情況下我們丟的東西初始狀態(tài)都很放松，所以都表現(xiàn)正常。這個(gè)slinky呢，在吊起來被釋放之前，由于重力，體內(nèi)儲(chǔ)存了大量的應(yīng)變能，也就是通常說的彈性勢能，那么，約束解除之后，應(yīng)變能去哪了？
先看看力學(xué)過程，這個(gè)仿真可以分兩個(gè)過程，第一個(gè)過程是求解彈簧在重力作用下變形的靜力學(xué)過程；第二個(gè)過程是求解變形后的彈簧在重力作用下自由落體的動(dòng)力學(xué)過程。
展開 兼職/遠(yuǎn)程】高薪招募ABAQUS仿真工程師 | 時(shí)間自由 | 長期穩(wěn)定合作
?? 你的核心職責(zé)
項(xiàng)目承接：承接平臺(tái)分發(fā)的各類ABAQUS仿真需求，涵蓋結(jié)構(gòu)靜力學(xué)/動(dòng)力學(xué)、非線性分析（接觸/材料非線性）、熱-力耦合、顯式動(dòng)力學(xué)（Explicit） 等方向。
技術(shù)支持：根據(jù)客戶提供的模型或圖紙，獨(dú)立完成幾何清理、網(wǎng)格劃分、求解設(shè)置、結(jié)果后處理及仿真報(bào)告撰寫。
專業(yè)背景：
本科及以上學(xué)歷（優(yōu)秀的在讀本碩博士亦可），力學(xué)、機(jī)械工程、車輛工程、材料科學(xué)與工程等相關(guān)專業(yè)。
具備扎實(shí)的有限元理論功底，熟悉材料力學(xué)、彈塑性力學(xué)、振動(dòng)理論等基礎(chǔ)學(xué)科 -
軟件技能：
精通 Abaqus/CAE 軟件，熟練掌握Standard與Explicit求解器。
熟練使用HyperMesh或ANSA等專業(yè)前處理軟件者優(yōu)先。
熟悉Python腳本進(jìn)行Abaqus二次開發(fā)（參數(shù)化建模、結(jié)果自動(dòng)處理）是加分項(xiàng)。
經(jīng)驗(yàn)與時(shí)間：
有完整的工程項(xiàng)目仿真經(jīng)驗(yàn)（學(xué)術(shù)論文案例不算，需解決實(shí)際工程問題）-
有充裕的業(yè)余時(shí)間（如晚上、周末），能夠靈活響應(yīng)項(xiàng)目需求，具備較強(qiáng)的責(zé)任心和契約精神。
軟性素質(zhì)：
具備良好的溝通能力和服務(wù)意識(shí)，能夠清晰地向非技術(shù)背景客戶解釋仿真結(jié)果的含義。
有技術(shù)鄰平臺(tái)接單經(jīng)驗(yàn)或已有技術(shù)鄰認(rèn)證者優(yōu)先 -
工作說明:我們每天都有在線發(fā)布任務(wù)，待遇優(yōu)厚、時(shí)間自由、地點(diǎn)自由，線上辦公 ，不耽誤您個(gè)人事宜。薪酬從幾百幾千不等，多勞多得，同時(shí)積累項(xiàng)目經(jīng)驗(yàn)。
添加備注來源，以便我們后續(xù)安排
符合以上信息并且感興趣的朋友可以發(fā)送郵件到: 19283933696@163.com網(wǎng)易郵箱
招聘信息長期有效，歡迎咨詢
展開 
abaqus的Fortran子程序中使用自由格式與固定格式
帖子Fortran語言的自由格式與固定格式指出Abaqus的Fortran子程序在默認(rèn)狀態(tài)下只能使用固定格式，而實(shí)際上自由格式更為靈活好用，其一行不受72個(gè)字符的限制，并且可以將多個(gè)語句寫在同一行，同時(shí)續(xù)行符相對(duì)于固定格式也更為好用。同時(shí)自由格式能夠采用很多現(xiàn)代Fortran的語法來編寫程序，簡化程序編寫，有相當(dāng)大的優(yōu)勢。另外一方面，目前現(xiàn)存很多代碼，課題組祖?zhèn)鞒绦虼蠖嗍枪潭ǜ袷骄帉?，這些代碼量多，質(zhì)量久經(jīng)考驗(yàn)，我們在新編代碼中又希望能使用他們，那么有沒有方法在一個(gè)文件中同時(shí)使用兩種風(fēng)格的代碼呢。答案是可以的，本文主要描述下如何實(shí)現(xiàn)在abaqus中采用Fortran自由格式編程以及自由格式和固定格式混編。
以for文件默認(rèn)是固定格式，可以通過!DIR$ FREEFORM 和 !DIR$ NOFREEFORM 分別控制代碼個(gè)編譯格式。以下兩個(gè)例子詳細(xì)描述。
1. 文檔中自由格式和固定格式混合。
計(jì)算結(jié)果：
2. 一個(gè)子程序中自由格式和固定格式混寫。
這種情況基本很少遇到，這是個(gè)騷操作，我也不過多介紹，直接搬運(yùn)IVF官方例子：
總結(jié)：Abaqus的for文件可以采用自由格式編寫，或者自由格式和固定格式混編。只需要在文件中加上!DIR$ FREEFORM 和 !DIR$ NOFREEFORM 進(jìn)行格式控制就行，不用修改環(huán)境文件和使用格式轉(zhuǎn)換鵝毛筆，這種方法簡單實(shí)用，希望對(duì)大家有所幫助。
展開 關(guān)于Abaqus UEL中RHS數(shù)組長度大于單元總自由度數(shù)的一些解釋
在Abaqus中，用戶自定義元素子程序（UEL）的開發(fā)需要遵循一些特定的規(guī)則和約定。其中一個(gè)關(guān)鍵約定是關(guān)于子程序中的RHS（右手邊）向量的維數(shù)，我最近在嘗試用UEL做一些二次開發(fā)，也發(fā)現(xiàn)了RHS向量的維數(shù)比單元的總自由度數(shù)多了4個(gè)這一現(xiàn)象，結(jié)合在站內(nèi)一些同行的猜測，我認(rèn)為這可能是由于Abaqus的內(nèi)部工作方式所導(dǎo)致的。
首先，有同行懷疑是因?yàn)橛昧怂墓?jié)點(diǎn)單元，所以多了四個(gè)，我開發(fā)的單元是12個(gè)節(jié)點(diǎn)的，依然多了四個(gè)維度，因此排除是單元內(nèi)節(jié)點(diǎn)個(gè)數(shù)導(dǎo)致的。
在Abaqus中，RHS向量的維數(shù)實(shí)際上包括了除了單元的位移自由度外的其他項(xiàng)，這些項(xiàng)用于處理多種情況，例如：
體積力和表面力的計(jì)算：RHS向量可能包括用于計(jì)算體積力和表面力的額外自由度。這些自由度用于存儲(chǔ)單元內(nèi)的體積力和表面力的貢獻(xiàn)。
約束和邊界條件：Abaqus可能需要額外的自由度來處理約束條件和邊界條件，以確保數(shù)值穩(wěn)定性和正確的求解結(jié)果。
內(nèi)部狀態(tài)變量：某些材料模型和非線性分析可能需要存儲(chǔ)和更新一些內(nèi)部狀態(tài)變量，這些變量也可以占用RHS向量中的額外位置。
因此，RHS向量的維數(shù)不僅僅包括單元的位移自由度，還包括其他與分析和模型特性相關(guān)的項(xiàng)。這是Abaqus設(shè)計(jì)的一部分，旨在確保通用性和可擴(kuò)展性，以處理各種復(fù)雜的問題。其中我認(rèn)為可擴(kuò)展性是一項(xiàng)比較重要的應(yīng)用，我目前所做的工作可能會(huì)利用到這一點(diǎn)，等有結(jié)果了會(huì)繼續(xù)更新。
展開 Abaqus鼓風(fēng)機(jī)葉輪模態(tài)(自由+重力+離心力)及諧響應(yīng)分析
[圖片]
萬類霜天競自由——Abaqus任意移動(dòng)熱源插件 焊接 ￥600
通過本款FreeWeld插件可以自由定義焊接路徑，傻瓜式操作，只需在窗口界面選取幾何邊特征作為移動(dòng)路徑，就能自動(dòng)生成相應(yīng)的DFLUX子程序。程序中的熱源采用高斯面熱源，參數(shù)Rh為高斯熱源的特征半徑。
經(jīng)過簡單修改可以實(shí)現(xiàn)移動(dòng)載荷DLOAD路徑的自定義。
（本插件支持單熱源生成，如需多熱源插件請見：https://www.yqgqt.org.cn/content/post/1283087）
實(shí)例1效果及插件操作過程：
實(shí)例2效果及插件操作過程：
對(duì)于空間曲線路徑同樣支持：
tips：
1. 運(yùn)行程序生成子程序文件之后后，為了獲得更好的網(wǎng)格質(zhì)量，可以刪除移動(dòng)路徑的邊線特征，不會(huì)影響移動(dòng)熱源程序運(yùn)行；
2. 移動(dòng)熱源使用注意事項(xiàng)：
① 應(yīng)選擇溫度-位移耦合分析步或傳熱分析步；
② 在需要加載移動(dòng)熱源的面上施加自定義表面熱流載荷，如下圖：
③ 材料屬性應(yīng)涵蓋密度、比熱容、熱導(dǎo)率、彈性模量、泊松比和熱膨脹系數(shù)；
展開 單自由度彈簧振子abaqus的實(shí)例設(shè)置，初始條件為初速度。
彈簧振子的振動(dòng)是我們在學(xué)習(xí)工作振動(dòng)力學(xué)和可靠性理論中最簡單，最基礎(chǔ)的單元，掌握其理論基礎(chǔ)，并用abaqus仿真。
現(xiàn)將實(shí)例概述如下：質(zhì)量1噸的物體用彈簧和固定的甲板相連，彈簧的剛度k=2000N/m，阻尼系數(shù)c=17.89N·s/m,求最大位移。
1振動(dòng)力學(xué)的原理：
2利用MATLAB4級(jí)龍格庫塔得出以下結(jié)論：
最大位移位x=338mm
3利用ANSYS得出的結(jié)論
最大位移x=338mm
4用abaqus仿真設(shè)置和結(jié)論
這個(gè)例仿真對(duì)于初學(xué)者有三 個(gè)難點(diǎn)，第一個(gè)是彈簧的設(shè)置，第二個(gè)是初始速度的施加。第三個(gè)是分析步的設(shè)置。
a.裝備。將模型簡化成兩個(gè)點(diǎn)，可以在裝配中直接做出來，兩點(diǎn)的距離不影響結(jié)果，為了觀察方便，距離要適當(dāng)?shù)脑龃蟆?b.分析步。分析步要設(shè)成兩步，第一步靜力學(xué)分析，第二步隱式動(dòng)力學(xué)分析。也可以適當(dāng)?shù)恼{(diào)節(jié)增量步的大小，也可以不改變。
c.相互作用。設(shè)置彈簧的剛度、阻尼和慣性。其中彈簧的設(shè)置有兩種這里只介紹一種特殊設(shè)置（另一種也不是很復(fù)雜）。
d.載荷。這里要將左側(cè)的甲板固定住，限制他的六個(gè)自由度。并且施加初速度，這里要注意初學(xué)者容易將初速度在邊界條件中施加，這并不是正確的。要在預(yù)定場當(dāng)中施加。
c.網(wǎng)格。因?yàn)槟Ｐ鸵呀?jīng)簡化成了兩個(gè)點(diǎn)，因此不必要在進(jìn)行網(wǎng)格劃分，直接提交作業(yè)。
d.可視化后處理
彈簧振子的位移曲線如圖，最大位移為338mm
5結(jié)論
最后，MATLAB數(shù)值仿真，ansys與abaqus結(jié)論相同。
本實(shí)例主要針對(duì)abaqus的初學(xué)者的彈簧振子的相關(guān)問題提供思路，由于水平有限，歡迎批評(píng)指正。Q:1035863272
展開 Abaqus+PyQt+Python平面變形歐拉角計(jì)算
下面以簡單例子介紹平面變形、指向歐拉角的計(jì)算，包括絕對(duì)歐拉角、相對(duì)歐拉角。
1 簡化模型
下面的六面體為表面殼模型，下面由三段梁支持，三段梁分別沿X、Y、Z軸向。六個(gè)面的厚度不同，在上側(cè)3個(gè)面施加不同的壓力，如下左圖所示。位移云圖如下右圖所示。
2 計(jì)算要求
計(jì)算六面體上面3個(gè)面的變形歐拉角，包括3個(gè)面的絕對(duì)歐拉角，平面2、3相對(duì)與平面1的相對(duì)歐拉角。平面1、2、3如下圖所示。
3 數(shù)據(jù)處理
使用平面節(jié)點(diǎn)坐標(biāo)、位移數(shù)據(jù)計(jì)算平面變形歐拉角?？梢允褂肞ython腳本輸出平面節(jié)點(diǎn)編號(hào)、節(jié)點(diǎn)坐標(biāo)（X、Y、Z）、節(jié)點(diǎn)位移（U1、U2、U3），如下圖所示。下圖為平面1的10個(gè)工況的數(shù)據(jù)文件，打開的文本文件中7列數(shù)據(jù)為節(jié)點(diǎn)編號(hào)、坐標(biāo)、位移。
三個(gè)平面10個(gè)工況的節(jié)點(diǎn)數(shù)據(jù)文件如下圖所示。每個(gè)文件中包含一個(gè)工況一個(gè)平面的節(jié)點(diǎn)編號(hào)、坐標(biāo)、位移數(shù)據(jù)。
4 絕對(duì)歐拉角計(jì)算
使用PyQt+Python開發(fā)了一個(gè)簡單的小軟件，計(jì)算絕對(duì)歐拉角、相對(duì)歐拉角。
首先計(jì)算各平面的絕對(duì)歐拉角。
計(jì)算平面1的10個(gè)工況的絕對(duì)歐拉角。
平面1變形的絕對(duì)歐拉角計(jì)算結(jié)果如下圖所示。
伴隨絕對(duì)歐拉角計(jì)算結(jié)果，軟件同時(shí)寫出了平面變形前后的坐標(biāo)系數(shù)據(jù)，如下圖。每行18個(gè)數(shù)據(jù)，每3個(gè)數(shù)據(jù)為一個(gè)坐標(biāo)軸向量，變形前后2個(gè)坐標(biāo)系，6個(gè)坐標(biāo)軸，18個(gè)數(shù)據(jù)。
5 相對(duì)歐拉角計(jì)算
利用計(jì)算絕對(duì)歐拉角時(shí)得到的坐標(biāo)系文件，計(jì)算平面變形相對(duì)歐拉角，如下圖所示，計(jì)算平面2相當(dāng)于平面1、平面3相對(duì)與平面1的相對(duì)歐拉角。
計(jì)算結(jié)果如下圖所示。
6 小結(jié)
上述軟件用的算法申請了發(fā)明專利，軟件申請了軟著。
展開 Abaqus模擬橡膠大變形
Abaqus為用戶提供了多種本構(gòu)關(guān)系來模擬超彈性材料，這種材料具有高度非線性，當(dāng)Abaqus進(jìn)行模擬時(shí)假設(shè)這種材料是具有彈性、各向同性，并且同時(shí)考慮幾何非線性效應(yīng)。與材料的剪切柔度相比，對(duì)于大多數(shù)類似橡膠的固體材料，其可壓縮性非常小，當(dāng)分析對(duì)象為平面應(yīng)力問題、殼、薄膜、梁、桁架、或者鋼筋等，這個(gè)問題不值得關(guān)注。但是對(duì)于固體、平面應(yīng)變或者軸對(duì)稱問題卻不能忽略。對(duì)此，Abaqus/Standard提供了雜交單元來模擬超彈性材料中完全的不可壓縮行為。
橡膠材料力學(xué)性能的描述方法主要為兩類：一類是認(rèn)為橡膠為連續(xù)介質(zhì)的現(xiàn)象學(xué)描述；另一類是基于熱力學(xué)統(tǒng)計(jì)的方法?；谶B續(xù)介質(zhì)力學(xué)的本構(gòu)模型主要有Polynomial、Reduce Polynomial、Ogden模型等，其中Mooney-Rivlin模型是 Polynomial的特殊形式，Neo-Hookean 模型是Reduce Polynomial的特殊形式。基于熱力學(xué)統(tǒng)計(jì)主要有Arruda-Boyce和Van der Waals等本構(gòu)模型。本文利用Abaqus模擬大變形的橡膠，具體步驟如下。
1、在Abaqus/CAE Sketch模塊中作出模型草圖，如圖1所示，然后在Part模塊中分別建立Push、Rubber、Base三個(gè)部件。其中Push為解析剛體，Base為離散剛體。
圖1 草圖
2、在Property模塊中定義橡膠的屬性，采用Mooney-Rivlin模型，參數(shù)如圖2所示，然后賦給Rubber部件。
圖2 橡膠參數(shù)設(shè)置
3、裝配，定義分析步，采用默認(rèn)的場輸出和歷史輸出。為了保證剛開始能夠較容易收斂，設(shè)置分析步初始增量步為0.01，打開幾何非線性。
圖3 分析步定義
4、定義接觸對(duì)：Push下表面和橡膠表面，Base上表面和橡膠表面。
展開 
ABAQUS 單向拉伸大變形模擬
靜態(tài)模擬一種軟材料POE的單向拉伸，拉伸應(yīng)變希望到300%，但是總是在100%就失敗了。不知道哪里出了問題，有沒有高手幫幫忙。
HyperMesh導(dǎo)入Abaqus的變形后模型
HyperMesh導(dǎo)入Abaqus的變形后模型
經(jīng)過Abaqus的操作可以將Abaqus計(jì)算后的變形單元和節(jié)點(diǎn)輸出為inp文件，這個(gè)inp文件可以通過HyperMesh讀取。
打開HyperMesh后設(shè)置user files為Abaqus，如圖1所示。
圖1
之后通過import Solve Deck導(dǎo)入Abaqus生成的inp文件，如圖2所示。
圖2
之后點(diǎn)擊import即可在窗口顯示導(dǎo)入的變形模型，如圖3所示。
圖3
在此基礎(chǔ)上，可以對(duì)網(wǎng)格進(jìn)行編輯重劃分，修改結(jié)束后同樣可以導(dǎo)出修改后的模型inp文件，如圖4所示。
圖4
此時(shí)的inp文件記錄了在HyperMesh中修改后的模型。
再之后，可以通過Abaqus導(dǎo)入Hm生成的inp文件，如圖5所示。
圖5
導(dǎo)入后的模型在Abaqus中的顯式如圖6所示。
圖6
在HyperMesh編輯三維單元網(wǎng)格不是那么方便，建議的思路是先生成實(shí)體再進(jìn)行網(wǎng)格編輯。
Abaqus6.12以后的版本據(jù)說可以由孤立的網(wǎng)格生成實(shí)體模型，這個(gè)還是比較方便的，生成實(shí)體模型后在Abaqus里面就可以進(jìn)行網(wǎng)格重劃分。
展開 Abaqus薄板彎曲變形分析實(shí)
ABAQUS提供了業(yè)內(nèi)領(lǐng)先的接觸建模能力，接觸中各種表面間的各類摩擦性質(zhì)可以建立相應(yīng)的模型模擬，來符合不同接觸行為的要求。
本文采用Abaqus/Standard求解器，進(jìn)行薄板彎曲變形分析，用以簡單展示ABAQUS接觸建模及其分析功能。
1、 計(jì)算模型
如圖1所示，懸臂梁左端受剛性模具固定，右端受移動(dòng)模具下壓產(chǎn)生變形。
2、 有限元模型
建立有限元模型，創(chuàng)建穩(wěn)態(tài)分析步，分析薄板和剛性表面間的接觸，平板使用實(shí)體平面應(yīng)變單元CPE4I, 該單元沿板厚方向只需要一個(gè)單元即可以準(zhǔn)確模擬彎曲行為。剛性表面以解析剛性面模擬。
3、 接觸建立
ABAQUS中，接觸的一般需要三個(gè)步驟。
首先定義接觸表面。剛性表面一般作為接觸對(duì)的主面，本例中將剛性模具的面定義為主面，薄板面為從面。
進(jìn)而定義接觸對(duì)。選擇發(fā)生接觸的主從面定義為接觸對(duì)。
最后定義接觸屬性。包括接觸類型，以及摩擦系數(shù)等相關(guān)接觸參數(shù)。本例選擇無摩擦的光滑接觸屬性。
本案例共包括三個(gè)接觸對(duì)，分別為三個(gè)剛性模具與薄板之間的接觸。
完成接觸設(shè)定后，對(duì)模型設(shè)定相關(guān)邊界條件：上下模具完全固定，沖頭向下移動(dòng)60mm。薄板左端固定。
在此邊界條件下，沖頭向下移動(dòng)時(shí)，薄板上的三個(gè)接觸對(duì)發(fā)生作用，使得薄板右端發(fā)生彎曲。
4、 接觸輸出
接觸設(shè)定中，對(duì)于多有表面的接觸信息，可以設(shè)定接觸應(yīng)力、接觸位移等接觸輸出信息。
5、 分析結(jié)果
如圖所示，計(jì)算完成后薄板發(fā)生預(yù)想彎曲。案例設(shè)定了接觸應(yīng)力輸出，接觸應(yīng)力包括接觸壓力、摩擦剪切力的輸出，均可以在后處理中進(jìn)行相應(yīng)結(jié)果顯示。圖中所示云圖所示為接觸壓力云圖。
展開 ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬 ￥70
ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬
一、建模技術(shù)
地震工況下邊坡可能失穩(wěn)進(jìn)而出現(xiàn)滑坡現(xiàn)象，為避免模擬滑坡時(shí)網(wǎng)格產(chǎn)生的畸變問題，采用耦合歐拉拉格朗日法(CEL)進(jìn)行滑坡的大變形模擬；土體本構(gòu)采用摩爾庫倫模型；采用模型底部小范圍內(nèi)的周期性荷載模擬地震荷載。
二、模型及部分結(jié)果展示
圖1：藍(lán)色為邊坡；紅色為空氣層
圖2：網(wǎng)格的劃分
圖3：賦予模型初始應(yīng)力
圖4：土體達(dá)到地應(yīng)力平衡時(shí)的應(yīng)力分布
圖5：土體底部的地震荷載施加區(qū)域
圖6：所施加的周期性荷載（地震荷載）
圖7：邊坡因地震荷載產(chǎn)生的位移
圖8：地震波產(chǎn)生的區(qū)域
展開 