!!!!圓筒焊接熱應力及殘余應力分析命令流
/prep7
/units,si
!///////////////////////定義材料屬性
ET,1,SOLID70
ET,2,PLANE55
ET,3,SURF152
KEYOPT,3,3,0
KEYOPT,3,4,1
KEYOPT,3,5,0
KEYOPT,3,6,0
KEYOPT,3,7,0
KEYOPT,3,8,5
KEYOPT,3,9,0
TOFFST,273    !!!!指定初始零度與絕對零度偏差值
!///////////////////////定義材料屬性
MPTEMP,1,0,100,200,300,400,500,600,1200,3000  !!!!!溫度
MPDATA,KXX,1,1,52,50.7,48.6,46.1,42.3,38.9,35.6,30,30  !!!!導熱系數
MPTEMP,1,0,100,200,400,500,1000,1480,1485,1500,1505,3000
MPDATA,C,1,1,450,469,481,536,569,700,700,8700,8700,700,700  !!!!定壓比熱
MPTEMP,1,0,100,200,400,500,1600,1750,2500
MPDATA,DENS,1,1,7820,7820,7800,7800,7800,7800,7800,7800    !!!!密度
MPTEMP,1,20,400,800,1200,1500,2500
MPDATA,HF,1,1,100,480,1050,2000,3000,3500    !!!!對流系數

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,1,20

!!!!!!!!material num of surface elememt
MPTEMP,1,0,100,200,300,400,500,600,1200,3000
MPDATA,KXX,3,1,52,50.7,48.6,46.1,42.3,38.9,35.6,30,30
MPTEMP,1,0,100,200,400,500,1000,1480,1485,1500,1505,3000
MPDATA,C,3,1,450,469,481,536,569,700,700,8700,8700,700,700
MPTEMP,1,0,100,200,400,500,1600,1750,2500
MPDATA,DENS,3,1,7820,7820,7800,7800,7800,7800,7800,7800
MPTEMP,1,20,400,800,1200,1500,2500
MPDATA,HF,3,1,100,480,1050,2000,3000,3500
!!!!!!!!!!!!!!定義材料屬性
!!!!!!!!!!!!!!!!!!!!!!!Create the model
CYL4,,,0.05,,0.047   !!!!在工作面上生成一個圓面或圓柱體
TYPE,2      !!!!!!單元類型
MAT,1       !!!!!!單元材料屬性指示器
REAL,       !!!!!!單元實常數屬性指示器
ESYS,0      !!!!!!單元坐標系統屬性指示器
SECNUM,     !!!!!!單元剖面屬性指示器
LESIZE,ALL,,,25,,1,,,1,      !!!!!!對所選線設置網格單元大小，25為分割角度數
MSHKEY,0    !!!!!!指定網格劃分方式，0為自由網格劃分
CM,_Y,AREA  !!!!!!根據所選實體生成一個元件
ASEL,,,,1
CM,_Y1,AREA
CHKMSH,'AREA'      !!!!!檢查劃分實體類型
CMSEL,S,_Y         !!!!!選擇一組元件或部件子集，s是選擇一個新設置
AMESH,_Y1          !!!!!在面上生成節點或點單元
CMDELE,_Y          !!!!!刪除一個現存的部件
CMDELE,_Y1
CMDELE,_Y2
TYPE,1
EXTOPT,ESIZE,25,0,    !!!!!面單元生成體單元的相關選項  
EXTOPT,ACLEAR,1       !!!!!ACLEAR生成體時面單元網格是否清除
EXTOPT,ATTR,0,0,0     !!!!!ATTR將面單元上的特殊屬性傳給生成的體單元
MAT,1
REAL,_Z4
ESYS,0
FLST,2,1,5,ORDE,1      !!GUI拾取命令，所選操作所需的特殊參數，Flst命令是GUI操作的拾取命令，總是與FITEM命令一起用，舉例明： 
                       !!!!!!!!FLST,2,4,4,ORDE,2 
                       !!第一個2表示拾取項作為后面命令的第一個條件，第一個4 表示拾取4項 
                       !!第三個4 表示拾取直線號 最后一個2 表示有2項FITEM 
                       !!!!!!FITEM,2,1 
                       !!!!!!FITEM,2,-4 !負號表示與上面同類，即拾取1，2，3，4四條線 
                       !!!!!!LCCAT,P51X !拾取的線作為ＬＣＣＡＴ的第一個條件

