程序?qū)崿F(xiàn)
關注創(chuàng)建者：胖子愛學習 創(chuàng)建時間：2021-01-04
程序?qū)崿F(xiàn)的視頻教程
（未完結(jié)，謹慎下單）ABAQUS用戶單元子程序（UEL）理論推導及程序實現(xiàn)
本視頻主要介紹了ABAQUS用戶單元子程序（UEL）的理論推導和程序?qū)崿F(xiàn)過程，具體內(nèi)容如下： （一）介紹了UEL學習的一些方法和參考資料； （二）結(jié)合自己學習UEL的理解和總結(jié)簡述了什么叫做UEL以及UEL在ABAQUS中起到的作用，讓大家了解UEL實現(xiàn)的原理； （三）基于有限元理論對UEL實現(xiàn)過程中的一些公式進行推導，推導的過程基于矩陣形式，簡明易懂； （四）基于理論推導的結(jié)果，以平面應力問題的四邊形單元為例演示程序的實現(xiàn)
￥80 3小時7分鐘 2513播放
查看
程序?qū)崿F(xiàn)的實例教程
newmark迭代的MATLAB程序?qū)崿F(xiàn)
% newmark方法的程序?qū)崿F(xiàn)
clc
clear
m=[6/12,1/12;1/12,2/12];
k=[2*3,-1*2;-1*2,1*2];
f1=[0;1];
d2=inv(m)*f1;
b=0.25;
r=0.5;
dt=0.12;
a0=1/(b*(dt)^2);
a1=r/(b*dt);
a2=1/(b*dt);
a3=1/(2*b)-1;
a4=r/b-1;
a5=0.5*dt*((r/b)-2);
a6=dt*(1-r);
a7=r*dt;
pk=k+a0*m;
pk=inv(pk);
d0=[0;0];
d1=[0;0];
for i=1:10
t=0.12*i;
if t>0.5
f1(2)=0;
end
ff=0.0;
d00=0.0;
ff=ff+m*(a0*d0+a2*d1+a3*d2);
pf1=f1+ff;
d00=d00+pk*pf1;
d22=a0*(d00-d0)-a2*d1-a3*d2;
d11=d1+a6*d2+a7*d22;
d0=d00;
d1=d11;
d2=d22;
a(i)=d00(1);
b(i)=d00(2);
end
a
b
這個程序其實較好理解，先搞懂newmark方法，一般的動力學有關的書上會有newmark方法的編程思路，結(jié)合著看一下，應該十分鐘內(nèi)您就會明白．
第一段是一些初值設置；第二段就是迭代算法！
批評指正，相互學習！
轉(zhuǎn)自 振動論壇
展開 結(jié)合論壇上各位大蝦的經(jīng)驗和自帶幫助,前段時間動手試了一下用戶子程序的實現(xiàn)過程,現(xiàn)將自己的一些體會與大家共享,相信對做二次開發(fā)的新手有用
N;x
用戶子程序?qū)崿F(xiàn)過程:
第一步：所需文件準備：把DEF_SIM目錄下的def_usr.f , DEF_SIM_USR_ABsoft70.gui , DEF_SIM_P4_USR_ABsoft70.gui , DEF_SIM_P4_USR_LIB.lib , DEF_SIM_USR_LIB.lib和lib目錄下的所有.lib文件拷到工作目錄下。備份DEF_SIM.exe ，DEF_SIM_P4.exe兩個文件。
第二步：編寫源程序。對于新手可以使用模版自帶的子程序，對于材料流動應力子程序模版里頭本身就自帶有一個，所以可以不需要自己編寫。
第三步：用absoft7.0或以上版本打開DEF_SIM_USR_ABsoft70.gui 文檔，指定好library files，用工作目錄下的lib文件來替換原來默認的lib文件（這一步很關鍵，試了好久才發(fā)現(xiàn)），要不然編譯可以通過，可是不能生成執(zhí)行文件。上述做完后直接點build就ok了，自動生成了DEF_SIM.exe。重復上面的過程打開DEF_SIM_P4_USR_ABsoft70.gui 文檔生成DEF_SIM_P4.exe。
第四步：用生成的DEF_SIM.exe ，DEF_SIM_P4.exe替換原來的這兩個文件。
第五步，運行模型。對于流動應力子程序，只要在material那里選擇子程序就ok了。
展開 （1）在Umat子程序中選用三維Hashin準則進行纖維束的失效判斷，選用最大應力準則進行基體的失效判斷。
（2）確定單胞疲勞失效載荷、極限應變，和纖維束平均應力。
（3）計算時，計算模型需對界面和邊界條件采用一定的假設與簡化。提取加載點位移載荷數(shù)據(jù)。
圖1 單胞靜力分析載荷位移曲線
2.2 疲勞分析
（1）建立纖維束、界面、基體的疲勞折減模型，定義損傷因子，這些均在Umat子程序中加以實現(xiàn)。同樣的，在建立疲勞折減模型時也需引入一些的假設以簡化模型。在本此研究模型中，首先利用以下文獻中提供的單向復合材料疲勞壽命預測公式確定不同應力水平下纖維束的疲勞壽命：
（2）依據(jù)計算所得的疲勞壽命，在Umat子程序中實現(xiàn)不同應力水平下纖維束、界面和基體的剛度折減、強度折減。
（3）最后，利用子程序?qū)崿F(xiàn)載荷步的循環(huán)加載，同時在每次循環(huán)中對單胞的疲勞狀態(tài)及失效狀態(tài)進行判斷，直至單胞達到最大疲勞壽命狀態(tài)。
2.3 結(jié)果展示
（1）單胞應力水平與疲勞壽命的關系
本文分別依據(jù)以上方法，進行了不同應力水平的單胞疲勞壽命預測，單胞應力水平與疲勞壽命的關系如圖2所示。
圖2 單胞疲勞壽命預測結(jié)果
（2）損失演化
其中，纖維束的損傷演化如圖3所示。
展開 1.用hypela2用戶子程序編寫的彈塑性材料本構(gòu)關系計算程序，其中
分別計算了彈性應變和塑性應變，該子程序的形參中有e(ngens)存放
總應變，可以在后處理中顯示(Total strain)。但是對于內(nèi)部的彈性
應變和塑性應變不知怎么樣更新讓它在后處理界面中顯示（elastic
strain和plastic strain），也就是說不知怎么樣在子程序中將計算
出來的彈性應變和塑性應變存放到系統(tǒng)中保存這兩個量的地方，可以
在每一步中更新；
2.用plotv也不行啊，
subroutine plotv(v, s, sp, etot, eplas, ecreep, t, m, nn, layer,ndi, nshear,
jpltcd)
etot：積分點處總塑性應變；
eplas：積分點處的塑性應變；
用了自己的本構(gòu)后這里的eplas就是0，如果用計算出的塑性應變直接
在這里以user defined variable來顯示的話則全部單元顯示的塑性
應變都是同一個值。
3.在c卷table3-3的單元后處理代碼(Element post codes)中，代碼
301代表總應變張量(Total strain tensor)，代碼321代表塑性應變
張量(plastic strain tensor)，401代表彈性應變張量(elastic
strain tensor)。
后處理代碼1-6說是應變分量(components of strain)，21-26代表
塑性應變分量(components o f plastic strain)，那么這些應變分量
跟前面的對應的應變張量有什么區(qū)別呢？
4.用hypela子程序?qū)崿F(xiàn)材料彈塑性本構(gòu)關系，根據(jù)小變形彈塑性來
推導的公式，計算了應力應變。
展開 利用MSC/PCL語言將該算法在MSC.Nastran上進行了程序?qū)崿F(xiàn)，開發(fā)出了優(yōu)化模塊的新版本。將每一步優(yōu)化迭代過程的初始值引入到用戶程序，經(jīng)過有限次迭代即可得到最優(yōu)設計結(jié)果，多個算例表明了程序的可靠性和精確性。
有需要的朋友可以去這里下載:
http://www.caenet.cn/paper/Paper.aspx?ID=381
程序?qū)崿F(xiàn)的相關專題、標簽、搜索
程序?qū)崿F(xiàn)的最新內(nèi)容
這個思想到今天看仍然非常高明，因為它在“物理真實性”和“程序可實現(xiàn)性”之間找到了很好的平衡。
第三，這篇文章并沒有急著把硬化寫得非常復雜。恰恰相反，作者采用了一個相對簡潔的 non-hardening 假設，認為材料宏觀上表現(xiàn)出來的許多“硬化感”，很大程度上其實來自取向變化和孿晶導致的晶格重排，而不只是每個滑移系的阻力不斷增加。
定制化服務：針對折疊屏手機、異形耳機等特殊產(chǎn)品，定制專用夾具與測試程序，實現(xiàn)折疊態(tài)、半開態(tài)、展開態(tài)多場景測試，覆蓋行業(yè)前沿產(chǎn)品檢測需求。
五、結(jié)語
在電子行業(yè)競爭白熱化的當下，抗摔性已成為產(chǎn)品核心競爭力，而非附加屬性。從試樣準備的嚴謹規(guī)范，到參數(shù)設置的科學精準，再到數(shù)據(jù)解讀的深度分析，跌落試驗機是企業(yè)把控產(chǎn)品質(zhì)量、降低售后率、守護品牌口碑的核心裝備。
程序員有責任實現(xiàn)能夠以某種方式正確反映這些愿景的代碼。
?InputField表示到達探測器的光場（等距采樣）。 遵循VirtualLab的模式概念，它由一組完全自相關的電磁模式組成，它們之間可以表現(xiàn)出不同的相干特性，如實地模擬物理場的相干特性。
?ParentLightPath是指包含相關探測器的光學系統(tǒng)。
多物理場耦合方便：DAMASK 的熱—機等耦合可通過 Marc 的 flux 類用戶子程序接入，實現(xiàn)塑性耗散生熱等效應的有機整合。
版本路線清晰：Abaqus 在 3.x 被官方棄用，Marc 則是 3.x 文檔中明確支持與維護的 FEM 選項，減少“接口過期”的不確定性與維護成本。
1 vumat與umat的區(qū)別
從程序?qū)崿F(xiàn)的角度，我們重點關注以下幾點區(qū)別：
? vumat不需要輸出一致性切線剛度矩陣
? vumat中應力應變存儲順序與umat不同
? vumat中存儲的應變值為張量應變值，而umat中為工程應變
? vumat的應力和狀態(tài)變量的更新方式不同，其分為old和new兩個數(shù)組
Abaqus/Explicit在啟動計算前，會進行數(shù)據(jù)檢查
需要注意的是，受壓不會造成粘結(jié)損傷（PS：若要考慮受壓粘結(jié)損傷，則需通過編寫子程序?qū)崿F(xiàn)）。</p><p class="ql-align-center"><strong>06</strong>二次開發(fā)接口</p><p>內(nèi)聚力單元與內(nèi)聚力接觸的二次開發(fā)接口也不同。
通過利用程序設計的力量，程序自動化實現(xiàn)了任務的自動執(zhí)行，從而釋放人力資源以從事更復雜和戰(zhàn)略性的活動。
程序自動化的好處：
1. 提高效率
2. 減少人工錯誤
3. 節(jié)省時間和成本
4. 可擴展性和可客制化
Moldex3D Studio Report Wizard 可以自動生成模擬結(jié)果的報告，可以自定義模板和結(jié)果項目。
通過利用程序設計的力量，程序自動化實現(xiàn)了任務的自動執(zhí)行，從而釋放人力資源以從事更復雜和戰(zhàn)略性的活動。
程序自動化的好處：
1. 提高效率
2. 減少人工錯誤
3. 節(jié)省時間和成本
4. 可擴展性和可客制化
Moldex3D Studio Report Wizard 可以自動生成模擬結(jié)果的報告，可以自定義模板和結(jié)果項目。
用機器學習破解新能源領域的“棄風”難題11個月前
其中：
數(shù)據(jù)清理部分，項目團隊基于軟件的數(shù)據(jù)處理模塊，開發(fā)自編程清洗程序，實現(xiàn)數(shù)據(jù)高效批量化處理。同時考慮風機工作特征和數(shù)據(jù)異常原因，將建模過程結(jié)合業(yè)務背景。
敏感性分析部分篩選掉“不重要”的自變量，提高建模精度。
模型訓練部分，先根據(jù)風電特征將數(shù)據(jù)分區(qū)，然后分別用不同的訓練算法做模型訓練。
模型對比階段，用測試集對基于不同算法訓練得到的模型進行精度對比。
圖2 材料屬性構(gòu)建
3、 激光熱源子程序開發(fā)
（1） 熱源特性：采用高斯分布模擬圓形激光束，功率密度函數(shù)為：
其中，P 為激光功率，r0為光斑半徑，r 為徑向坐標
（2） 子程序?qū)崿F(xiàn)：基于ABAQUS的用戶子程序接口（如DFLUX或HETVAL），編寫 Fortran/Python 程序生成動態(tài)加載的圓形激光熱源，通過時間 - 空間函數(shù)控制熱源移動軌跡
