三維結冰程序的案例
—Spring3D-Drop結冰程序研發記
1 Spring3D應該長啥樣
我在之前一篇文章中提到過，要啟動三維結冰程序的研發，這是20年6月份的事情。
結冰程序，有三個關鍵的技術：流場計算、水滴收集系數計算、冰形計算。而在三維結冰這個方面，一般能完成前兩個，就已經達到了相當高的程度，并具備了很強的實用性。
在研發之前，一個非常重要的問題就是要確定，這個軟件應該長啥樣？這個問題倒不是說界面炫不炫，而是整個軟件的設計理念是什么？也就是所謂的概念設計。
基于前面二維結冰程序開發的經驗，還有我自己的從業經歷，我給這個三維軟件一些要求如下：
（1）快。真男人就是要快。好軟件就是要快?？炀鸵馕吨杀镜?，意味著效率高，意味著設計迭代短。快就是正義。部件級別的完成計算，用時2小時以內，整機級別的計算，控制在4-6小時。
（2）具備可用的精度。工程上能夠接受的精度，這個是軟件的底線。
（3）程序模塊化。大到流場計算模塊，小到一個判斷單元位置的算法，通通模塊化，設置為自定義函數，在主程序里面調用。這樣主程序簡潔、清晰，子程序可以隨時優化更新，不斷升級。
（4）設計輸入規范。根據軟件自身特點，確定輸入文件的格式，并具有普適性。就是說，只要按照這個格式和要求提供輸入，就能成功計算。
（5）軟件輸出可直接用于設計。比如像Fensap三維結冰軟件，計算完水滴撞擊后，如果你想在數模上畫出撞擊范圍，就需要自己后處理，取多個剖面的結果。Spring3D軟件則要根據設計研發需求，直接給出可用的結果，避免過多的人工后處理操作。
20年6月份啟動的研發，那時我仍在西北，搞飛機試飛，當月就采用三維面元法解決了流場計算問題。
展開 ABAQUS三維hill48彈塑性模型VUmat子程序（彈性為正交各向異性） ￥388
1.ABAQUS三維hill48彈塑性模型VUmat子程序
2.彈性階段為正交各項異性材料
3.hill48和正交各項異性材料參數參考ABAQUS靜力模塊自帶的模型參數
4.發貨方式為百度網盤鏈接，包含子程序及上面跑的兩個模型相關文件，包含Cae，inp文件，odb文件等
5.ABAQUS版本為2024，低版本可以利用導入inp文件的方式運行及修改
6.可以免費答疑三次，后續添加你自己的模型或者相關參數等輔導都可以優惠。
計算三維問題的應力強度因子的程序
/COM,ANSYS MEDIA REL. 60 (090601) REF. VERIF. MANUAL: REL. 60
/VERIFY,VM143
*CREATE,FRACT,MAC
/NOPR
NSEL,ALL
*GET,N,NODE,,NUM,MAX ! CURRENT MAXIMUM NODE NUMBER
CMSEL,S,CRACKTIP ! SELECT THE TIP NODES
ESLN ! ANY ELEMENTS ATTACHED
*GET,ELMAX,ELEM,,NUM,MAX ! CURRENT MAXIMUM ELEMENT NUMBER
*DO,IEL,1,ELMAX ! LOOP ON MAX ELEMENT
ELMI=IEL
*IF,ELMI,LE,0,EXIT ! NO MORE SELECTED
*GET,ELTYPE,ELEM,ELMI,ATTR,TYPE ! GET ELEMENT TYPE
*IF,ELTYPE,NE,ARG1,CYCLE ! CHECK FOR SELECTED ELEMENT
N3 = NELEM(ELMI,3) ! GET NODE 3 (K)
*IF,NSEL(N3),LE,0,CYCLE ! IT MUST BE SELECTED
N7 = NELEM(ELMI,7) ! GET NODE 7 (L)
*IF,NSEL(N7),LE,0,CYCLE ! IT MUST ALSO BE SELECTED
N1 = NELEM(ELMI,1) ! GET NODE 1 (I)
N2 = NELEM(ELMI,2) ! GET NODE 2 (J)
N5 = NELEM(ELMI,5) ! GET NODE
展開 ABAQUS 隨機生成三維柱形纖維python程序 ￥10
打開txt文件，在txt調整參數后復制黏貼到abaqus最下面即可
木材三維彈塑性本構子程序開發
子程序編寫流程
本工作室在三維hashing模型的基礎上，利用Abaqus軟件平臺，開發了完整的木材的彈塑性本構umat子程序，包含木材完整的彈性、塑性、強化以及軟化階段。編寫子程序的流程如下：
三。結果驗證
通過如下圖的木材模型進行驗證：
該模型在受壓、受剪及受拉的工況下，應力應變曲線如下所示：
該子程序還有以下特征：
能計算靜力非線性
2. 收斂性好
3. 能計算復雜應力狀態
附件為本子程序參考的文獻，供大家學習探討~
2. 木材的力學性質試驗研究及數值模擬方法.pdf
最后，大家有關于Abaqus二次開發的相關需求可以添加管理員扣扣：3045552826，微信：CAE320，同時也歡迎大家關注“320科技工作室”的微信公眾號，掃一掃二維碼即可關注~~
展開 二維三維彈塑性斷裂相場------uel子程序 ￥69.9
包含二維三維彈塑性斷裂相場程序，原始代碼，以及對應的程序公式，job文件，計算收斂性較好，運算穩定，可使用二維的三角形，四邊形單元，以及三維四面體，六面體單元的程序計算案例的結果：
二維：
二維斷裂相場分布：
三維斷裂相場模型：
一個完整的三維移動高斯分布表面熱源焊接例子（附子程序） ￥2.4
完整的三維移動高斯分布表面熱源焊接例子，完整的子程序在“付費后”附件中。
基于Abaqus的umat子程序實現三維編織復合材料疲勞性能預測
后記
最后，筆者在這邊強調下，本文提到的方法適用面較廣，不僅僅只針對于三維編織復合材料，對于一般的纖維增強復合材料均可適用，只需在Umat子程序中調整對應參數即可。同時，大家若有更好的想法、建議或相關需求，可通過公眾號或站內私信聯系我們。
歡迎大家關注“320科技工作室”微信公眾號，
參考文獻：
馮繼強. 三維四向編織復合材料疲勞性能研究[D]. 2017.
ABAQUS焊接模擬-空間三維多路徑串行焊接（Fortran子程序二次開發）
本文通過ABAQUS熱傳導方式講解空間三維多路徑順序焊接建模過程，多路徑焊接重點在于子程序編寫上面。
詳細操作視頻講解請查看：http://www.yqgqt.org.cn/college/video/c10113?nagivator=training
本文使用的是熱傳導分析步，只獲得溫度場，溫度結果如圖所示：
如需應力場或變形，可將分析步改成熱力耦合分析步，如下圖所示，單元類型也得改成熱力耦合。
多路徑的實現可以通過多個分析步，然后子程序里面在對應的分析步里面編寫路徑，也可以使用一個分析步，通過時間控制，這些都用條件語句就可以實現。
路徑的運動可以用參數方程來表示，將參數方程寫入子程序里對應的坐標中。實現起來還是挺簡單的。
本例僅供參考，如若有錯誤，歡迎指正。本人QQ:289328659，歡迎交流。
版權所有，轉載請注明出處?。。?/span>
展開 