焊接.txt

关于焊接模拟的ANSYS程序

FINISH
/CLEAR,START
/BATCH  
/CWD,'d:\Ansyswork\EBM FEM Simulation for Metals\3D Temperture distrbution' 
/FILNAME,MatCha_db,0
/TITLE,Change material performance db Method
!选择分析问题的类型：应力
!结构计算，在计算前应删除表面载荷并转换单元类型
/UNITS,SI
/PREP7                                                     !进入前处理
!0参数定义
!模型参数
!几何模型尺寸me-3=mm
*SET,Block_x,7e-3    
*SET,Block_y,3e-3
*SET,Block_z,0.2e-3
!划分网格的大小
*SET,LMeshx,1e-4  !line1 
*SET,LMeshy,1e-4  !line4
*SET,LMeshz,1e-4  !line9
!边界条件参数
*SET,Tinit,30        !初始温度
!1定义单元类型
ET,1,SOLID45       !or solid185
!2实常数
!Nothing
!3材料性能
toffst,273    !从绝对温度改为摄氏度

*DO,MNUM,1,2
MPTEMP  
MPTEMP,1,24,400,1200,1600,2000,2500 
MPDATA,DENS,MNUM,1,3983,4688.4,5406.8,5948.8,6908,6331          !密度
MPDATA,KXX,MNUM,1,1.6278,2.3812,3.6396,29.8,30.6,32        !热导率 
mpdata,c,MNUM,1,472,535,679,859,945,1034                !比热

MPTEMP,,,,,,,,  
MPTEMP,1,0,200,400,800,1200 
MPTEMP,6,1400,1454,1800,2000 
MPDATA,ENTH,MNUM,1,0,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9 
MPDATA,ENTH,MNUM,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10        !热焓

MPTEMP,1,20,300,600,750,900,1200,1500,1800,2000
MPDATA,HF,MNUM,1,6,50,120,180,200,250,378,700,850 !对流系数

mp,murx,MNUM,1                    !***输入磁场相对渗透系数
mp,perx,MNUM,1                    !介质常数(PERY,PERZ)

MPTEMP,,,,,,,,  
MPTEMP,1,20,200,400,600,800  
MPTEMP,6,900,1000,1100,1200,1420,1460,3500
MPDATA,ALPX,MNUM,1,15.24e-6,16.43e-6,17.44e-6,18.21e-6,18.83e-6         !热膨胀系数
MPDATA,alpx,MNUM,6,19.11e-6,19.38e-6,19.66e-6,19.95e-6,20.7,20.7,20.7 

MPTEMP  
MPTEMP,1,20,100,200,300,400 
MPTEMP,6,500,600,700,800,1000   
MPTEMP,11,1200,1400,3500
MPDATA,EX,MNUM,1,2.04E+11,1.96E+11,1.87E+11,1.81E+11,1.71E+11  !   
MPDATA,EX,MNUM,6,1.6E+11,1.52E+11,1.44E+11,1.29E+11,6E+10  
MPDATA,EX,MNUM,11,2E+10,1E+10,1E+10         !杨氏（弹性）模量 


MPTEMP 
MPTEMP,1,20,200,400,600,800
MPTEMP,6,900,1000,1100,1200,1420,1460,3500
MPDATA,PRXY,MNUM,0.267,0.29,0.322,0.296,0.262,0.24,0.229,0.223,0.223,0.223,0.223,0.223 !泊松比

!**应力应变求解设置**
TB,MISO,MNUM,13,4                ! Activate a data table
TBTEMP,20                   ! Temperature = 20
TBPT,DEFI,9.5931e-4,1.957e8    ! Strain, stress at temperature =2 0
TBPT,DEFI,2E-3,        2.86e8                !原来是2.06
TBPT,DEFI,1E-2,        3.11e8
TBPT,DEFI,0.2,        5.2e8

TBTEMP,100                   ! Temperature =100
TBPT,DEFI,7.4643e-4,1.463e8    ! Strain, stress at temperature =100
TBPT,DEFI,2E-3,        2.04e8                !原来是1.54
TBPT,DEFI,1E-2,        2.52e8
TBPT,DEFI,0.2,        4.47e8

TBTEMP,200                   ! Temperature =200
TBPT,DEFI,6.5535e-4,1.2255e8    ! Strain, stress at temperature =200
TBPT,DEFI,2E-3,        1.54e8                
TBPT,DEFI,1E-2,        2.19e8
TBPT,DEFI,0.2,        4.00e8


TBTEMP,300                   ! Temperature =300
TBPT,DEFI,6.0359e-4,1.0925e8    ! Strain, stress at temperature =300
TBPT,DEFI,2E-3,        1.15e8
TBPT,DEFI,1E-2,        2.09e8
TBPT,DEFI,0.2,        3.97e8

TBTEMP,400                   ! Temperature =400
TBPT,DEFI,5.8889e-4,1.007e8   ! Strain, stress at temperature =400
TBPT,DEFI,2E-3,        1.06e8
TBPT,DEFI,1E-2,        2.03e8
TBPT,DEFI,0.2,        3.97e8

TBTEMP,500                   ! Temperature =500
TBPT,DEFI,5.9375e-4,9.5e7    ! Strain, stress at temperature =500
TBPT,DEFI,2E-3,        1.00e8
TBPT,DEFI,1E-2,        1.90e8
TBPT,DEFI,0.2,        3.71e8

