24小时热门版块排行榜

返回列表

liulingsu

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 106.2
帖子: 26
在线: 12.9小时
虫号: 34032179
注册: 2023-08-07
专业: 海洋工程

[交流] apdl热分析遇到一些问题已有1人参与

大家好!我在弄平板对接多道焊的热结构耦合分析，但是遇到一些问题，第一道焊缝和第二道焊缝冷却后最低温度低于室温，有的还变成了负的，第三道焊缝冷却后温度是正常的，请问这是什么原因呢？下面是我的命令流，可以指导一下吗？

/FILNAME,fuxian,0
/PREP7
/UNITS,SI

!定义焊接参数
ET,1,SOLID70 !!热分析单元

TOFFST,273
MPTEMP,,,,,,,,    !定义温度
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,210
MPTEMP,4,330
MPTEMP,5,420
MPTEMP,6,540
MPTEMP,7,660
MPTEMP,8,780
MPTEMP,9,985
MPTEMP,10,1320
MPTEMP,11,1420
MPTEMP,12,1500

MPDATA,KXX,1,,51  !定义导热系数
MPDATA,KXX,1,,50
MPDATA,KXX,1,,49
MPDATA,KXX,1,,46
MPDATA,KXX,1,,41
MPDATA,KXX,1,,38
MPDATA,KXX,1,,34
MPDATA,KXX,1,,28
MPDATA,KXX,1,,25
MPDATA,KXX,1,,32
MPDATA,KXX,1,,42
MPDATA,KXX,1,,42

MPDATA,C,1,,450 !定义比热容
MPDATA,C,1,,475
MPDATA,C,1,,530
MPDATA,C,1,,560
MPDATA,C,1,,630
MPDATA,C,1,,720
MPDATA,C,1,,830
MPDATA,C,1,,910
MPDATA,C,1,,1055
MPDATA,C,1,,2000
MPDATA,C,1,,2100
MPDATA,C,1,,2150

MPDATA,DENS,1,,7860 !定义密度
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860
MPDATA,DENS,1,,7860

MPDATA,EX,1,,210.0E9 !定义弹性模量
MPDATA,EX,1,,195.0E9
MPDATA,EX,1,,195.0E9
MPDATA,EX,1,,185.0E9
MPDATA,EX,1,,168.0E9
MPDATA,EX,1,,118.0E9
MPDATA,EX,1,,52.0E9
MPDATA,EX,1,,12.0E9
MPDATA,EX,1,,11.8E9
MPDATA,EX,1,,10.4E9
MPDATA,EX,1,,10.2E9
MPDATA,EX,1,,10E9

MPDATA,PRXY,1,,0.260  !定义泊松比
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260
MPDATA,PRXY,1,,0.260

UIMP,1,REFT,,,20
MPDATA,ALPX,1,,11.2E-6  !定义膨胀系数
MPDATA,ALPX,1,,11.8E-6
MPDATA,ALPX,1,,12.4E-6
MPDATA,ALPX,1,,13.1E-6
MPDATA,ALPX,1,,13.6E-6
MPDATA,ALPX,1,,14.1E-6
MPDATA,ALPX,1,,14.6E-6
MPDATA,ALPX,1,,14.6E-6
MPDATA,ALPX,1,,14.6E-6
MPDATA,ALPX,1,,14.6E-6
MPDATA,ALPX,1,,14.6E-6
MPDATA,ALPX,1,,14.6E-6

TB,BISO,1,6
TBTEMP,20
TBDATA,,380E+6,210.0E8
TBTEMP,210
TBDATA,,320E+6,195.0E8
TBTEMP,420
TBDATA,,190E+6,168.0E8
TBTEMP,660
TBDATA,,75E+6,52.0E8
TBTEMP,1320
TBDATA,,28E+6,10.4E8
TBTEMP,1500
TBDATA,,20E+6,10E8

!!建立几何模型
L=300E-3 !焊件的长度
W=250E-3  !焊件的宽度
H=6E-3    !焊件的高度
H1=4E-3    !第一道焊缝高度
H2=2E-3    !第二道和第三道焊缝高度
Ww=24E-3 !焊缝宽度

BLOCK,0,Ww/2,0,-H,0,W,
BLOCK,0,-Ww/2,0,-H,0,W,

