abaqus子程式技術文件
我們將有限元軟體（以abaqus為例）分割來看，分為：前處理，單元部分，求解器。在這三個模組當中，每個主流商業軟體手冊中都會註明各個單元的理論採用了哪種理**式，但是只是提一下用什麼方法修正，很多沒有介紹具體實現的公式。這也會讓許多使用人員感覺在用黑盒子，對於分析問題的本質會造成很大的阻礙。下面我們主要講解uel使用者子程式的使用方法。
uel使用者子程式開發步驟
使用者子程式主要是用來實現特定的材料本構模型和單元演算法的有效介面，便於商用軟體和求解器之間的互動迭代。
程式語言主要採用fortran，當然也可以修改子程式介面，使用matlab或者c++等語言進行實現。
1. uel的關鍵輸入輸出引數
|關鍵引數名稱| 型別 |描述|
rhs輸出 第一列是殘餘力向量；第二列是節點外部載荷增量|
amatrx 輸出 剛度/質量/阻尼矩陣，由lflags決定
svars 輸入/輸出 狀態相關變數，大小為nsvars
energy 輸入/輸出 單元能量
ndofel 輸入 單元自由度數量
nrhs 輸入 載荷向量個數rhs的列數，非線性演算法一般為1
nsvars 輸入 單元狀態變數的數量
props 輸入 使用者定義的單元屬性資料，浮點陣列
nprops 輸入 使用者定義的單元屬性資料的數量
coords 輸入 單元節點的座標
mcrd 輸入 任意節點使用者所需定義的最大座標數的最大值
nnode 輸入 單元節點數
u, du 輸入 當前增量步單元節點的位移，位移增量
v 輸入 相關變數的時間變化率
a 輸入 相關變數的加速度
jtype 輸入 定義單元型別的整型值
time 輸入 當前分析步時間和時間總量
dtime 輸入 時間增量
kstep 輸入 當前分析步的數
kinc 輸入 當前增量步數
jelem 輸入 使用者指派的單元號
params 輸入 和求解過程相關的引數
ndload 輸入 作用在當前單元上的熱力學力
jdltyp 輸入 用來定義單元分布載荷型別的整型陣列
adlmag 輸入 當前增量步結束時分布載荷相關陣列
predef 輸入 預定義場變數陣列
npredf 輸入 預定義場變數的數目
lflags 輸入 當前求解步的標識陣列
mlvarx 輸入 當使用多個位移或右手邊向量時的量綱引數
ddlmag 輸入 關於分布載荷大小的增量
mdload 輸入 單元上熱力學力的總數量
pnewdt 輸入 建議的新時間增量與當前使用的時間增量之比
jprops 輸入 單元中的整型屬性值
njprop 輸入 整型屬性值的數量
period 輸入 當前分析步的時間段
2. 基於fortran的uel開發步驟
在inp中定義必要的單元引數，網格資訊，分析步，輸出變數等資訊
子程式字尾為.for格式
編譯採用abaqus6.14+intel fortran2013+visual studio2013的整合開發環境
其餘版本也可匹配整合
在環境配置完成之後，開啟abaqus command，輸入命令abaqus make library=***.for，即開始編譯，編譯過程中的警告和錯誤都會列印在命令提示框內。（這一步為非必要步驟）
執行一種方式是通過command輸入命令列，另一種方式則是gui方式在建立job時進行選擇。
除錯（debug）
採用在程式中加入特定程式段的方式
3. 小結
uel開發需要打下堅實的力學基礎，熟悉有限元分析過程，掌握一定的演算法實現和程式設計能力。
Abaqus熱源子程式Dflux教程
有限元法模擬焊接過程，要根據熱源分布數學模型來描述焊接過程的輸入熱源，常用的有兩種 熔池深度小，可近似認為熱量從表面傳入，面熱源 假設熱源沿x方向以速度v移動，熱源表示式則為 q為熱源有效功率 熔池深度大，考慮熱源在深度方向的變化，體熱源 假設熱源沿x方向以速度v移動，熱源表示式則為 其中，f1，f...
微信小程式開發技術文件
一 目錄結構 乙個小程式主體部分由三個檔案組成，必須放在專案的根目錄，如下 2.乙個小程式頁面由四個檔案組成，分別是 注意 為了方便開發者減少配置項，描述頁面的四個檔案必須具有相同的路徑與檔名。二 配置小程式 pages index index window 三 小程式框架 註冊頁面 對於小程式中的...
小議程式設計師編寫技術文件
一提到寫文件，可能很多程式設計師可能會不屑一顧，但是，無論處於規範開發流程，還是就於逃避嫌責的目的，能夠將自己所從事的工作用文件描述記錄下來，還是一件很有成就感的事情，拋開其功用不談，就個人的成長程序看，也是乙個循序漸進式的好習慣，還是值得大家稍微關注一下的。昨天在和同事的一次交流過程中，就自己編寫...