abaqus計算編程的案例
Abaqus,Ansys等模型計算完成后自動提醒 無需編程
前言
你在工作站上批量算模型的時候是否會頻頻去檢查計算進度？
你是否有過信心滿滿提交計算作業，結果過段時間回來看第一步就不收斂？
你在趕ddl時是不是有著“人可以休息，電腦不可以休息”的心態？
如果您曾遇到過以上的煩惱，
TaskReminder_v1.0或許可以幫助你更加高效地進行計算任務。
軟件用途
監測程序運行情況，在程序完成或中斷時本軟件會通過郵件發送提醒。
使用場景
1、在工作站或其它電腦上運行程序
2、程序運行時間長
現有功能
1、 監測Abaqus單任務運行情況：在當前運行的Job計算完成或中斷時發送提醒消息；
2、 監測Abaqus多任務運行情況：在批量模型均計算完成或中斷時發送提醒消息；
3、 監測文件夾大小：當使用其他軟件如Ansys、Python、Matlab工作時，本軟件可以通過工作目錄文件夾大小來對程序運行進度監測，當工作目錄文件大小不變時發送提醒消息。
使用方法
將本軟件安裝在被檢測程序運行的電腦，選擇需要被監測的工作目錄和任務模式，輸入接受提醒的郵箱號（微信可綁定QQ郵箱接受消息），憑使用碼進行使用。
開啟TaskReminder軟件，點擊選擇路徑按鈕，選擇Abaqus工作目錄，或需要監測狀態的文件夾路徑。
點選任務模式，選擇需要進行的任務類型。
輸入接收消息提醒的郵箱地址，推薦采用qq郵箱（在微信中可以幾乎無延遲地收到消息提醒），輸入軟件使用碼。
點擊測試按鈕，測試軟件與服務器的連接狀態。
點擊開始任務，即可實時監測軟件運行狀態。任務完成后收到提醒。
軟件使用界面
測試后，連接正常，可以啟動任務
軟件正在運行中，可以監測到Abaqus計算進度。
展開 大規模工程計算解決方案-并行計算編程高級培訓班
中國科學院計算技術研究所教育中心
隨著CAE技術的高速發展，工程計算的規模逐漸擴大，以致多核計算成為發展趨勢。目前各研究機構多核并行程序開發不足，隨著多核硬件技術發展，在硬件設備上存在不同程度的利用不足。通過加強并行程序的設計與開發，研究機構可以充分的利用現有設備資源提高工作效率，降低時間成本，實現更大規模的計算。中國科學院計算技術研究所教育中心依托自身強大優勢，特舉辦“大規模工程計算解決方案-并行計算編程”高級培訓班，相關具體事宜通知如下 ：
一、培訓目標：
1、現有硬件資源最大化利用
2、提高并行程序設計與開發水平；
3、實際體驗超級計算環境；
4、提供符合自身特點的工程計算解決方案：
三、課程實操：
學員體驗曙光6000高性能計算機環境，在改并行機群上編譯、調試和運行教學例子，并在最多達128個核規模下測試和分析程序的擴展性和性能。
四、時間地點:2012年11月29日
— 2012年12月02日
北 京（時間安排：第一天報到、授課三天）
課
程
內
容
第一章：基本的MPI編程技術
1．高性能計算概述
2．簡單的MPI程序設計
3．如何運行MPI程序
上機實操課程
應用實例
1．MPICH的安裝、程序編譯與運行。
2．Hello world程序的編寫
3. 循環消息傳遞程序的編寫。
第二章：高級MPI編程技術
1．定義新的數據類型
2．非阻塞通信
3．各種組通信
上機實操課程
應用實例
1．4種產生新數據類型的方法
；
2．各種非阻塞的發送與接收操作
3 最常見的幾種組通信操作。
展開 Julia：高效易用的數值計算/優化編程語言
三、交互性編程
我們指出，類似于IPython，Julia中的IJulia package也可以讓Julia的所有編譯過程在Jupyter notebook里進行。我個人是很喜歡在這個環境里進行Julia編程的（詳見：JuliaLang/IJulia.jl）。當然，我知道也有不少人喜歡Juno的：Juno, the Interactive Development Environment 可能這個更有碼農的感覺吧hhh
進一步的，我們可以有很多交互式的操作，這在Julia中主要通過Interact package實現。比如，我們可以自定義slider，按鈕等對一個參數曲線進行互動。一個例子見如下視頻，或者JuliaGizmos/Interact.jl。
Interact in Julia
vimeo.com
至于各種畫圖，我傾向于使用Plots這個package。入門可以見：Plots - powerful convenience for visualization in Julia.
四、寫在最后
自然，本文給出的只是很少的一些例子和對Julia這門編程語言的最基本的介紹。無論你只是希望有個方便的語言調用solver，或者做數值計算居多，還是比較高級的優化算法專家，多實踐，一邊"get hands dirty"一邊學習我覺得總是最有效率的。
首當其沖的是Julia官網上提供的大量學習資源：包括視頻，具體的算例，和各種pdf教程。
Learning Julia3 ulialang.org
作為Julia cofounder之一的Prof.
展開 matlab編程計算矩形薄板的固有頻率
程序目錄：
rectangularplate.m——主程序
rectangularMeshKirchhoff.m——網格劃分子程序
formStiffnessMassKirchhoff.m——剛度、質量陣計算及組裝子程序
EssentialBCKirchhoff.m——施加邊界條件子程序
本程序可以計算各種邊界條件下的矩形板薄板固有頻率和對應振型。主程序例子為懸臂板，
其中邊界條件可更改，即第23行的'cfff' 可更改為：’ssss’、’cccc’、’scsc’、’cccf’、’cfff’
s為簡支，c為固支，f為自由，四邊的順序為左下右上。
程序已經驗證正確。
rectangularplate.rar
展開 
CNC數控編程加工常用計算公式，有學到嗎？
一、 三角函數計算
1.tanθ=b/a θ=tan-1b/a
2.Sinθ=b/c Cos=a/c
二、切削速度的計算
Ｖc＝(π*Ｄ*Ｓ)／１０００
Vc：線速度(m/min) π：圓周率(3.14159) D：刀具直徑(mm) S：轉速(rpm)
例題. 使用Φ25的銑刀Vc為(m/min)25求S=?rpm
Vc=πds/1000
25=π*25*S/1000
S=1000*25/ π*25
S=320rpm
三、進給量(F值)的計算
Ｆ＝Ｓ＊Ｚ＊Ｆz
F：進給量(mm/min) S：轉速(rpm) Z：刃數
Ｆz：(實際每刃進給)
例題.一標準2刃立銑刀以2000rpm)速度切削工件，求進給量(F值)為多少？(Ｆz=0.25mm)
Ｆ＝Ｓ＊Ｚ＊Ｆz
Ｆ＝2000*2*0.25
Ｆ＝1000(mm/min)
四、殘料高的計算
Scallop＝（ae＊ae）／8R
Scallop：殘料高(mm) ae：XYpitch(mm) R刀具半徑(mm)
例題.Φ20R10精修2枚刃，預殘料高0.002mm，求Pitch為多
少？mm
Scallop=ae2/8R
0.002=ae2/8*10
ae=0.4mm
五、逃料孔的計算
Φ＝√２Ｒ2 Ｘ、Ｙ＝Ｄ／４
Φ：逃料孔直徑(mm) R刀具半徑(mm) D：刀具直徑(mm)
例題. 已知一模穴須逃角加工(如圖)，
所用銑刀為ψ10；請問逃角孔最小
為多少？圓心坐標多少？
展開 數控編程計算公式，趕快收藏，保證會用得上！
一、 三角函數計算
1.tanθ=b/a θ=tan-1b/a
2.Sinθ=b/c Cos=a/c
二、切削速度的計算
Ｖc＝(π*Ｄ*Ｓ)／１０００
Vc：線速度(m/min) π：圓周率(3.14159) D：刀具直徑(mm)
S：轉速(rpm)
例題. 使用Φ25的銑刀Vc為(m/min)25求S=?rpm
Vc=πds/1000
25=π*25*S/1000
S=1000*25/ π*25
S=320rpm
三、進給量(F值)的計算
Ｆ＝Ｓ＊Ｚ＊Ｆz
F：進給量(mm/min) S：轉速(rpm) Z：刃數
Ｆz：(實際每刃進給)
例題.一標準2刃立銑刀以2000rpm)速度切削工件，求進給量(F
值)為多少？(Ｆz=0.25mm)
Ｆ＝Ｓ＊Ｚ＊Ｆz
Ｆ＝2000*2*0.25
Ｆ＝1000(mm/min)
四、殘料高的計算
Scallop＝（ae＊ae）／8R
Scallop：殘料高(mm) ae：XYpitch(mm) R刀具半徑(mm)
例題.Φ20R10精修2枚刃，預殘料高0.002mm，求Pitch為多
少？mm
Scallop=ae2/8R
0.002=ae2/8*10
ae=0.4mm
五、逃料孔的計算
Φ＝√２Ｒ2 Ｘ、Ｙ＝Ｄ／４
Φ：逃料孔直徑(mm) R刀具半徑(mm) D：刀具直徑(mm)
例題.
展開 數控編程課堂：關于螺紋螺距的含義和計算方式，太全了！
例2：銑15齒的齒輪計算。
代入公式：n=40/15
計算n=2(2/3)
是轉2整圈再選被3整除的分度孔如24，30，39，42.51.54.57，66等孔板上再向前轉過16，20，26，28，34，36，38，44加1孔即17，21，27，29，35，37，39，45孔作為本輪的定位孔。
例3：銑82齒的分度計算。
代入公式：n=40/82
計算n=20/41
即：只要選41孔的分度板，在上輪定位孔上再轉過20+1即21孔作為本輪的定位孔便是。
例4：銑51齒的分度計算
代入公式n=40/51由于此時分數無法計算則只能直接選孔，即選51孔的分度板，在上輪定位孔上再轉過51+1即52孔作為本輪的定位孔即是。
例5：銑100齒的分度計算。
代入公式n=40/100
計算n=4/10=12/30
即選30孔的分度板，在上輪定位孔上再轉過12+1即13孔作為本輪的定位孔即是。
如所有分度板無計算所需的孔數時則應采用復式分度法計算，不在本計算方法之列，實際生產時一般采用滾齒，因復式分度計算后的實際操作極為不便。
展開 有限元基礎編程 ｜高階單元計算環形區域慣性矩
但若是模型較大時，一味的為了追求精度而增加網格數量，勢必會引起計算成本的增加，計算成本與數值精度的取舍，是個值得思考的問題。
本文的主要參考內容及Matlab代碼均來自Amar Khennane編寫的《Introduction to Finite Element Analysis Using MATLAB and Abaqus》，想要進一步了解有限元編程的小伙伴可以入手，強烈推薦！
完整代碼可在公眾號后臺回復：8節點單元，即可自動獲取。
本次分享僅限于此了，歡迎大家點贊收藏轉發！
謝謝你看完木木同學的分享，今日份閱讀花費的流量+1M哈哈哈哈哈哈。
-End-
易木木響叮當
想陪你一起度過短暫且漫長的科研生活
展開 mathcad14編程實例補充之6—“橋梁板式橡膠支座選用計算”
這是mathcad14編程實例補充之6—“橋梁板式橡膠支座選用計算”
為的是補充《編程淺說》中的“編程實例”部分，使其包含的領域更廣泛
一些。
這是編程實例部分的第16節，“橋梁板式橡膠支座選用計算”。
橋梁板式橡膠支座選用計算1.jpg
橋梁板式橡膠支座選用計算2.jpg
橋梁板式橡膠支座選用計算3.jpg
橋梁板式橡膠支座選用計算4.jpg
板式橡膠支座選用計算.rar
PLC編程必備：模擬量計算，結合實際案例詳解，這操作一看就懂
那么TIA Portal博途中模擬量輸入/輸出的編程
1）SCALE/UNSCALE在哪里？
在TIA Portal中指令SCALE/UNSCALE的功能和STEP7中的FC105/FC106功能相同，用于模擬量輸入/輸出的編程。
編程界面下，在 "指令" 中的"基本指令 ——> 轉換操作"中可以找到，
abaqus盤式制動器的py編程語言 ￥10
<p>abaqus盤式制動器的py編程語言</p><p>可以快速導入abaqus中，無需建模和設置參數。</p><p>附帶py編程語言的講解</p>
基于python編程操作ABAQUS輸入文件生成PD3D單元顆粒
在空間中生成剛性顆粒（注意是剛性顆粒）有下列幾種方法：
1.修改關鍵字，構建粒子生成器模型生成隨機分布剛性顆粒
2.使用python語言直接在ABAQUS中生成顆粒，并進行剛體綁定，使其成為剛性顆粒，或者直接生成解析剛體或離散剛體。
方法1生成顆粒的隨機性較好，操作簡單。方法2直接在ABAQUS界面生成顆粒，當所需顆粒數量以萬為計量單位時，在前處理界面時就會卡死，對顯卡要求極高。因此，在僅考慮到這些弊端情況下，就已經使研究人員頭皮發麻，無從下手。
在一些特定應用場合下，比如所需顆粒數量數以萬計，我們只能采用方法1生成顆粒，但我們不僅僅是需要顆粒，還需將這些顆粒與其它模型進行耦合求解計算，這個時候粒子生成器就會有局限性。此外，考慮到顆粒在空間中排布的多樣性，比如最典型的高斯分布，那么粒子生成器很難做到一步到位生成所需分布特征的顆粒。
我們今天介紹的通過python編程操作ABAQUS輸入文件生成PD3D單元顆粒，其可操作性更強，我們可以不采用粒子生成器內部定義的隨機算法生成顆粒，用戶可以根據需求自定義顆粒分布算法，以契合實際工況。此外，可省去粒子生成顆粒的分析步，直接進行工況建模求解計算。
本貼只是個人興趣，只提供思路，不提供源碼，用戶需了解ABAQUS的inp文件的書寫規則、python操作文件語法和生成顆粒的底層邏輯(分布模型)。感興趣的可以私信，提供編寫思路。
下面我們采用這一方法生成直徑2mm、3mm、4mm和5mm的混合顆粒，數量為1000。具體生成結果如下圖所示。
展開 技術鄰周報Q16：CAE編程/Abaqus/傅里葉/Python/螺紋/NVH/結構/Fluent...
5、基于OpenSees和SAP2000靜力動力計算案例分析
作者：
建源之光
鏈接：https://www.yqgqt.org.cn/content/post/1824247
面對復雜的建筑的設計，它是由無數的梁、柱、板、墻最基本的單元組成，其中最基本的結構是一榀門式框架，在我們對結構仿真中，軟件其實只是工具，真正對我們設計負責的還是對力學、結構變形機理的認識，這是我們最應該掌握的，這樣才能明白軟件是怎么算的，對我們的仿真結果才可以加以正確的判斷。
6、Abaqus預應力模態分析
作者：裴一鳴
鏈接：https://www.yqgqt.org.cn/content/post/1824278
模態分析是一個線性攝動分析，只能進行線性求解。在動力學方程中，其載荷矩陣和阻尼矩陣為0，特征值的提取只取決于剛度矩陣和質量矩陣。而結構在外載荷的作用下剛度矩陣會發生變化，也就間接影響了結構的固有頻率。而預應力狀態下，我們不清楚剛度矩陣的變化對模態頻率的影響時，便需要進行預應力模態分析。
7、iSolver案例分享：波紋鋼腹板簡支梁受彎分析
作者：
Infinite_9882
鏈接：https://www.yqgqt.org.cn/content/post/1824652
結構有限元求解器iSolver已發展到一定階段，現采用結構有限元軟件iSolver進行結構分析，iSolver可使用Abaqus作為前后處理工具，本帖以波紋鋼腹板簡支梁彎曲分析為例，將iSolver求解器和Abaqus計算結果進行對比，比對兩種有限元軟件的計算結果。
展開 模塊化編程在Abaqus UMAT編寫中的應用實例
3、模塊化編程在UMAT編程中的應用實例
首先我們需要安裝好版本匹配的VS+IVF+Abaqus，并完成關聯。
下面我們來看一下, 一個UMAT模塊化編程例子的大體框架，我們把整個UMAT分成三大模塊，分別為前處理模塊，更新模塊，后處理模塊，我們這里僅看前處理模塊，希望能給大家一些啟發。
接著我們再來簡單介紹UMAT前處理模塊化編程常用子函數
模型參數傳遞子函數——用于將UMAT中的模型參數數組的值提取出來，并將其賦值給有物理含義的變量，這樣一來方便在UMAT編寫公式時使用，二來可以防止我們不小心在UMAT中修改模型參數的值造成錯誤。
二三維判斷子函數——在UMAT中在應力應變的向量的維度為NTENS，這就使得我們UMAT中應力應變向量的維度取決于模型是三維問題還是平面應力（變）的問題，這不便于我們編寫具有更加普遍適用性的UMAT，因此，按照規則填滿STRESS和STRAN向量，使其維度均為1*6，這也為后面將其轉為張量形式提供便利。
向量轉張量子函數——由于UMAT中默認的應力應變是存儲為向量的形式，但在本構關系式中，我們通常采用的是張量的形式，因此將應力應變向量轉為為張量形式，可以讓我們更加方便將本構關系式子用代碼實現，但這里特別要在處理應變時候要記得應變向量存儲的是工程剪切應變，需要除以2才可以存入到應變張量中。
雅克比張量轉換子函數——UMAT中DDSDDE(NTENS,NTENS)，但我們在本構關系中的常用四維的雅克比矩陣，按照一定規則實現二者的形式上的轉換也是十分有必要的。
應力（應變）Invariant 計算子函數——應力應變張量不變量是我們在UMAT中經常需要求解的量，因此將其編寫為統一接口的子函數，我們在其他UMAT中只需要輕松調用即可，節省不少工作量。
展開 abaqus子程序開發：fortran-c-python混合編程開發
abaqus子程序開發語言不局限fortran和c了，還能支持python！
這是一套自研的開發框架，源于客戶需要獲取分析過程的中間數據作為python機器學習程序的輸入這樣的需求。這種需求乍一聽確實非主流、非常規。
后來仔細做了一些了解，這種需求有其合理性。第一，python在數據分析、機器學習方面有豐富的資源；第二，python程序嵌入abaqus求解器一起運行，那開發就會很方便高效；第三，借助python的pdb可以隨心所欲的設置斷點，調試程序。在此之前，python一般作為獨立程序，依賴數據文件的順序傳遞，實現與abaqus的協同。直白的講，abaqus算完了，形成數據文件再導入python。
于是火力全開研究解決方案，結果可謂“大快人心”：實現了abaqus的子程序對python的支持，而且還實現了fortran、c、python三種語言的混合編程，真正讓不同的編程語言在它最擅長的領域發揮作用。
以上介紹了abaqus子程序fortran-c-python混合編程，如有這方面的需求，歡迎私信聯系開展合作。
展開 