weld_pipe.txt

圆筒焊接热应力及残余应力分析的ansys有限元命令流

!!!!圆筒焊接热应力及残余应力分析命令流
/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
!!!!!!!!!!!!!!!应力场分析