Abaqus子程序之UTEMP
1.簡介
用戶子程序 UTEMP ：
允許您規(guī)定模型節(jié)點(diǎn)的溫度；
每當(dāng)用戶定義的溫度場定義下列出的節(jié)點(diǎn)需要當(dāng)前溫度值時(shí)，將被調(diào)用；
忽略為用戶子程序之外的關(guān)聯(lián)溫度場定義提供的任何溫度；和
可用于修改從結(jié)果文件中讀取的任何溫度。
2.本文討論：
用戶子程序接口
要定義的變量
傳遞信息的變量
3.調(diào)用方法
產(chǎn)品Abaqus/Standard
4.用戶子程序接口
SUBROUTINE UTEMP(TEMP,NSECPT,KSTEP,KINC,TIME,NODE,COORDS)
C? ? ?
INCLUDE 'ABA_PARAM.INC'C
DIMENSION TEMP(NSECPT), TIME(2), COORDS(3)
C user coding to define TEMP
RETURN
END
5.要定義的變量
TEMP(NSECPT)
節(jié)點(diǎn)號(hào) NODE 處的溫度值數(shù)組。如果節(jié)點(diǎn)未連接到梁或殼單元，則僅需返回一個(gè)溫度值 (NSECPT=1)。否則，要返回的溫度數(shù)取決于為梁或殼截面選擇的溫度模式和場變量輸入。以下情況是可能出現(xiàn)的：
梁截面的溫度和場變量作為梁截面描述中顯示的點(diǎn)處的值給出。所需值的數(shù)量 NSECPT 由指定的特定截面類型確定，如梁截面庫中所述。
溫度和場變量作為殼截面每一層的 n 個(gè)等距點(diǎn)的值給出。所需值的數(shù)量 NSECPT 等于 n。
梁截面的溫度和場變量以截面原點(diǎn)處的值以及相對(duì)于 2 方向的梯度給出，對(duì)于三維梁，截面的 1 方向；或殼截面的溫度和場變量作為參考表面的值以及相對(duì)于厚度的梯度給出。所需值的數(shù)量 NSECPT 對(duì)于三維梁為 3，對(duì)于二維梁為 2，對(duì)于殼為 2。如梁單元和殼單元中所述，首先給出中間曲面值，然后是第一個(gè)和（如果需要的話）第二個(gè)梯度。
還可以請(qǐng)求以其他兩種方式之一設(shè)置溫度：從先前生成的結(jié)果文件或通過直接數(shù)據(jù)輸入。當(dāng)數(shù)組 TEMP 被傳遞到用戶子程序 UTEMP 時(shí)，它將僅包含從結(jié)果文件中獲得的溫度。您可以在此例程中修改這些值。任何作為直接數(shù)據(jù)輸入給出的值都將被忽略。
6.傳遞信息的變量
NSECPT
模型中任何節(jié)點(diǎn)所需的最大截面值數(shù)。
KSTEP
分析步號(hào)
KINC
增量步號(hào)
TIME(1)
步驟時(shí)間的當(dāng)前值。
TIME(2)
總時(shí)間的當(dāng)前值。
NODE
節(jié)點(diǎn)號(hào)
COORDS
包含該點(diǎn)當(dāng)前坐標(biāo)的數(shù)組。
7.案例
這里為了便于學(xué)習(xí)理解，給出一個(gè)小案例吧，首先UTEMP子程序自定義了節(jié)點(diǎn)溫度，接著使用UFIELD子程序定義了與溫度有關(guān)的三個(gè)場變量。
結(jié)語
前段時(shí)間太忙，所以更新的很慢，看到有同學(xué)私信我寫個(gè)Vufield的教程，實(shí)際上V開頭的子程序是提供給顯式動(dòng)力分析求解器用的，也就是說ufield和Vufield差異性很小，如果有問題歡迎私信一起交流學(xué)習(xí)（有空的話一定會(huì)回復(fù)的）。需要案例的公眾號(hào)回復(fù)UTEMP 即可。
