非線性用戶單元的開發(fā)

非線性用戶單元的開發(fā)第1頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系ABAQUS的用戶子程序第2頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系用戶子程序用戶子程序為分析提供了一個非常強大、靈活的工具。用戶子程序：可以提高幾個ABAQUS選項的功能，僅用數(shù)據(jù)行的方式使用這些選項可能太過嚴格；用FORTRAN語言寫成；在某些情況下可以利用ABAQUS的功能程序。第3頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系ABAQUS中的用戶子程序：CREEP－定義時間相關的粘塑性行為（蠕變和膨脹）的用戶子程序DFLOW－在壓實分析中定義非均勻孔隙流體速度的用戶子程序DFLUX－在熱傳導或質量擴散分析中定義非均勻分布流的用戶子程序DISP－定義邊界條件的用戶子程序第4頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系DLOAD－定義非均勻分布載荷的用戶子程序FILM－在熱傳導分析中定義非均勻膜系數(shù)和相關的下沉溫度的用戶子程序FLOW－在壓實分析中定義非均勻滲流系數(shù)和相關的下沉孔隙壓力的用戶子程序FRIC－定義接觸面摩擦行為的用戶子程序GAPCON－在溫度－位移耦合或者純熱傳導分析中定義接觸面或節(jié)點之間的熱傳導系數(shù)的用戶子程序第5頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系GAPELECTR－在熱－電耦合分析中定義表面間導電系數(shù)的用戶子程序HARDINI－定義初始等效塑性應變和初始背應力張量的用戶子程序HETVAL－在熱傳導分析中定義內部熱產(chǎn)生的用戶子程序MPC－定義多點約束的用戶子程序ORIENT－為運動耦合約束定義局部方向或者定義局部材料取向的用戶子程序第6頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系RSURFU－定義剛性表面的用戶子程序SDVINI －定義依賴于解的狀態(tài)變量場初始值的用戶子程序SIGINI－定義初應力場的用戶子程序UCORR－定義隨機響應載荷的交叉 －相關屬性的用戶子程序UEL－定義單元的用戶子程序第7頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系非線性用戶單元的開發(fā)動機定義用戶單元UEL接口例：具有非線性截面行為的平面梁單元例：新開發(fā)的兩種用戶單元第8頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系概要ABAQUS/Standard提供接口，允許用戶實現(xiàn)線性和非線性有限單元。非線性有限單元在用戶子程序UEL中實現(xiàn)。接口允許定義任意復雜的單元。如果編寫適當，用戶單元可以在ABAQUS/Standard的大部分程序中應用?？梢栽谝粋€UEL子程序中實現(xiàn)多個用戶單元，這些單元可以同時使用。這里只討論非線性用戶單元的實現(xiàn)，并采用實例說明。第9頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系動機ABAQUS/Standard是一個功能強大的分析工具，擁有一個大單元庫，有500多種單元－可以進行絕大部分復雜結構問題的分析。然而，在有些情況下需要開發(fā)用戶單元來補充ABAQUS的單元庫：所需的單元幾何信息與ABAQUS已有的單元不同模擬與結構行為耦合的非結構物理過程施加依賴于解的載荷模擬主動控制機理第10頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系在分析程序如ABAQUS中開發(fā)用戶單元而不是編寫一個完全的分析程序的優(yōu)點是很明顯的：ABAQUS提供大量的結構單元、分析程序和建模工具。ABAQUS提供前處理和后處理。許多第三方軟件提供帶有ABAQUS接口的前后處理器。維護和移植子程序比維護和移植一個完全的有限元程序容易的多。第11頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系定義用戶單元用戶單元的主要特征在編寫UEL子程序前，必須先定義單元的以下主要特征：單元的節(jié)點數(shù)節(jié)點的坐標數(shù)每一個節(jié)點處的自由度數(shù)第12頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系其它重要的單元性質另外，必須定義下面的單元性質：需要在UEL外定義的單元性質數(shù)目每一個單元需要存儲的依賴于解的狀態(tài)變量數(shù)目（SDVs）可能會作用于單元上的（分布）載荷類型數(shù)目這些數(shù)目不必立即定義：它們可以很容易地在完成基本的UEL子程序后加入。第13頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系定義用戶單元行為在一般的分析步中單元對模型的主要貢獻是提供作用于節(jié)點上的“力”FN－依賴于節(jié)點的自由度uN的值。如果自由度是位移，那么相關的“力”為節(jié)點力。類似地，旋轉對應于轉矩，溫度對應于熱流。在非線性用戶單元中“力”常常依賴于自由度的增量uN和內部狀態(tài)變量H。狀態(tài)變量必須在用戶子程序中進行更新。第14頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系在一般分析步中求解（非線性）方程組需要你定義單元的剛度矩陣：剛度矩陣應包含F(xiàn)N對uN所有的直接和間接的依賴性，例如在一般分析步中剛度矩陣定義的越精確，收斂性就越好。剛度矩陣確定線性擾動分析步的解，所以其定義必須精確。剛度矩陣可以為對稱或非對稱。第15頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系用戶單元的復雜程度可以變化很大。在由常規(guī)單元構成的分析中可以采用簡單的用戶單元來達到“控制”和“反饋”的功能。復雜的非線性結構單元的編制常常需要很大的工作量。如果單元由非線性材料構成，你最好用單獨的子程序描述材料行為。如果材料模型是在子程序UMAT中實現(xiàn)的，那么可以在UEL中包含對UMAT的調用。前面對UMAT中積分的討論同樣適用于UEL中的材料模型。第16頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系UEL接口ABAQUS選項用戶單元是用*USERELEMENT選項定義的。在輸入文件中這個選項必須出現(xiàn)在引用用戶單元的*ELEMENT選項的前面。輸入文件中用于UEL的語法如下：

