abaqus網(wǎng)格編程的案例
自編程實(shí)現(xiàn)不同單元類型的網(wǎng)格映射（包括像素網(wǎng)格）-原創(chuàng)帖
前面帖子講解的都是基于第三方軟件進(jìn)行像素網(wǎng)格或者界面自適應(yīng)網(wǎng)格的生成方法，在本帖簡單展示下通過自編程插件或程序?qū)崿F(xiàn)網(wǎng)格映射，這樣做的優(yōu)點(diǎn)是：不針對某類單元、不針對幾何形狀、不針對幾何空間、整體效率高、方式多樣、操作簡單、節(jié)約時間，缺點(diǎn)是：不能處理真實(shí)形貌SEM圖片的網(wǎng)格映射。
具體思路在一個model中建立一個目標(biāo)part1，進(jìn)行多區(qū)域劃分，然后復(fù)制建立一個與目標(biāo)part相同尺寸的無區(qū)域切分part2（模型樹里刪除切分操作即可），然后進(jìn)行網(wǎng)格劃分，最后通過自編程插件或程序?qū)崿F(xiàn)無區(qū)域切分part2網(wǎng)格到目標(biāo)part1的映射。
下面給出幾個不同的例子;
1 規(guī)則形狀純四面體網(wǎng)格
2 規(guī)則模型純?nèi)切?em>網(wǎng)格
3 規(guī)則模型四邊形和三角形混合網(wǎng)格
4 規(guī)則模型純六面體網(wǎng)格
5 規(guī)則模型純四面體網(wǎng)格
6 規(guī)則模型純鍥形體網(wǎng)格
7 規(guī)則模型六面體和鍥形體混合網(wǎng)格
8 非規(guī)則形狀二維模型網(wǎng)格映射
純四邊形
純?nèi)切?四邊形和三角形混合
9 非規(guī)則形狀三維模型網(wǎng)格映射
純六面體
純四面體
純鍥形體
六面體和鍥形體混合
最后給出總結(jié)如下圖所示;
ABAQUS斷裂模擬收徒 ，快速學(xué)會各種ABAQUS斷裂模擬方法 **/人（將有機(jī)會享有各種插件以及程序，價值**、專門定制視頻、全程親自教學(xué)、各種模型調(diào)試及解答問題等等，傾囊相教）
展開 ansys 仿真模擬 hypermesh 網(wǎng)格劃分 matlab 編程 代做/咨詢/培訓(xùn)
服務(wù)價格：
至少300及以上，具體項(xiàng)目根據(jù)難易程度、工作量具體商談
服務(wù)范圍：
Ansys 靜力模擬，振動模擬，沖擊模擬，流體模擬，復(fù)合材料分析，疲勞壽命，結(jié)構(gòu)/熱/流體/電場/磁場多場耦合、結(jié)構(gòu)優(yōu)化二次開發(fā)，培訓(xùn)；Hypermesh 復(fù)雜結(jié)構(gòu)分網(wǎng)，二次開發(fā)，培訓(xùn)；Matlab 編程，simulation 仿真等項(xiàng)目咨詢，非誠勿擾！
結(jié)構(gòu)一膠應(yīng)力.png
Out.bmp
abaqus盤式制動器的py編程語言 ￥10
<p>abaqus盤式制動器的py編程語言</p><p>可以快速導(dǎo)入abaqus中，無需建模和設(shè)置參數(shù)。</p><p>附帶py編程語言的講解</p>
基于python編程操作ABAQUS輸入文件生成PD3D單元顆粒
在空間中生成剛性顆粒（注意是剛性顆粒）有下列幾種方法：
1.修改關(guān)鍵字，構(gòu)建粒子生成器模型生成隨機(jī)分布剛性顆粒
2.使用python語言直接在ABAQUS中生成顆粒，并進(jìn)行剛體綁定，使其成為剛性顆粒，或者直接生成解析剛體或離散剛體。
方法1生成顆粒的隨機(jī)性較好，操作簡單。方法2直接在ABAQUS界面生成顆粒，當(dāng)所需顆粒數(shù)量以萬為計量單位時，在前處理界面時就會卡死，對顯卡要求極高。因此，在僅考慮到這些弊端情況下，就已經(jīng)使研究人員頭皮發(fā)麻，無從下手。
在一些特定應(yīng)用場合下，比如所需顆粒數(shù)量數(shù)以萬計，我們只能采用方法1生成顆粒，但我們不僅僅是需要顆粒，還需將這些顆粒與其它模型進(jìn)行耦合求解計算，這個時候粒子生成器就會有局限性。此外，考慮到顆粒在空間中排布的多樣性，比如最典型的高斯分布，那么粒子生成器很難做到一步到位生成所需分布特征的顆粒。
我們今天介紹的通過python編程操作ABAQUS輸入文件生成PD3D單元顆粒，其可操作性更強(qiáng)，我們可以不采用粒子生成器內(nèi)部定義的隨機(jī)算法生成顆粒，用戶可以根據(jù)需求自定義顆粒分布算法，以契合實(shí)際工況。此外，可省去粒子生成顆粒的分析步，直接進(jìn)行工況建模求解計算。
本貼只是個人興趣，只提供思路，不提供源碼，用戶需了解ABAQUS的inp文件的書寫規(guī)則、python操作文件語法和生成顆粒的底層邏輯(分布模型)。感興趣的可以私信，提供編寫思路。
下面我們采用這一方法生成直徑2mm、3mm、4mm和5mm的混合顆粒，數(shù)量為1000。具體生成結(jié)果如下圖所示。
展開 
技術(shù)鄰周報Q16：CAE編程/Abaqus/傅里葉/Python/螺紋/NVH/結(jié)構(gòu)/Fluent...
10、Abaqus+PyQt+Python平面變形歐拉角計算
作者：
jianghu
鏈接：https://www.yqgqt.org.cn/content/post/1824899
在剛體運(yùn)動學(xué)、飛機(jī)飛行、衛(wèi)星姿態(tài)等領(lǐng)域，歐拉角是一個非常重要的概念和控制參數(shù)。
11、如何利用自適應(yīng)網(wǎng)格加速Fluent仿真
作者：
安世亞太
鏈接：https://www.yqgqt.org.cn/content/post/1824847
大多數(shù)CFD模擬都是采用生成具有局部區(qū)域細(xì)化和粗化的網(wǎng)格來計算的。這些經(jīng)過細(xì)化或粗化的區(qū)域確保有足夠高的分辨率，以準(zhǔn)確捕獲重要區(qū)域位置的結(jié)果，同時也使得總網(wǎng)格數(shù)量在可控范圍內(nèi)。
12、汽車結(jié)構(gòu)開發(fā)中的常見的CAE優(yōu)化方法
作者：
luck露
鏈接：https://www.yqgqt.org.cn/content/post/1825044
當(dāng)下，汽車行業(yè)面臨巨大挑戰(zhàn)，各個車企之間的競爭，已經(jīng)由之前粗暴的增量競爭模式，演變到現(xiàn)在更加殘酷的存量競爭模式，這關(guān)系到每個車企的生死存亡，也對每個車企提出了更高的要求。一方面，車企需要不斷適應(yīng)市場的需求，加快產(chǎn)品開發(fā)的速度；同時，車企還必須提升產(chǎn)品的品質(zhì)，增強(qiáng)產(chǎn)品競爭力。這也對每個參與其中的汽車工程師提出了更高的要求。
技術(shù)鄰鼓勵創(chuàng)作者發(fā)布優(yōu)質(zhì)的文章/視頻/問答/文檔，快來發(fā)布內(nèi)容上周報吧~
展開 模塊化編程在Abaqus UMAT編寫中的應(yīng)用實(shí)例
3、模塊化編程在UMAT編程中的應(yīng)用實(shí)例
首先我們需要安裝好版本匹配的VS+IVF+Abaqus，并完成關(guān)聯(lián)。
下面我們來看一下, 一個UMAT模塊化編程例子的大體框架，我們把整個UMAT分成三大模塊，分別為前處理模塊，更新模塊，后處理模塊，我們這里僅看前處理模塊，希望能給大家一些啟發(fā)。
接著我們再來簡單介紹UMAT前處理模塊化編程常用子函數(shù)
模型參數(shù)傳遞子函數(shù)——用于將UMAT中的模型參數(shù)數(shù)組的值提取出來，并將其賦值給有物理含義的變量，這樣一來方便在UMAT編寫公式時使用，二來可以防止我們不小心在UMAT中修改模型參數(shù)的值造成錯誤。
二三維判斷子函數(shù)——在UMAT中在應(yīng)力應(yīng)變的向量的維度為NTENS，這就使得我們UMAT中應(yīng)力應(yīng)變向量的維度取決于模型是三維問題還是平面應(yīng)力（變）的問題，這不便于我們編寫具有更加普遍適用性的UMAT，因此，按照規(guī)則填滿STRESS和STRAN向量，使其維度均為1*6，這也為后面將其轉(zhuǎn)為張量形式提供便利。
向量轉(zhuǎn)張量子函數(shù)——由于UMAT中默認(rèn)的應(yīng)力應(yīng)變是存儲為向量的形式，但在本構(gòu)關(guān)系式中，我們通常采用的是張量的形式，因此將應(yīng)力應(yīng)變向量轉(zhuǎn)為為張量形式，可以讓我們更加方便將本構(gòu)關(guān)系式子用代碼實(shí)現(xiàn)，但這里特別要在處理應(yīng)變時候要記得應(yīng)變向量存儲的是工程剪切應(yīng)變，需要除以2才可以存入到應(yīng)變張量中。
雅克比張量轉(zhuǎn)換子函數(shù)——UMAT中DDSDDE(NTENS,NTENS)，但我們在本構(gòu)關(guān)系中的常用四維的雅克比矩陣，按照一定規(guī)則實(shí)現(xiàn)二者的形式上的轉(zhuǎn)換也是十分有必要的。
應(yīng)力（應(yīng)變）Invariant 計算子函數(shù)——應(yīng)力應(yīng)變張量不變量是我們在UMAT中經(jīng)常需要求解的量，因此將其編寫為統(tǒng)一接口的子函數(shù)，我們在其他UMAT中只需要輕松調(diào)用即可，節(jié)省不少工作量。
展開 從編程角度闡述有限元最佳入門方法：以Abaqus 和OpenSees 為例
劃重點(diǎn)：通過本文，初步了解有限元軟件的計算內(nèi)核和前后處理器的對計算機(jī)性能和編程語言的要求，基于此，我們應(yīng)該選擇門檻最低的GUI方法去嫻熟有限元軟件。本文通過對比學(xué)習(xí)Abaqus 和OpenSees的過程, 建議各位可以嘗試使用STKO搭載OpenSees像使用具有友好GUI 的商業(yè)有限元一樣使用OpenSEES, 而無需任何編程能力。（本文由心塵軒博士執(zhí)筆，就讀于香港理工大學(xué)，如果各位對STKO 有興趣，可加入群752746704，相互交流學(xué)習(xí)）
1. 前言：
作者曾執(zhí)筆《無需TCL編程能力，STKO帶你輕松玩轉(zhuǎn)OpenSEES》，如果有興趣可以點(diǎn)擊下
面連接查看（https://www.jishulink.co m /content/post/1848287 ），在上一篇推文中，作
者翔實(shí)論述了OpenSees混合編程（Tcl+CPP）的本質(zhì)，并用一個案列論證了STKO 可以使我
們輕松便利的玩轉(zhuǎn)Opensees，而無需任何的編程能力。
2. 正文：
本文則站在更高的維度來橫向?qū)Ρ绕渌虡I(yè)有限元如Abaqus 的學(xué)習(xí)方法，以此來論述為什么建議初學(xué)者可以借助具有友好GUI界面的STKO 來掌握OpenSEES的學(xué)習(xí)和使用.
在學(xué)習(xí)有限元前，我們必然已經(jīng)掌握了基本的有限元理論，這對我們更好地了解有限元軟件時大有裨益的。本質(zhì)上，任何有限元應(yīng)該包括三部分內(nèi)容：與求解相關(guān)模塊，前處理模塊，后處理模塊。其中決定一個有限元優(yōu)劣的就是其求解模塊，因此可以說，有限元軟件的CPU就是其求解模塊。
展開 abaqus子程序開發(fā)：fortran-c-python混合編程開發(fā)
abaqus子程序開發(fā)語言不局限fortran和c了，還能支持python！
這是一套自研的開發(fā)框架，源于客戶需要獲取分析過程的中間數(shù)據(jù)作為python機(jī)器學(xué)習(xí)程序的輸入這樣的需求。這種需求乍一聽確實(shí)非主流、非常規(guī)。
后來仔細(xì)做了一些了解，這種需求有其合理性。第一，python在數(shù)據(jù)分析、機(jī)器學(xué)習(xí)方面有豐富的資源；第二，python程序嵌入abaqus求解器一起運(yùn)行，那開發(fā)就會很方便高效；第三，借助python的pdb可以隨心所欲的設(shè)置斷點(diǎn)，調(diào)試程序。在此之前，python一般作為獨(dú)立程序，依賴數(shù)據(jù)文件的順序傳遞，實(shí)現(xiàn)與abaqus的協(xié)同。直白的講，abaqus算完了，形成數(shù)據(jù)文件再導(dǎo)入python。
于是火力全開研究解決方案，結(jié)果可謂“大快人心”：實(shí)現(xiàn)了abaqus的子程序?qū)ython的支持，而且還實(shí)現(xiàn)了fortran、c、python三種語言的混合編程，真正讓不同的編程語言在它最擅長的領(lǐng)域發(fā)揮作用。
以上介紹了abaqus子程序fortran-c-python混合編程，如有這方面的需求，歡迎私信聯(lián)系開展合作。
展開 Abaqus,Ansys等模型計算完成后自動提醒 無需編程
使用場景
1、在工作站或其它電腦上運(yùn)行程序
2、程序運(yùn)行時間長
現(xiàn)有功能
1、 監(jiān)測Abaqus單任務(wù)運(yùn)行情況：在當(dāng)前運(yùn)行的Job計算完成或中斷時發(fā)送提醒消息；
2、 監(jiān)測Abaqus多任務(wù)運(yùn)行情況：在批量模型均計算完成或中斷時發(fā)送提醒消息；
3、 監(jiān)測文件夾大小：當(dāng)使用其他軟件如Ansys、Python、Matlab工作時，本軟件可以通過工作目錄文件夾大小來對程序運(yùn)行進(jìn)度監(jiān)測，當(dāng)工作目錄文件大小不變時發(fā)送提醒消息。
使用方法
將本軟件安裝在被檢測程序運(yùn)行的電腦，選擇需要被監(jiān)測的工作目錄和任務(wù)模式，輸入接受提醒的郵箱號（微信可綁定QQ郵箱接受消息），憑使用碼進(jìn)行使用。
開啟TaskReminder軟件，點(diǎn)擊選擇路徑按鈕，選擇Abaqus工作目錄，或需要監(jiān)測狀態(tài)的文件夾路徑。
點(diǎn)選任務(wù)模式，選擇需要進(jìn)行的任務(wù)類型。
輸入接收消息提醒的郵箱地址，推薦采用qq郵箱（在微信中可以幾乎無延遲地收到消息提醒），輸入軟件使用碼。
點(diǎn)擊測試按鈕，測試軟件與服務(wù)器的連接狀態(tài)。
點(diǎn)擊開始任務(wù)，即可實(shí)時監(jiān)測軟件運(yùn)行狀態(tài)。任務(wù)完成后收到提醒。
軟件使用界面
測試后，連接正常，可以啟動任務(wù)
軟件正在運(yùn)行中，可以監(jiān)測到Abaqus計算進(jìn)度。
任務(wù)完成后收到郵箱提醒
軟件試用方式
軟件正在處于開放試用階段，歡迎嘗試體驗(yàn)。私信獲取軟件下載鏈接與使用碼。
展開 Abaqus利用python編程進(jìn)行仿真結(jié)果(odb)數(shù)據(jù)圖片批量處理的案例講解
續(xù)集(一維彈簧單元的直接剛度方法)Python編程和ABAQUS結(jié)果對比
我們編程的結(jié)果如下所示
ABAQUS結(jié)果和編程結(jié)果對比
所以ABAQUS提取的整體剛度矩陣實(shí)際是經(jīng)過補(bǔ)償后的剛度矩陣，嚴(yán)格來說是不正確的，但是并不影響力和位移的結(jié)果。
正確的剛度矩陣如下
如果疑問，歡迎交流和指正。
補(bǔ)償法的部分代碼如下
(一維彈簧單元的直接剛度方法)Python編程和ABAQUS結(jié)果對比
在學(xué)習(xí)《有限元方法基礎(chǔ)教程》過程中，通過自己編程實(shí)現(xiàn)有限單元法的數(shù)值解答，加深對理論的理解。這個在我去年的時候就已經(jīng)發(fā)過帖子了，但是沒有講解代碼，沒有和ABAQUS有限元軟件對比。這次發(fā)帖子就是想講解代碼以及和ABAQUS結(jié)果對比，發(fā)現(xiàn)我又更深層次理解了有限元方法。
這次寫得是最簡單的模型：一維線彈簧單元。采用的是直接剛度法求解。
直接剛度法的求解思路如下所示，其中粉紅色的是輸入，淡藍(lán)色是輸出。主要是要區(qū)分齊次邊界條件和非齊次邊界條件，非齊次邊界條件的話就要修改【F】。
下面將貼出我用python寫得一維彈簧單元的直接剛度法：
例子計算：
如圖是一個彈簧系統(tǒng)，單元節(jié)點(diǎn)信息如下，5節(jié)點(diǎn)受到一個強(qiáng)制位移20mm，明顯這是一個非齊次邊界條件問題。
Python編程輸入信息如下：
結(jié)果如下：
可以看到，輸出結(jié)果和書上的答案一致。
下面進(jìn)行ABAQUS模擬：
添加彈簧單元
添加邊界條件
顯示單元編號、節(jié)點(diǎn)編號如下所示，紅色表示單元編號，黃色是節(jié)點(diǎn)編號
ABAQUS計算結(jié)果如下：
首先是變形圖前后對比
反力云圖如下所示，基本和直接剛度法計算的結(jié)果一致
位移云圖如下所示，基本和直接剛度法計算的結(jié)果一致
整體剛度矩陣如下所示，因?yàn)?em>ABAQUS彈簧單元是三維的，每一個節(jié)點(diǎn)有3個自由度，15x15，原味的剛度矩陣如下
我們把剛度矩陣轉(zhuǎn)化為一維的，方便和編程的結(jié)果對比
從結(jié)果可以看出，ABAQUS的整體剛度矩陣和直接剛度法計算出來的整體剛度矩陣有些差異，如圖標(biāo)紅的所示。
那么在整體剛度矩陣上為什么ABAQUS會和直接剛度法的整體剛度矩陣有差異呢，到底ABAQUS的整體剛度矩陣對不對呢？答案將在下一期揭曉。歡迎大家積極討論。
展開 Abaqus:利用python編程進(jìn)行仿真結(jié)果（odb）數(shù)據(jù)圖片批量處理案例講解（下）
abaqus做DCB模擬時如何導(dǎo)出應(yīng)變能釋放率曲線，可以不編程嗎？
寫論文需要R曲線（斷裂韌性-裂紋長度），但是不會編程
基于ABAQUS子程序UAMP編程實(shí)現(xiàn)水平井分段多簇壓裂流量動態(tài)分配
利用ABAQUS平臺的二次開發(fā)功能，通過Fortran語言進(jìn)行用戶子程序UAMP的編程，求解流量分配控制方程。UAMP可用來定義當(dāng)前幅值隨任意函數(shù)的變化，每一個增量步開始前，縫內(nèi)流體壓力將通過編程指令傳入子程序中，用來求解方程（7）和（8）得到每條裂縫注入速度，求解之后將結(jié)果傳遞到ABAQUS主程序中繼續(xù)求解流—固耦合方程。
在最初始的階段，各條裂縫都會注入等量的壓裂液。但隨著裂縫的擴(kuò)展，在應(yīng)力干擾的作用下壓力阻力發(fā)生變化，下一階段會吸收不同分量的壓裂液，隨著泵送的持續(xù)進(jìn)行，大裂縫的擴(kuò)展速度將會加快，小裂縫的擴(kuò)展速度將會相對減慢，導(dǎo)致它們的形態(tài)出現(xiàn)差異。
三、有限元模型：
本模型基于前述理論分析，根據(jù)平面應(yīng)變假設(shè)，利用擴(kuò)展有限元法（XFEM）建立三條裂縫同步擴(kuò)展數(shù)值模型，見圖4。模型尺寸為400m×200m，四周位移邊界固定，地層參數(shù)以吉木薩爾蘆草溝組某井?dāng)?shù)據(jù)為例。
圖4水平井多裂縫擴(kuò)展數(shù)值模型
計算分為兩步，第一步平衡地應(yīng)力，模擬儲層初始的賦存狀態(tài)，第二步，以總排量為9m3/min的速度注入壓裂液模擬水力壓裂過程。
四、計算結(jié)果分析
儲層的孔隙壓力分布及裂縫擴(kuò)展形態(tài)如圖5所示，同時提取了壓力—流量曲線如圖6所示。由圖可知在初始階段，內(nèi)側(cè)裂縫的縫內(nèi)壓力稍大于外側(cè)裂縫，對于壓裂液的阻力較大，使外側(cè)裂縫的注入流量大于內(nèi)側(cè)裂縫，隨著壓裂液注入，內(nèi)外裂縫的壓力出現(xiàn)明顯差異，內(nèi)側(cè)裂縫壓力遠(yuǎn)大于外側(cè)裂縫，而注入流量遠(yuǎn)小于外側(cè)裂縫，在注入流量差異和縫間干擾的作用下，內(nèi)側(cè)裂縫受到抑制，外側(cè)裂縫擴(kuò)展較快而發(fā)生偏轉(zhuǎn)，三縫同步擴(kuò)展時壓力與流量成反比例關(guān)系，與流量分配控制的原理規(guī)律相符，因此證明了子程序編程的正確性。
當(dāng)裂縫近間距擴(kuò)展時，由于縫間干擾應(yīng)力，使得地應(yīng)力場發(fā)生不同程度的偏轉(zhuǎn)，中間裂縫受到抑制，外側(cè)裂縫發(fā)生偏轉(zhuǎn)現(xiàn)象。
展開 