FITEM,2,1
VEXT,P51X,,,0,0,0.1,,,,    !!!!通過給定偏移量由面生成體

!!!!!!!!!!!!!!!!!!創建表面效應單元
FLST,5,4,5,ORDE,2
FITEM,5,3
FITEM,5,-6
ASEL,S,,,P51X
NSLA,S,1                !!!!!!選擇與選中面相關的節點,s選一套新節點,1選所有和面相聯系的節點（如面內線，關鍵點處的節點）
TYPE,3
MAT,2
REAL,
ESYS,0
SECNUM,
TSHAP,ARC              !!!!!!為目標單元指定2D或3D幾何表面，ARC-2D順時針圓弧
!*
FLST,5,2600,1,ORDE,8
FITEM,5,1
FITEM,5,-100
FITEM,5,201
FITEM,5,-300
FITEM,5,401
FITEM,5,-496
FITEM,5,593
FITEM,5,-2896
CM,_Y,node
NSEL,,,,P51X
CM,_Y1,node
CMSEL,,_Y
!*
CMSEL,,_Y1
ESURF,0              !!!!!!在存在的已選單元的自由面上生成重疊單元
CMSEL,,_Y
CMDELE,_Y
CMDELE,_Y1
ALLSEL,ALL
!!!!!!!!!!!!!!!!創建表面效應單元
!!!!!!!!




NSEL,S,LOC,Z,0,0.004
ESLN,S,1        !!!!!!選擇已選點上的單元
EMODIF,ALL,MAT,3,      !!!!!!改變選中的單元類型為所需要的類型
ALLSEL,ALL
!!!!!!!!!!!!!!求解
/SOLU
ANTYPE,4
TRNOPT,FULL      !!!!!!指定瞬態分析選項 
LUMPM,0          !!!!!!指定一個集中質量矩陣公式，0 使用與單元相關的質量矩陣公式
TIME,0.25        !!!!!!指定荷載步結束時間
AUTOTS,1         !!!!!!是否使用自動時間步長,1: 由程序選擇（當solcontrol為on且不發生autots命令時在 .log文件中紀錄“1”
                 !!!!!!注意：當使用自動時間步長時，也會使用步長預測器和二分步長
DELTIM,0.04,0.03,0.05,1      !!!!!!定義時間步長
KBC,1                        !!!!!!指定載荷增加方式,1為階躍荷載
TSRES,ERASE                  !!!!!!時間步模式改變時定義關鍵時間數組
TIMINT,1                     !!!!!!打開瞬態效應
TINTP,0.005,,,1,0.5,0        !!!!!!指定瞬態積分參數
OUTRES,ALL,LAST,             !!!!!!規定寫入數據庫的求解信息
FLST,2,5200,1,ORDE,2
FITEM,2,1
FITEM,2,-5200
IC,P51X,TEMP,20,             !!!!!!在節點上指定初始條件