BLOCK,Ww/2,Ww/2+30e-3,0,-H,0,W,
BLOCK,Ww/2+30e-3,Ww/2+40e-3, 0,-H,0,W,
BLOCK,Ww/2+40e-3,L/2,0,-H,0,W,

BLOCK,-Ww/2,-Ww/2-30e-3,0,-H,0,W,
BLOCK,-Ww/2-30e-3,-Ww/2-40e-3, 0,-H,0,W,
BLOCK,-Ww/2-40e-3,-L/2,0,-H,0,W,

wpro,,90.000000,
wpof,,,H2
RECTNG,-0.3,0.3,-0.3,0.3,
FLST,2,8,6,ORDE,2
FITEM,2,1
FITEM,2,-8
VSBA,P51X,    49
VGLUE,ALL
WPCSYS,,0

!!划分网格
LESIZE,12,,,W/2e-3
LESIZE,112,,,H2/2e-3
LESIZE,104,,,H1/2e-3
FLST,2,2,4,ORDE,2
FITEM,2,8
FITEM,2,103
LESIZE,P51X,,,Ww/2/2e-3

VMESH,4
VMESH,10
VMESH,1
VMESH,9

ESIZE,2E-3
FLST,2,4,6,ORDE,4
FITEM,2,2
FITEM,2,3
FITEM,2,5
FITEM,2,6
VMESH,P51X

ESIZE,3E-3
FLST,2,4,6,ORDE,4
FITEM,2,7
FITEM,2,8
FITEM,2,27
FITEM,2,28
VMESH,P51X

ESIZE,5E-3
FLST,2,4,6,ORDE,4
FITEM,2,25
FITEM,2,26
FITEM,2,29
FITEM,2,30
VSWEEP,P51X

!!温度场
!第一道焊缝参数
v1=0.22/60 !焊接速度
I1=190 !电流
U1=22    !电压
!第二道焊缝参数
v2=0.55/60 !焊接速度
I2=220 !电流
U2=27    !电压
!第三道焊缝参数
v3=0.52/60 !焊接速度
I3=220 !电流
U3=27    !电压

!双椭球热源参数
ff=0.6
fr=1.4
a1=4e-3
a2=3e-3
a3=3e-3
b1=4e-3
b2=3e-3
b3=3e-3
cf=4e-3
cr=8e-3
Q1=0.75*U1*I1/(100*v1)
Q2=0.75*U2*I2/(100*v2)
Q3=0.75*U3*I3/(100*v3)

local,11,0,0,-H2,0,,,,1,1
local,12,0,-Ww/4,0,0,,,,1,1
local,13,0,Ww/4,0,0,,,,1,1

/SOL
ANTYPE,4
TRNOPT,FULL                         !瞬态分析选项，完全分析
TUNIF,20,                                  !初始温度
NROPT,FULL,,ON                      !定义完全牛顿-拉普森方法，激活自适应下降（默认）
KBC,1
AUTOTS,ON                               !自动步长
LNSRCH,ON                               !线性搜索
OUTRES,ALL ,ALL
PRED,ON                                  !打开预测校正
TIMINT,ON                               !时间积分设置
TINTP,,,,1,,                                  !定义瞬态综合参数
TREF,20                                     !参考温度

!把焊缝区的单元杀死
FLST,5,4,6,ORDE,4
FITEM,5,1
FITEM,5,4
FITEM,5,9
FITEM,5,-10
VSEL,S, , ,P51X
NSLV,R,1
ESLN,R,1,
EKILL,ALL
ALLSEL,ALL