*USERELEMENT,TYPE=Un,NODES=,COORDINATES=, PROPERTIES=,IPROPERTIES=,VARIABLES=,UNSYMM Datalines(s) *ELEMENT,TYPE=Un,ELSET=UEL Dataline(s) *UELPROPERTY,ELSET=UEL Dataline(s) *USERSUBROUTINE,(INPUT=file_name)*USERSUBROUTINE選項在ABAQUS6.1版本中將不再用。第17頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系UEL接口用戶子程序UEL的接口如下：

SUBROUTINEUEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS,1PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,2KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF,3NPREDF,LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,4PERIOD)CINCLUDE'ABA_PARAM.INC'CDIMENSIONRHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),1SVARS(NSVARS),ENERGY(8),PROPS(*),COORDS(MCRD,NNODE),2U(NDOFEL),DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2),3PARAMS(3),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*),4DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAGS(*),JPROPS(*)“INCLUDE”聲明語句為浮點變量設置合適的精度（大部分機器上為REAL*8）。第18頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系UEL變量ABAQUS為用戶單元UEL提供下面的變量：坐標、位移、增量位移，對于動態(tài)，還有速度和加速度增量開始時刻的狀態(tài)變量SDVs總時間和增量時間、溫度以及用戶定義的場變量用戶單元性質載荷類型及當前總載荷值和載荷增量值單元類型和用戶定義的單元號分析類型標志，對于動態(tài)，積分算子值當前所在的分析步數(shù)和增量步數(shù)第19頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系必須定義下面的變量：右手端矢量（節(jié)點流或力）剛度矩陣依賴于解的狀態(tài)變量可以定義下面的變量：單元的能量（應變能、塑性耗散、動能等）建議的新時間增量所有參數(shù)的詳細描述請參見ABAQUS/Standard用戶手冊的第6.2.18節(jié)。第20頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系編寫和測試UEL編寫ABAQUS用戶子程序的基本規(guī)則：遵從FORTRAN77或C的語法。確保所有的變量都定義和初始化過。為狀態(tài)變量分配足夠的存儲空間。ABAQUS5.8-10版本要求FORTRAN編譯器的版本為5.0；從ABAQUS5.8-14開始，要求FORTRAN編譯器的版本為6.0。第21頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系復雜的UELs可能會存在很多潛在的問題，調試一個不收斂的包含UEL的大模型會帶來不必要的困難。用只有一個單元的輸入文件驗證UEL：首先，用一般分析步進行測試－將所有的解變量都賦值，用來驗證總流（力）。其次，用線性擾動步進行測試－將所有的載荷都賦值，用來驗證單元的剛度矩陣。最后，用一般分析步進行測試－將所有的載荷都賦值，用來驗證剛度矩陣和流矢量的匹配。逐漸增大測試問題的復雜程度。如果可能，與標準ABAQUS單元的結果進行比較。第22頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系例1：具有非線性截面行為的

平面梁單元目標分析一個平面混凝土框架結構。這個框架所受載荷很大，混凝土內已有明顯的非線性，但位移仍然很小，可以不考慮幾何非線性。建立一個模型：直接用軸向力和彎矩描述非線性截面行為。這與選項*BEAMSECTION,SECTION=NONLINEARGENERAL類似，但是允許軸向和彎曲項之間的耦合。忽略橫向剪切變形。第23頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系編程要求單元采用數(shù)值積分；因此，需要在UEL中定義下面的量：單元[B]矩陣，用于聯(lián)系軸向應變、曲率與單元位移{ue}:本構律矩陣[D]，用于聯(lián)系軸向力F、彎矩M與軸向應變、曲率：第24頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系單元剛度矩陣：單元力矢量：采用數(shù)值積分：其中n為積分點數(shù)，li為與積分點i相關聯(lián)的長度。第25頁，共30頁，2023年，2月20日，星期四ABAQUS6.1版本發(fā)布會暨99中國地區(qū)用戶會議 清華大學工程力學系評論這個UEL所采用的描述和簡單的B23單元