*DEL,_FNCNAME
*DEL,_FNCMTID
*SET,_FNCNAME,'cweld'
!/INPUT,cweld.func
*DIM,%_FNCNAME%,TABLE,6,27,1     !!!!!!定義數組參數及維數
!
!Begin of equation:91351697*2.718^(-99649*(({x}-0.05*Sin!(0))^2+({Y}-0.05*cos(0))^2+{Z}^2))
%_FNCNAME%(0,0,1)=0.0,-999
%_FNCNAME%(2,0,1)=0.0
%_FNCNAME%(3,0,1)=0.0
%_FNCNAME%(4,0,1)=0.0
%_FNCNAME%(5,0,1)=0.0
%_FNCNAME%(6,0,1)=0.0
%_FNCNAME%(0,1,1)=1.0,-1,0,0,0,0,0
%_FNCNAME%(0,2,1)=0.0,-2,0,1,0,0,-1
%_FNCNAME%(0,3,1)=0.0,-3,0,1,-1,2,-2
%_FNCNAME%(0,4,1)=0.0,-1,0,0,0,0,-3
%_FNCNAME%(0,5,1)=0.0,-2,0,1,-3,3,-1
%_FNCNAME%(0,6,1)=0.0,-1,9,0,0,0,0
%_FNCNAME%(0,7,1)=0.0,-3,0,0.05,0,0,-1
%_FNCNAME%(0,8,1)=0.0,-4,0,1,-3,3,-1
%_FNCNAME%(0,9,1)=0.0,-1,0,1,2,2,-4
%_FNCNAME%(0,10,1)=0.0,-3,0,2,0,0,-1
%_FNCNAME%(0,11,1)=0.0,-4,0,1,-1,17,-3
%_FNCNAME%(0,12,1)=0.0,-1,10,314.2-314.2,0,0,0
%_FNCNAME%(0,13,1)=0.0,-3,0,0.05,0,0,-1
%_FNCNAME%(0,14,1)=0.0,-5,0,1,-3,3,-1
%_FNCNAME%(0,15,1)=0.0,-1,0,1,3,2,-5
%_FNCNAME%(0,16,1)=0.0,-3,0,2,0,0,-1
%_FNCNAME%(0,17,1)=0.0,-5,0,1,-1,17,-3
%_FNCNAME%(0,18,1)=0.0,-1,0,1,-4,1,-5
%_FNCNAME%(0,19,1)=0.0,-3,0,2,0,0,4
%_FNCNAME%(0,20,1)=0.0,-4,0,1,4,17,-3
%_FNCNAME%(0,21,1)=0.0,-3,0,1,-1,1,-4
%_FNCNAME%(0,22,1)=0.0,-1,0,1,-2,3,-3
%_FNCNAME%(0,23,1)=0.0,-2,0,2.718,0,0,-1
%_FNCNAME%(0,24,1)=0.0,-3,0,1,-2,17,-1
%_FNCNAME%(0,25,1)=0.0,-1,0,91351697,0,0,-3
%_FNCNAME%(0,26,1)=0.0,-2,0,1,-1,3,-3
%_FNCNAME%(0,27,1)=0.0,99,0,1,-2,0,0
!END of equation:91351697*2.718^(-99649*(({x}-0.05*Sin!(0))^2+({Y}-0.05*cos(0))^2+{Z}^2))
FLST,2,4,5,ORDE,2
FITEM,2,3
FITEM,2,-6






/GO
SFA,P51X,,HFLUX,%CWELD%          !!!!!!在指定面上加荷載
!!!!!!!!!!!!!!!!!!!!!!!
FLST,2,4,5,ORDE,2
FITEM,2,7
FITEM,2,-10
/GO
!*
SFA,P51X,1,CONV,-1,20
!///////////////
ALLSEL,ALL





ESEL,S,MAT,,2


FLST,5,2500,2,ORDE,2
FITEM,5,2501
FITEM,5,-5000
CM,_Y,ELEM
ESEL,,,,P51X
CM,_Y1,ELEM
CMSEL,S,_Y
CMDELE._Y
!*



ALLSEL,ALL
/GO
!*
SFE,_Y1,_Z4,CONV,,-1         !!!!!!分布力作用單元上
SFE,_Y1,_Z4,CONV,2,20
CMDELE,_Y1
ALLSEL,ALL
solve
finish