!第一道焊缝
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*SET,_FNCNAME,'func1'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),ff
*SET,_FNC_C2(1),Q1
*SET,_FNC_C3(1),a1
*SET,_FNC_C4(1),b1
*SET,_FNC_C5(1),cf
*SET,_FNCCSYS,11
! /INPUT,func1.func,,,1
*DIM,%_FNCNAME%,TABLE,7,42,1,,,,%_FNCCSYS%
!
! Begin of equation: 6*sqrt(3)*ff*Q/({PI}*sqrt({PI})*a*b*cf)*exp(-3*{x}^2/a^2-3
! *{y}^2/b^2-3*{Z}^2/cf^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), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 16, 3.14159265358979312, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 3.14159265358979312, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, -4, 3, 19
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -1, 3, 20
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, -3, 3, 21
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -2, 4, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -2, 3, -4
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,22,1), 0.0, -4, 0, 1, 19, 17, -2
*SET,%_FNCNAME%(0,23,1), 0.0, -2, 0, 1, -1, 4, -4
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,25,1), 0.0, -4, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -4
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,29,1), 0.0, -4, 0, 1, 20, 17, -1
*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 1, -5, 4, -4
*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, -2, 2, -1
*SET,%_FNCNAME%(0,32,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 1, 4, 17, -1
*SET,%_FNCNAME%(0,34,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,37,1), 0.0, -2, 0, 1, 21, 17, -1
*SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, -5, 4, -2
*SET,%_FNCNAME%(0,39,1), 0.0, -2, 0, 1, -4, 2, -1
*SET,%_FNCNAME%(0,40,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,41,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,42,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 6*sqrt(3)*ff*Q/({PI}*sqrt({PI})*a*b*cf)*exp(-3*{x}^2/a^2-3*
! {y}^2/b^2-3*{Z}^2/cf^2)
!-->
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*SET,_FNCNAME,'func2'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),fr
*SET,_FNC_C2(1),Q1
*SET,_FNC_C3(1),a1
*SET,_FNC_C4(1),b1
*SET,_FNC_C5(1),cr
*SET,_FNCCSYS,11
! /INPUT,func2.func,,,1
*DIM,%_FNCNAME%,TABLE,7,42,1,,,,%_FNCCSYS%
!
! Begin of equation: 6*sqrt(3)*fr*Q/({PI}*sqrt({PI})*a*b*cr)*exp(-3*{x}^2/a^2-3
! *{y}^2/b^2-3*{Z}^2/cr^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), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 16, 3.14159265358979312, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 3.14159265358979312, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, -4, 3, 19
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -1, 3, 20
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, -3, 3, 21
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -2, 4, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -2, 3, -4
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,22,1), 0.0, -4, 0, 1, 19, 17, -2
*SET,%_FNCNAME%(0,23,1), 0.0, -2, 0, 1, -1, 4, -4
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,25,1), 0.0, -4, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -4
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,29,1), 0.0, -4, 0, 1, 20, 17, -1
*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 1, -5, 4, -4
*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, -2, 2, -1
*SET,%_FNCNAME%(0,32,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 1, 4, 17, -1
*SET,%_FNCNAME%(0,34,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,37,1), 0.0, -2, 0, 1, 21, 17, -1
*SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, -5, 4, -2
*SET,%_FNCNAME%(0,39,1), 0.0, -2, 0, 1, -4, 2, -1
*SET,%_FNCNAME%(0,40,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,41,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,42,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 6*sqrt(3)*fr*Q/({PI}*sqrt({PI})*a*b*cr)*exp(-3*{x}^2/a^2-3*
! {y}^2/b^2-3*{Z}^2/cr^2)
!-->

!第二道焊缝
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*SET,_FNCNAME,'func3'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),ff
*SET,_FNC_C2(1),Q2
*SET,_FNC_C3(1),a2
*SET,_FNC_C4(1),b2
*SET,_FNC_C5(1),cf
*SET,_FNCCSYS,12
! /INPUT,func3.func,,,1
*DIM,%_FNCNAME%,TABLE,7,42,1,,,,%_FNCCSYS%
!
! Begin of equation: 6*sqrt(3)*ff*Q/({PI}*sqrt({PI})*a*b*cf)*exp(-3*{x}^2/a^2-3
! *{y}^2/b^2-3*{Z}^2/cf^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), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 16, 3.14159265358979312, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 3.14159265358979312, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, -4, 3, 19
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -1, 3, 20
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, -3, 3, 21
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -2, 4, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -2, 3, -4
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,22,1), 0.0, -4, 0, 1, 19, 17, -2
*SET,%_FNCNAME%(0,23,1), 0.0, -2, 0, 1, -1, 4, -4
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,25,1), 0.0, -4, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -4
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,29,1), 0.0, -4, 0, 1, 20, 17, -1
*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 1, -5, 4, -4
*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, -2, 2, -1
*SET,%_FNCNAME%(0,32,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 1, 4, 17, -1
*SET,%_FNCNAME%(0,34,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,37,1), 0.0, -2, 0, 1, 21, 17, -1
*SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, -5, 4, -2
*SET,%_FNCNAME%(0,39,1), 0.0, -2, 0, 1, -4, 2, -1
*SET,%_FNCNAME%(0,40,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,41,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,42,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 6*sqrt(3)*ff*Q/({PI}*sqrt({PI})*a*b*cf)*exp(-3*{x}^2/a^2-3*
! {y}^2/b^2-3*{Z}^2/cf^2)
!-->
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*SET,_FNCNAME,'func4'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),fr
*SET,_FNC_C2(1),Q2
*SET,_FNC_C3(1),a2
*SET,_FNC_C4(1),b2
*SET,_FNC_C5(1),cr
*SET,_FNCCSYS,12
! /INPUT,func4.func,,,1
*DIM,%_FNCNAME%,TABLE,7,42,1,,,,%_FNCCSYS%
!
! Begin of equation: 6*sqrt(3)*fr*Q/({PI}*sqrt({PI})*a*b*cr)*exp(-3*{x}^2/a^2-3
! *{y}^2/b^2-3*{Z}^2/cr^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), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 16, 3.14159265358979312, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 3.14159265358979312, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, -4, 3, 19
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -1, 3, 20
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, -3, 3, 21
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -2, 4, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -2, 3, -4
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,22,1), 0.0, -4, 0, 1, 19, 17, -2
*SET,%_FNCNAME%(0,23,1), 0.0, -2, 0, 1, -1, 4, -4
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,25,1), 0.0, -4, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -4
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,29,1), 0.0, -4, 0, 1, 20, 17, -1
*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 1, -5, 4, -4
*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, -2, 2, -1
*SET,%_FNCNAME%(0,32,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 1, 4, 17, -1
*SET,%_FNCNAME%(0,34,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,37,1), 0.0, -2, 0, 1, 21, 17, -1
*SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, -5, 4, -2
*SET,%_FNCNAME%(0,39,1), 0.0, -2, 0, 1, -4, 2, -1
*SET,%_FNCNAME%(0,40,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,41,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,42,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 6*sqrt(3)*fr*Q/({PI}*sqrt({PI})*a*b*cr)*exp(-3*{x}^2/a^2-3*
! {y}^2/b^2-3*{Z}^2/cr^2)
!-->

!第三道焊缝
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*SET,_FNCNAME,'func5'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),ff
*SET,_FNC_C2(1),Q3
*SET,_FNC_C3(1),a3
*SET,_FNC_C4(1),b3
*SET,_FNC_C5(1),cf
*SET,_FNCCSYS,13
! /INPUT,func5.func,,,1
*DIM,%_FNCNAME%,TABLE,7,42,1,,,,%_FNCCSYS%
!
! Begin of equation: 6*sqrt(3)*ff*Q/({PI}*sqrt({PI})*a*b*cf)*exp(-3*{x}^2/a^2-3
! *{y}^2/b^2-3*{Z}^2/cf^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), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 16, 3.14159265358979312, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 3.14159265358979312, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, -4, 3, 19
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -1, 3, 20
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, -3, 3, 21
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -2, 4, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -2, 3, -4
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,22,1), 0.0, -4, 0, 1, 19, 17, -2
*SET,%_FNCNAME%(0,23,1), 0.0, -2, 0, 1, -1, 4, -4
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,25,1), 0.0, -4, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -4
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,29,1), 0.0, -4, 0, 1, 20, 17, -1
*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 1, -5, 4, -4
*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, -2, 2, -1
*SET,%_FNCNAME%(0,32,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 1, 4, 17, -1
*SET,%_FNCNAME%(0,34,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,37,1), 0.0, -2, 0, 1, 21, 17, -1
*SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, -5, 4, -2
*SET,%_FNCNAME%(0,39,1), 0.0, -2, 0, 1, -4, 2, -1
*SET,%_FNCNAME%(0,40,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,41,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,42,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 6*sqrt(3)*ff*Q/({PI}*sqrt({PI})*a*b*cf)*exp(-3*{x}^2/a^2-3*
! {y}^2/b^2-3*{Z}^2/cf^2)
!-->
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*SET,_FNCNAME,'func6'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),fr
*SET,_FNC_C2(1),Q3
*SET,_FNC_C3(1),a3
*SET,_FNC_C4(1),b3
*SET,_FNC_C5(1),cr
*SET,_FNCCSYS,13
! /INPUT,func6.func,,,1
*DIM,%_FNCNAME%,TABLE,7,42,1,,,,%_FNCCSYS%
!
! Begin of equation: 6*sqrt(3)*fr*Q/({PI}*sqrt({PI})*a*b*cr)*exp(-3*{x}^2/a^2-3
! *{y}^2/b^2-3*{Z}^2/cr^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), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 16, 3.14159265358979312, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 3.14159265358979312, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, -4, 3, 19
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -1, 3, 20
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, -3, 3, 21
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -2, 4, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -2, 3, -4
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,22,1), 0.0, -4, 0, 1, 19, 17, -2
*SET,%_FNCNAME%(0,23,1), 0.0, -2, 0, 1, -1, 4, -4
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,25,1), 0.0, -4, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -4
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,29,1), 0.0, -4, 0, 1, 20, 17, -1
*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 1, -5, 4, -4
*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, -2, 2, -1
*SET,%_FNCNAME%(0,32,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 1, 4, 17, -1
*SET,%_FNCNAME%(0,34,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,37,1), 0.0, -2, 0, 1, 21, 17, -1
*SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, -5, 4, -2
*SET,%_FNCNAME%(0,39,1), 0.0, -2, 0, 1, -4, 2, -1
*SET,%_FNCNAME%(0,40,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,41,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,42,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 6*sqrt(3)*fr*Q/({PI}*sqrt({PI})*a*b*cr)*exp(-3*{x}^2/a^2-3*
! {y}^2/b^2-3*{Z}^2/cr^2)
!-->

!!!第一道焊缝
t1=W/v1    !第一道焊缝时间
n1=W/2E-3 !载荷步数
n_L1=W/n1       !一个载荷步长度
T=n_L1/v1             !一个载荷步所需时间

!!随着热源移动，单元依次被激活
!一道焊时随着热源移动，单元依次被激活
*do,i,1,n1,1
TIME,T*i
NSUBST,1
local,11,0,0,-H2,n_L1*i,,,,1,1
vsel,s,,,1
vsel,a,,,9
nslv,r,1
nsel,r,loc,z,0,-n_L1
esln,r,1,all
ealive,all             !复活部分单元
allsel,all

vsel,s,,,1
vsel,a,,,9
nslv,r,1
nsel,r,loc,z,0,W
bf,all,hgen,%func1%

vsel,s,,,1
vsel,a,,,9
nslv,r,1
nsel,r,loc,z,-W,0
bf,all,hgen,%func2%
allsel

esel,s,live
nsle,s
nsel,r,ext
SF,all,CONV,15,25,       !!对流系数
SF,all,RDSF,0.9,1,       !!辐射系数
spctemp,1,25
stef,5.67e-8
toffst,273
ALLSEL,ALL
SOLVE
!删除移动热源载荷
BFDELE,ALL,HGEN
SFDELE,ALL,CONV
SFDELE,ALL,RDSF
*ENDDO
ALLSEL

!冷却
esel,s,live
nsle,s
nsel,r,ext
SF,all,CONV,15,25    !!对流系数
SF,all,RDSF,0.9,1,       !!辐射系数
spctemp,1,25
stef,5.67e-8
toffst,273
allsel,all
TIME,400
NSUBST,400,200,500
ALLSEL
SOLVE
SFDELE,ALL,CONV
SFDELE,ALL,RDSF
ALLSEL

!!!第二道焊缝
t2=W/v2    !第一道焊缝时间
n2=W/2E-3 !载荷步数
n_L2=W/n2       !一个载荷步长度
TT=n_L2/v2             !一个载荷步所需时间

!!随着热源移动，单元依次被激活
!一道焊时随着热源移动，单元依次被激活
*do,i,1,n2,1
TIME,400+TT*i
NSUBST,1
local,12,0,-WW/4,0,n_L2*i,,,,1,1
vsel,s,,,4
nslv,r,1
nsel,r,loc,z,0,-n_L2
esln,r,1,all
ealive,all             !复活部分单元
allsel,all

vsel,s,,,4
nslv,r,1
nsel,r,loc,z,0,W
bf,all,hgen,%func3%

vsel,s,,,4
nslv,r,1
nsel,r,loc,z,-W,0
bf,all,hgen,%func4%
allsel

esel,s,live
nsle,s
nsel,r,ext
SF,all,CONV,15,25,       !!对流系数
SF,all,RDSF,0.9,1,       !!辐射系数
spctemp,1,25
stef,5.67e-8
toffst,273
ALLSEL,ALL
SOLVE
!删除移动热源载荷
BFDELE,ALL,HGEN
SFDELE,ALL,CONV
SFDELE,ALL,RDSF
*ENDDO
ALLSEL

!冷却
esel,s,live
nsle,s
nsel,r,ext
SF,all,CONV,15,25    !!对流系数
SF,all,RDSF,0.9,1,       !!辐射系数
spctemp,1,25
stef,5.67e-8
toffst,273
allsel,all
TIME,800
NSUBST,400,200,500
ALLSEL
SOLVE
SFDELE,ALL,RDSF
SFDELE,ALL,CONV
ALLSEL

!!!第三道焊缝
t3=W/v3    !第一道焊缝时间
n3=W/2E-3 !载荷步数
n_L3=W/n3       !一个载荷步长度
TTT=n_L3/v3             !一个载荷步所需时间

!!随着热源移动，单元依次被激活
!一道焊时随着热源移动，单元依次被激活
*do,i,1,n3,1
TIME,800+TTT*i
NSUBST,1
local,13,0,WW/4,0,n_L3*i,,,,1,1
vsel,s,,,10
nslv,r,1
nsel,r,loc,z,0,-n_L2
esln,r,1,all
ealive,all             !复活部分单元
allsel,all

vsel,s,,,10
nslv,r,1
nsel,r,loc,z,0,W
bf,all,hgen,%func5%

vsel,s,,,10
nslv,r,1
nsel,r,loc,z,-W,0
bf,all,hgen,%func6%
allsel

esel,s,live
nsle,s
nsel,r,ext
SF,all,CONV,15,25,       !!对流系数
SF,all,RDSF,0.9,1,       !!辐射系数
spctemp,1,25
stef,5.67e-8
toffst,273
ALLSEL,ALL
SOLVE
!删除移动热源载荷
BFDELE,ALL,HGEN
SFDELE,ALL,CONV
SFDELE,ALL,RDSF
*ENDDO
ALLSEL

!冷却
esel,s,live
nsle,s
nsel,r,ext
SF,all,CONV,15,25    !!对流系数
SF,all,RDSF,0.9,1,       !!辐射系数
spctemp,1,25
stef,5.67e-8
toffst,273
allsel,all
TIME,1500
NSUBST,400,200,500
ALLSEL
SOLVE
SFDELE,ALL,RDSF
SFDELE,ALL,CONV
ALLSEL

发自小木虫Android客户端

回复此楼

» 猜你喜欢

1楼 2023-09-09 00:43:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

应助: 668 (博士)
金币: 9199
红花: 146
帖子: 2205
在线: 346小时
虫号: 728202
注册: 2009-03-22
性别: GG
专业: 动力学与控制

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

程序太长，智能简单看了一下。提几个建议，可以修改一下：
1 命令 KBC, 1 改为 KBC,0
   将阶梯式加载改为斜坡方式加载，更符合实际情况。
2  所有的函数定义，移到使用该函数的命令之前。因为按照 ANSYS 的规定，函数中的所有参数 (包括局部坐标系) 都是使用该函数定义之前的数据，你的函数定义是在求解部分的最前面，其中的参数将固定为该处的值，不会随着热源的移动而变化，与实际情况不符合。因而所有三道焊缝的计算结果都是错误的。
3  多处 nsubst 命令，最好修改一下。该命令的格式是：
      NSUBST,  常规子步数,  最大子步数,  最小子步数
   你的命令流中：
   第一个：  NSUBST,1 建议改为：
      NSUBST,5,5,1
   其它的：  NSUBST,400,200,500  ，不符合该命令的用法，建议改为：
      NSUBST, 10,1000,10
4  你的热辐射的定义可能有问题，由于我手头没有资料，暂时说不出哪里有问题，你可以先不施加热辐射载荷，试试看是否冷却后还出现你说的问题。
以上仅供参考，你不妨试试。

赞一下(1人)