TBTEMP,600                   ! Temperature =600
TBPT,DEFI,5.8125e-4,8.835e7    ! Strain, stress at temperature =600
TBPT,DEFI,2E-3,        9.30e7
TBPT,DEFI,1E-2,        1.66e8
TBPT,DEFI,0.2,        3.11e8

TBTEMP,700                   ! Temperature =700
TBPT,DEFI,5.3438e-4,7.695e7    ! Strain, stress at temperature =700
TBPT,DEFI,2E-3,        8.10e7
TBPT,DEFI,1E-2,        1.27e8
TBPT,DEFI,0.2,        2.18e8

TBTEMP,800                   ! Temperature =800
TBPT,DEFI,4.4186e-4,5.7e7    ! Strain, stress at temperature =800
TBPT,DEFI,2E-3,        6e7
TBPT,DEFI,1E-2,        7.63e7
TBPT,DEFI,0.2,        1.09e8


TBTEMP,1000                   ! Temperature =1000
TBPT,DEFI,4.75e-4,2.85e7    ! Strain, stress at temperature =1000
TBPT,DEFI,2E-3,        3e7
TBPT,DEFI,1E-2,        4e7
TBPT,DEFI,0.2,        6e7

TBTEMP,1200                   ! Temperature =1200
TBPT,DEFI,4.75e-4,9.5e6    ! Strain, stress at temperature =1200
TBPT,DEFI,2E-3,        10e6
TBPT,DEFI,1E-2,        13.3e6
TBPT,DEFI,0.2,        2e7

TBTEMP,1400                   ! Temperature =1400
TBPT,DEFI,4.75e-4,4.75e6    ! Strain, stress at temperature =1400
TBPT,DEFI,2E-3,        5e6
TBPT,DEFI,1E-2,        6.67e6
TBPT,DEFI,0.2,        1e7

TBTEMP,3500                   ! Temperature =2000
TBPT,DEFI,4.75e-4,4.75e6    ! Strain, stress at temperature =3000
TBPT,DEFI,2E-3,        4.55e6
TBPT,DEFI,1E-2,        6.e6
TBPT,DEFI,0.2,        0.9e7
*ENDDO        
mp,reft,1,24
mp,reft,2,1400
!4几何模型
BLOCK,0,Block_x,0,Block_y,0,Block_z,
!5网格划分控制
!ESIZE,1e-4,
lesize,1,LMeshx,
lesize,4,LMeshy,
lesize,9,LMeshz,

MSHKEY,1
VMESH,ALL


*set,F,10e-3            !10mm/s
*SET,dt,LMeshx/F          !时间步长 即每步加热的时间        

/prep7
etchg,stt

ldread,temp,,,,,,rth

/post1
*dim,TB,char,120
TB(1)='TB1','TB2','TB3','TB4','TB5','TB6','TB7','TB8','TB9','TB10'
!TB(11)='TB11','TB12','TB13','TB14','TB15','TB16','TB17','TB18','TB19','TB20'
!TB(21)='TB21','TB22','TB23','TB24','TB25','TB26','TB27','TB28','TB29','TB30'
!TB(31)='TB31','TB32','TB33','TB34','TB35','TB36','TB37','TB38','TB39','TB40'
!TB(41)='TB41','TB42','TB43','TB44','TB45','TB46','TB47','TB48','TB49','TB50'
!TB(51)='TB51','TB52','TB53','TB54','TB55','TB56','TB57','TB58','TB59','TB60'
!TB(61)='TB61','TB62','TB63','TB64','TB65','TB66','TB67','TB68','TB69','TB70'
!TB(71)='TB71','TB72','TB73','TB74','TB75','TB76','TB77','TB78','TB79','TB80'
!TB(81)='TB81','TB82','TB83','TB84','TB85','TB86','TB87','TB88','TB89','TB90'
!TB(91)='TB91','TB92','TB93','TB94','TB95','TB96','TB97','TB98','TB99','TB100'
!TB(101)='TB101','TB102','TB103','TB104','TB105','TB106','TB107','TB108','TB109','TB110'
!TB(111)='TB111','TB112','TB113','TB114','TB115','TB116','TB117','TB118','TB119','TB120'
*do,J,1,120
set,J,last,1
etable,TB(J),temp
*enddo
finish

/VIEW,1,1,2,3   
/ANG,1  
/REP,FAST 

/prep7
etchg,tts

/solu
ANTYPE,TRANS
nlgeom,on
tref,24
cnvtol,u,,0.5
CUTCONTROL,PLSLIMIT, 10,
lnsrch,on
nropt,full
neqit,2500
EQSLV,ICCG,,2
kbc,0        
TIMINT,ON,STRUC
autots,on
outres,all
outpr,all,last
NSUBST,10,25,5
!NSUBST,1,1,1

local,11,0,0,0,,,,,,
nsel,s,loc,z,0
nsel,r,loc,y,Block_y/2
nsel,r,loc,x,Block_x/2
d,all,all,0
allsel                                        !约束自由度

/solu
*set,tim,0
*set,i,0
!进入主循环
*DO,Y,1,3           !三道扫描线
  *do,X,1,40       !扫描线长度5mm

tim=tim+dt
i=i+1
time,tim

allsel
ESEL,S,ETAB,TB(i),1395
mpchg,2,all
ekill,all                        !杀高温单元
esel,inve
ealive,all
eplot
allsel

ldread,temp,i,,,,,rth  !读入时间tim时的温度分布 ,（i步时）
/pbf,temp,,1                !Shows magnitude of body force loads on displays.
solve

eplot
/auto,1
/replot

  *enddo
*ENDDO 

save