*do,i,1,100
/SOLU
ANTYPE,,REST
TIME,0.25*(i+1)
AUTOTS,1
DELTIM,0.04,0.03,0.05,
KBC,1
TSRES,ERASE
*DEL,cweld
*DEL,_FNCNAME
*DEL,_FNCMTID
*SET,_ FNCNAME,'cweld'
!/INPUT,cweld.func
*DIM,%_FNCNAME%,TABLE,6,27,1
!
!Begin of equation:91351697*2.718^(-99649*(({X}-0.05*sin !(3.6*i))^2+({Y}-0.05*cos (3.6*i))^2+{Z}^2))
%_FNCNAME%(0,0,1)=0.0,-999
%_FNCNAME%(2,0,1)=0.0
%_FNCNAME%(3,0,1)=0.0
%_FNCNAME%(4,0,1)=0.0
%_FNCNAME%(5,0,1)=0.0
%_FNCNAME%(6,0,1)=0.0
%_FNCNAME%(0,1,1)=1.0,-1,0,0,0,0,0
%_FNCNAME%(0,2,1)=0.0,-2,0,1,0,0,-1
%_FNCNAME%(0,3,1)=0.0,-3,0,1,-1,2,-2
%_FNCNAME%(0,4,1)=0.0,-1,0,99649,0,0,-3
%_FNCNAME%(0,5,1)=0.0,-2,0,1,-3,3,-1
%_FNCNAME%(0,6,1)=0.0,-1,9,3.6/180*3.14159265*i,0,0,0
%_FNCNAME%(0,7,1)=0.0,-3,0,0.05,0,0,-1
%_FNCNAME%(0,8,1)=0.0,-4,0,1,-3,3,-1
%_FNCNAME%(0,9,1)=0.0,-1,0,1,2,2,-4
%_FNCNAME%(0,10,1)=0.0,-3,0,2,0,0,-1
%_FNCNAME%(0,11,1)=0.0,-4,0,1,-1,17,-3
%_FNCNAME%(0,12,1)=0.0,-1,10,3.6/180*3.14159265*i,0,0,0
%_FNCNAME%(0,13,1)=0.0,-3,0,0.05,0,0,-1
%_FNCNAME%(0,14,1)=0.0,-5,0,1,-3,3,-1
%_FNCNAME%(0,15,1)=0.0,-1,0,1,3,2,-5
%_FNCNAME%(0,16,1)=0.0,-3,0,2,0,0,-1
%_FNCNAME%(0,17,1)=0.0,-5,0,1,-1,17,-3
%_FNCNAME%(0,18,1)=0.0,-1,0,1,-4,1,-5
%_FNCNAME%(0,19,1)=0.0,-3,0,2,0,0,4
%_FNCNAME%(0,20,1)=0.0,-4,0,1,4,17,-3
%_FNCNAME%(0,21,1)=0.0,-3,0,1,-1,1,-4
%_FNCNAME%(0,22,1)=0.0,-1,0,1,-2,3,-3
%_FNCNAME%(0,23,1)=0.0,-2,0,2.718,0,0,-1
%_FNCNAME%(0,24,1)=0.0,-3,0,1,-2,17,-1
%_FNCNAME%(0,25,1)=0.0,-1,0,91351697,0,0,-3
%_FNCNAME%(0,26,1)=0.0,-2,0,1,-1,3,-3
%_FNCNAME%(0,27,1)=0.0,99,0,1,-2,0,0
!END of equation:91351697*2.718^(-99649*(({X}-0.05*sin 
!(314.2))^2+({Y}-0.05*cos (314.2))^2+{Z}^2))
FLST,2,4,5,ORDE,2
FITEM,2,3
FITEM,2,-6
/GO
SFA,P51X,,HFLUX, %CWELD%
solve
finish
*enddo
finish
!///////////
*do,i,101,180
/SOLU
ANTYPE,,REST
TIME,0.25*(i+1)
AUTOTS,1
DELTIM,0.015,0.012,0.02,
KBC,1
TSRES,ERASE
FLST,2,4,5,ORDE,2
FITEM,2,3
FITEM,2,-6
SFADELE,P51X,1,HFLUX
SOLVE
FINISH
*enddo
finish





!!!!!!!!!!!!!!!!!!!!!!!!應力場分析






/PREP7
ETCHG,TTS       !!!!!!改變單元類型到相應類型TTS熱到結構
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,3,,5.00e9
MPDATA,PRXY,3,,0.45      !!!!!!??????
TB,BISO,3,1,2,           !!!!!!定義非線性材料特性表Biso: 雙線性等向強化,材料號,數據的溫度數,對某一給定溫度數據的點數
TBTEMP,0                 !!!!!!為材料表定義溫度值
TBDATA,,306e6,0.26e9,,,, !!!!!!給當前數據表定義數據,所要輸入數據在數據表中的初始位置，缺省為上一次的位置加1,
!!!!!!!!!!每重新發生一次tbtemp命令上一次位置重設為1，（發生tb后第一次用空閑此項，則c1賦給第一個常數）

MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,3,REFT,,,           !!!!!!求解過程中修改材料特性,3材料編號REFT,參考溫度，必需定義為常量
MPDELE,ALPX,3            !!!!!!刪除線性材料特性，ALP熱膨脹割線模量
MPDELE,ALPY,3
MPDELE,ALPZ,3
MPDATA,ALPX,3,,13.16e-6
MPTEMP,1,20,100,200,300,400,500,700,750,1200,3000
MPDATA,EX,1,1,2.02e11,1.87e11,1.79e11,1.70e11,1.61e11,1.37e11,0.5e11,0.3e11,5.00e9,5.00e9
MPDATA,NUXY,1,,,0.30,0.30,0.32,0.32,0.32,0.35,0.4,0.5
TB,BISO,1,6,2,
TBTEMP,20
TBDATA,,305e6,0.25e9,,,,
TBTEMP,200
TBDATA,,259e6,0.22e9,,,,
TBTEMP,450
TBDATA,,168e6,0.205e9,,,,
TBTEMP,750
TBDATA,,20e6,0.03e9,,,,
TBTEMP,1200
TBDATA,,3e6,0.10e7,,,,
TBTEMP,3000
TBDATA,,3e6,0.10e7,,,,
MPTEMP,1,20,100,200,300,400,500,600,700,800,900,1000,3000
MPDATA,ALPX,1,1,11.00e-6,11.16e-6,12.12e-6,12.78e-6,13.83e-6,13.93e-6,14.38e-6,14.81e-6,12.93e-6,12.48e-6,13.16e-6,13.16e-6
!!!!!!!!!!!!!!!!!!!!!!!
ETCHG,STT     !!!!!!結構到熱
finish
/POST1
SET,1,LAST,1,     !!!!!!從結果文件中讀出所指定的數據集
NSEL,S,TEMP,,1600,3000,,0      !!!!!!選擇一組節點子集
ESLN,S                         !!!!!!選擇已選節點上的單元
finish
/PREP7
EMODIF,ALL,MAT,1,              !!!!!!調整單元坐標系方向
ALLSEL,ALL
ETCHG,TTS   !!!熱到結構
finish
/SOLU
TREF,20,      !!!!!!定義參考溫度
DA,1,SYMM     !!!!!!在面上定義約束條件,面1，SYMM 對稱
FLST,2,1,5,ORDE,1
FITEM,2,2
!*
/GO
DA,P51X,ALL,
LDREAD,TEMP,,,0.25,,'conv','rth',''    !!!!!!施加耦合場載荷
ANTYPE,4
NLGEOM,1      !!!!!!在靜態分析或完全瞬態分析中(1)包含大變形效應,幾何非線性 
NSUBST,5,0,0  !!!!!!直接指定子步數
KBC,0         !!!!!!荷載連續增加
LNSRCH,1      !!!!!!激活一個線性搜索與Newton-Raphson一起使用
TIME,0.25
NLGEOM,1
SSTIF,1       !!!!!!在非線性分析中激活應力剛度效應
NROPT,AUTO,,OFF      !!!!!!指定Newton-Raphson法求解的選項,AUTO:程序選擇,OFF：不使用自適應下降因子
EQSLV,,,0,           !!!!!!指定一個方程求解器,0?  MULT：在收斂極端中，用來控制所完成最大迭代次數的乘數，取值范圍為1到3，1是表                     !!!!!!示關閉求解控制。一般取2 
PRECISION,0          !!!!!!對求解器指定機器的精度，0為雙精度
MSAVE,0              !!!!!!設置求解器內存節省選項，0對整個模型的剛度矩陣使用總體組裝
TOFFST,273,          !!!!!!0度
OUTRES,ALL,LAST,     !!!!!!控制數據庫和結果文件中的記錄內容
solve
finish
!
*do,i,2,5
/PREP7
ETCHG,STT
finish
/POST1
SET,i,LAST,1,
NSEL,S,TEMP,,2000,3000,,0
ESLN,S
finish
!PREP7
EMODIF,ALL,MAT,1,
ALLSEL,ALL
ETCHG,TTS
finish
/SOLU
ANTYPE,,REST
LDREAD,TEMP,,,0.25*i,,'conv','rth',''
TIME,0.25*i
OUTRES,ALL,LAST,
solve
finish
*enddo





*do,i,6,101
/PREP7
ETCHG,STT
finish
/POST1
SET,i,LAST,1,
NSEL,S,TEMP,,2300,3000,,0
ESLN,S
finish
EMODIF,ALL,MAT,1,
/PREP7
ALLSEL,ALL
ETCHG,TTS
finish
/SOLU
ANTYPE,,REST
LDREAD,TEMP,,,0.25*i,,'conv','rth',''
TIME,0.25*i
OUTRES,ALL,LAST,
solve
finish
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!
*do,i,102,181
/SOLU
ANTYPE,,REST
LDREAD,TEMP,,,0.25*i,,'conv','rth',''
TIME,0.25*i
OUTRES,ALL,LAST,
solve
finish
*enddo
!!!!!!!!!!!!!!!應力場分析