模型：2D軸對稱T-CPTU模型，先貫入，再給熱源，實現加熱-自然散熱（取消熱源）的熱傳導過程。土體：MCC本構，探頭：剛體。涉及大變形以及重啟動操作，因此有兩個模型：PENE、TEMP，TEMP是在PENE的基礎上建立的。

First model：PENE

Part&Property&Assembly：為了后續進行傳熱，探頭要采用變形體，但是顯然不對的，這里我做了兩個改變：（1）把彈性模量拉到很大；（2）給探頭整體一個剛體約束，避免探頭本身的變形。土體采用多孔彈性+黏土塑性，通過設置參數將其等價于MCC，需要注意的是：在設置滲透系數的時候，后面的孔隙比代表的意思是滲透系數隨孔隙比的變化，如果設置滲透系數為常參數的話，孔隙比可以是任意值，在后續的預定義場施加初始孔隙比即可。探頭或者土體的熱物性參數都要給：熱導率、比熱、熱膨脹系數。主包給土體設置了一個damping參數為了增加收斂性，Alpha-0.05，Beta-0.0005（不保證啊不保證，也不知道有沒有用）。

Step：2步====分析步均采用耦合溫度-位移分析。(1)geo，地應力平衡，transient /1s ，打開大變形，增量步選擇automatic 采用非對稱求解器；（2）pene，貫入分析步（要考的，記清楚），前兩個分析步均未采用automatic stabilization，但是一定要打開大變形選項防止網格過度扭曲。最重要的操作：重啟動。Step界面——Output——restart requests——在geo/pene分析步勾選frequency以及overlay。即每個增量步讀取一次數據以及后續在該步驟可以重新啟動計算（是這個意思嗎？不知道啊再找找資料吧主包）。

Interaction：建立了7個接觸，探頭的各個分區與土體左邊界之間。探頭分為金屬區域和特氟龍隔熱區，兩個區域的接觸屬性不同，主要是比熱、熱導率、熱擴散系數的區別。另外就是前面提到的剛體約束rigid body（給探頭的）。

Load：在PENE模型的load僅有一個上覆荷載P=50kpa，三個邊界條件：土體底部位移全固定，右邊界水平位移=0，探頭設置一個參考點，初始固定全部位移，pene分析步設置幅值勻速貫入。由于預定義場要在initial分析步創建，因此在PENE模型中的Predefined Field要定義：初始溫度場、孔隙比場、應力場。土體：初始溫度+孔隙比+應力，探頭：初始溫度，一共4個預定義場。

Mesh：網格劃分的很丑好在能跑。不建議學我的。

Job：終于來到了作業，建立一個test-pene的job文件（要考的，記清楚），CPU拉滿開始算。結束可以收獲test-pene.odb。OK啊朋友們，PENE模型結束掉了。