/FILNAME,Welding residual stress,0
!进入前处理
/PREP7
ET,1,SOLID226,11 !定义热-结构耦合单元
ET,2,SURF152 !定义结构的表面效应单元
KEYOPT,2,8,1 !设置表面效应单元关键字
U=22 !焊接电压
I=150 !焊接电流
V=0.01 !焊接速度
YITA=0.7 !焊接热效率
R=0.008 !电弧有效加热半径
Q=U*I*YITA !电弧热功率
Qm=3/3.1415/R**2*Q !加热斑点中心最大热流密度
!设置材料参数
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1700
MPTEMP,8,2500
MPDATA,KXX,1,,50
MPDATA,KXX,1,,47
MPDATA,KXX,1,,40
MPDATA,KXX,1,,27
MPDATA,KXX,1,,30
MPDATA,KXX,1,,35
MPDATA,KXX,1,,40
MPDATA,KXX,1,,55
MPDATA,DENS,1,,7820
MPDATA,DENS,1,,7700
MPDATA,DENS,1,,7610
MPDATA,DENS,1,,7550
MPDATA,DENS,1,,7490
MPDATA,DENS,1,,7350
MPDATA,DENS,1,,7300
MPDATA,DENS,1,,7090
MPDATA,C,1,,460
MPDATA,C,1,,480
MPDATA,C,1,,530
MPDATA,C,1,,675
MPDATA,C,1,,670
MPDATA,C,1,,660
MPDATA,C,1,,780
MPDATA,C,1,,820
MPDATA,EX,1,,2.05E11
MPDATA,EX,1,,1.87E11
MPDATA,EX,1,,1.5E11
MPDATA,EX,1,,0.7E11
MPDATA,EX,1,,0.2E11
MPDATA,EX,1,,0.19E2
MPDATA,EX,1,,0.18E2
MPDATA,EX,1,,0.12e2
MPDATA,PRXY,1,,0.28
MPDATA,PRXY,1,,0.29
MPDATA,PRXY,1,,0.31
MPDATA,PRXY,1,,0.35
MPDATA,PRXY,1,,0.4
MPDATA,PRXY,1,,0.45
MPDATA,PRXY,1,,0.48
MPDATA,PRXY,1,,0.5
UIMP,1,REFT,,,20
MPDATA,ALPX,1,,1.1e-5
MPDATA,ALPX,1,,1.22e-5
MPDATA,ALPX,1,,1.39e-5
MPDATA,ALPX,1,,1.48e-5
MPDATA,ALPX,1,,1.34e-5
MPDATA,ALPX,1,,1.33e-5
MPDATA,ALPX,1,,1.32e-5
MPDATA,ALPX,1,,1.31e-5
TB,BISO,1,6,2,
TBTEMP,20
TBDATA,,225e6,0,,,,
TBTEMP,250
TBDATA,,185e6,0,,,,
TBTEMP,500
TBDATA,,85e6,0,,,,
TBTEMP,750
TBDATA,,45E6,0,,,,
TBTEMP,1000
TBDATA,,15E6,0,,,,
TBTEMP,1500
TBDATA,,1E-5,0,,,,
!建立模型
BLOCK,0,0.18,0,0.005,0,0.15,
BLOCK,0.18,0.23,0,0.005,0,0.15,
BLOCK,0.23,0.41,0,0.005,0,0.15,
!划分网格
ESIZE,0.005
VMESH,ALL
!生成表面效应单元
NSEL,ALL
NSEL,S,LOC,Y,0.005
NPLOT
TYPE, 2
MAT, 1
REAL,
ESYS, 0
SECNUM,
TSHAP,LINE
ESURF,0
ALLSEL,ALL
!定义焊件之间的接触对
!*
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,
MAT,1
MP,EMIS,1,
R,3
REAL,3
ET,3,170
ET,4,174
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,2E6,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,0
KEYOPT,4,10,2
KEYOPT,4,11,0
KEYOPT,4,12,0
KEYOPT,4,2,0
KEYOPT,3,5,0
KEYOPT,4,1,1
! Generate the target surface
ASEL,S,,,6
CM,_TARGET,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
NSLE,A,CT2
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,11
CM,_CONTACT,AREA
TYPE,4
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,3
ESEL,A,TYPE,,4
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,3
ESEL,A,TYPE,,4
ESEL,R,REAL,,3
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
!*
!*
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0
MAT,1
R,4
REAL,4
ET,5,170
ET,6,174
KEYOPT,6,9,0
KEYOPT,6,10,2
R,4,
RMORE,
RMORE,,2E6
RMORE,0
KEYOPT,6,1,1
! Generate the target surface
ASEL,S,,,12
CM,_TARGET,AREA
TYPE,5
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
NSLE,A,CT2
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,17
CM,_CONTACT,AREA
TYPE,6
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,5
ESEL,A,TYPE,,6
ESEL,R,REAL,,4
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,5
ESEL,A,TYPE,,6
ESEL,R,REAL,,4
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/MREP,EPLOT
)/GOP ! Resume printing after UNDO process
)! We suggest a save at this point
R,3,0,0,1, 0.1,0,0
RMORE,0,0,1e+020, 0,1,0
RMORE,0,2000000,1, 0,1,0.5
RMORE,,,1, 0,0,
RMORE,,,,,,
RMORE,,,,,1000
R,4,0,0,0, 0,0,0
RMORE,0,0,0, 0,0,0
RMORE,0,2000000,0, 0,0,0
RMORE,,,0, 0,0,
RMORE,,,,,,
RMORE,,,,,1000
!进入求解
/SOL
ANTYPE,4 !设置分析类型为瞬态热-结构分析
TRNOPT,FULL !使用完全法完成瞬态分析
NLGEOM,1 !激活大变形分析
!设置求解器
EQSLV,SPAR
TUNIF,20, !设置模型的初始均匀温度
TREF,20, !设置模型的参考温度
!定义对流边界条件
SFA,4, ,CONV,30,20
SFA,10, ,CONV,30,20
SFA,16, ,CONV,30,20
SFA,2, ,CONV,30,20
SFA,8, ,CONV,30,20
SFA,14, ,CONV,30,20
SFA,5, ,CONV,30,20
SFA,18, ,CONV,30,20
SFA,1, ,CONV,30,20
SFA,7, ,CONV,30,20
SFA,13, ,CONV,30,20
SFA,3, ,CONV,300,20
SFA,9, ,CONV,300,20
SFA,15, ,CONV,300,20
!定义结构约束
DA,5,UX,
DA,5,UY,
DA,5,UZ,
DA,18,UX,
DA,18,UY,
DA,18,UZ,
DA,3,UY,
DA,9,UY,
DA,15,UY,
DA,2,UZ,
DA,8,UZ,
DA,14,UZ,
!定义函数
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS
*SET,_FNCNAME,'gaosi1'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),qm
*SET,_FNC_C2(1),v
*SET,_FNC_C3(1),r
*SET,_FNCCSYS,0
! /INPUT,gaosi1.func,,,1
*DIM,%_FNCNAME%,TABLE,6,23,1,,,,%_FNCCSYS%
!
! Begin of equation: Qm*exp(-3*(({X}-0.18)^2+(0.15-V*{TIME}-{Z})^2)/R^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.18, 0, 0, 2
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 2, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -3, 17, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, 18, 3, 1
*SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 0.15, 0, 0, -1
*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, -3, 2, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -5, 2, 4
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, -1, 17, -3
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,17,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 1, 19, 17, -1
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 4, -2
*SET,%_FNCNAME%(0,21,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,23,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: Qm*exp(-3*(({X}-0.18)^2+(0.15-V*{TIME}-{Z})^2)/R^2)
!-->
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS
*SET,_FNCNAME,'GAOSI2'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),QM
*SET,_FNC_C2(1),V
*SET,_FNC_C3(1),R
*SET,_FNCCSYS,0
! /INPUT,GAOSI2.func,,,1
*DIM,%_FNCNAME%,TABLE,6,23,1,,,,%_FNCCSYS%
!
! Begin of equation: Qm*exp(-3*(({X}-0.23)^2+(V*({TIME}-15)-{Z})^2)/R^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.23, 0, 0, 2
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 2, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -3, 17, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 15, 0, 0, 1
*SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, 1, 2, -1
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, 18, 3, -3
*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, -1, 2, 4
*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, -3, 17, -1
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,17,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 1, 19, 17, -1
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 4, -2
*SET,%_FNCNAME%(0,21,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,23,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: Qm*exp(-3*(({X}-0.23)^2+(V*({TIME}-15)-{Z})^2)/R^2)
!-->
!定义第1个载荷步文件
ESEL,ALL
ESEL,S,TYPE,,2
EPLOT
SFE,all,1,HFLUX, , %GAOSI1%
ALLSEL,ALL
OUTRES,ALL,ALL,
TIME,15
AUTOTS,1
NSUBST,75,150,50,1
KBC,0
LSWRITE,1
!定义第2个载荷步文件
ESEL,ALL
ESEL,S,TYPE,,2
EPLOT
SFEDELE,ALL,1,HFLUX
SFE,all,1,HFLUX, , %GAOSI2%
ALLSEL,ALL
TIME,30
AUTOTS,1
NSUBST,75,150,50,1
KBC,0
LSWRITE,2
!定义第3个载荷步文件
ESEL,ALL
ESEL,S,TYPE,,2
EPLOT
SFEDELE,ALL,1,HFLUX
ALLSEL,ALL
TIME,200
AUTOTS,1
NSUBST,600,,50,1
KBC,0
LSWRITE,3
!定义第4个载荷步文件
TIME,3600
AUTOTS,1
NSUBST,100,150,50,1
KBC,0
LSWRITE,4
!使用载荷步文件求解
LSSOLVE,1,4,1 |
回复此楼
为什么 国际(地区)合作与交流项目 没有放榜？
10+3