ANSYS 14.0 实例 搅拌摩擦焊接(焊接)模拟

/FILNAME,Friction Stir Welding,0
/PREP7
/COM, ***** 模型参数 ********

L=76.2E-03 !每一个平板的长度
W=31.75E-03 !每一个平板的宽度
T=3.18E-03 !每一个平板的厚度
R1=7.62E-03 !焊机工具的轴肩半径
H=15.24E-03 !轴肩的高度H
L1=R1 !轴肩中心线在焊接线上的起始位置
L2=L-L1


NR1=3.141593*11 !第二载荷步的转速
NR2=3.141593*45 !第三载荷步的转速
UY1=60.96E-03 !轴肩沿焊接线的平动距离


ET,1,SOLID226,11 !COUPLED-FIELD SOLID ELEMENT,KEYOPT(1) IS
!SET TO 11 FOR A STRUCTURAL-THERMAL ANALYSIS

/COM, * MATERIAL PROPERTIES FOR 304L STAINLESS STEEL PLATES
MP, EX,1,193E9 !弹性模量 (N/M^2)
MP, NUXY,1,0.3 !泊松比
MP, ALPX,1,1.875E-5 !热膨胀系数
MP,QRATE,1,0.8 !摩擦塑性功转换成热的效率80%

/COM, *BISO MATERIAL MODEL
TB,BISO, 1, 1 !双线性等向强化模型
TBDATA,1,290E6,2.8E9 !屈服应力和切线模量

MPTEMP,1, 0, 200, 400, 600, 800,1000
MPDATA,KXX,1,1, 16, 19, 21, 24, 29,30 !热传导系数THERM COND.(W/M'C)
MPDATA,C,1,1, 500, 540, 560, 590, 600,610 !比热(J/KG'C)
MPDATA,DENS,1,1,7894,7744,7631,7518,7406,7406 !KG/M^3

/COM, 定义PCBN焊接工具的材料参数
MP, EX,2,680E9 !弹性模量(N/M^2)
MP,NUXY,2,0.22 !泊松比
MP, KXX,2,100 !热传导系数(W/M'C)
MP, C,2,750 !比热(J/KG'C)
MP,DENS,2,4280 !密度,KG/M^3


TB,FRIC,3,6,
TBTEMP,25
TBDATA,1,0.4
TBTEMP,200.0
TBDATA,1,0.4
TBTEMP,400.0
TBDATA,1,0.4
TBTEMP,600.0
TBDATA,1,0.3
TBTEMP,800.0
TBDATA,1,0.3
TBTEMP,1000.0
TBDATA,1,0.2


BLOCK,0, W,-L1,L2,0,-T
BLOCK,0,-W,-L1,L2,0,-T

CYLIND,R1,0,0,H,0,360,

LESIZE,17, , ,17,0.16, , , ,1
LESIZE,5, , ,17,6.25, , , ,1

LESIZE,19, , ,17,6.25, , , ,1
LESIZE,7, , ,17,0.16, , , ,1

LESIZE,21, , ,2, , , , ,1
LESIZE,10, , ,2, , , , ,1

LESIZE,8, , ,37,
LESIZE,18, , ,37,

MSHMID,2
VSWEEP,1
VSWEEP,2

MAT,2
ESIZE,0.002
VSWEEP,3


!定义焊接平板之间的接触
!*
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,2,170
ET,3,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
RMORE,10.0
KEYOPT,3,4,3
KEYOPT,3,5,0
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,2
KEYOPT,3,11,0
KEYOPT,3,12,0
KEYOPT,3,2,0
KEYOPT,2,5,0
KEYOPT,3,1,1
! Generate the target surface
ASEL,S,,,5
CM,_TARGET,AREA
TYPE,2
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,,,12
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,

0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
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
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


!定义轴肩底面与焊接工件表面的接触

!*
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MAT,3
MP,EMIS,3,0
R,4
REAL,4
ET,4,170
ET,5,174
R,4,,,1.0,0.1,0,
RMORE,,,500e6,0.0,1.0,
RMORE,0.0,10,0.8,,1,0.95
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,5,4,0
KEYOPT,5,5,3
KEYOPT,5,7,0
KEYOPT,5,8,0
KEYOPT,5,9,1
KEYOPT,5,10,2
KEYOPT,5,11,0
KEYOPT,5,12,0
KEYOPT,5,2,0
KEYOPT,4,5,0
KEYOPT,5,1,1
! Generate the target surface
ASEL,S,,,13
ASEL,A,,,15
ASEL,A,,,16
CM,_TARGET,AREA
TYPE,4
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,,,2
ASEL,A,,,8
CM,_CONTACT,AREA
TYPE,5
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,4
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
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


!定义轴肩的刚性面约束
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
! Define surface-based constraint type of pair
MAT,1
R,5
REAL,5
ET,6,170
ET,7,174
KEYOPT,7,12,5
KEYOPT,7,4,2
KEYOPT,7,2,2
KEYOPT,6,2,1
KEYOPT,6,4,111111
KEYOPT,7,1,1
TYPE,6
! Create a pilot node
! At center of mass of contact geometric entities
N,9793, -8.61577111833e-019,4.49281288333e-019,0.01524
TSHAP,PILO
E,9793
! Generate the contact surface
ASEL,S,,,14
CM,_CONTACT,AREA
TYPE,7
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,6
ESEL,A,TYPE,,7
ESEL,R,REAL,,5
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,6
ESEL,A,TYPE,,7
ESEL,R,REAL,,5
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


/SOLU
ANTYPE,4 !激活瞬态动力学分析
NLGEOM,ON !激活大变形分析
!NROPT,UNSYM !使用非对称的完全牛顿-拉弗森算法


!力学边界条件
NSEL,S,LOC,X,-0.8*W,-W
NSEL,A,LOC,X,0.8*W,W
D,ALL,UZ,0 !DISPLACEMENT CONSTRAINT IN X-DIRECTION
D,ALL,UY,0 !DISPLACEMENT CONSTRAINT IN Y-DIRECTION
D,ALL,UX,0 !DISPLACEMENT CONSTRAINT IN Z-DIRECTION
ALLSEL,ALL


DA,1,UZ,
DA,7,UZ,


!热边界条件
SFA,11,1,CONV,30,25
SFA,8,1,CONV,30,25
SFA,2,1,CONV,30,25
SFA,6,1,CONV,30,25
SFA,4,1,CONV,30,25
SFA,10,1,CONV,30,25
SFA,9,1,CONV,30,25
SFA,3,1,CONV,30,25
SFA,15,1,CONV,30,25
SFA,16,1,CONV,30,25
SFA,14,1,CONV,30,25
SFA,1,1,CONV,300,25
SFA,7,1,CONV,300,25


TREF,25 !REFERENCE TEMPERATURE 25'C
TUNIF,25, !INITIAL CONDITION AT NODES,TEMP 25'C

OUTRES,ALL,ALL
LNSRCH,1
CUTCONTROL,PLSLIMIT,0.15 !

TIMINT,0,STRUC !关闭结构动力学效应
pscontrol,all,on

!*
RESCONTROL,DEFINE,ALL,ALL,0

!LOAD STEP1
D,9793,UZ,-T/4000 !轴肩工具的压入工作面的深度TOOL PLUNGES INTO THE WORKPIECE
D,9793,UX,0
D,9793,UY,0
D,9793,ROTX,0
D,9793,ROTY,0
D,9793,ROTZ,0
TIME,1
AUTOTS,ON !激活自动时间步
KBC,0
NSUBST,10,1000,10
LSWRITE,1


!LOAD STEP2
TIME,6.5
D,9793,ROTZ,11*3.141592 !轴肩工具的转速, 60RPM
DELTIM,0.01,0.004,0.2
LSWRITE,2


!LOAD STEP3
TIME,29
D,9793,ROTZ,45*3.141592 !ROTATION OF TOOL,60RPM
D,9793,UY,60.96E-03 !DISPLACEMENT OF TOOL ALONG WELDLINE
DELTIM,0.01,0.004,0.2
LSWRITE,3

LSSOLVE,1,3,1