ANSYS 14.0 实例 搅拌摩擦焊接(焊接)模拟
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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
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
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
! 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
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
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
!定义轴肩的刚性面约束
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
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
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
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
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
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
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
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)
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
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,
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)
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
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
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,
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
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
H=15.24E-03 !轴肩的高度H
L1=R1 !轴肩中心线在焊接线上的起始位置
L2=L-L1
NR1=3.141593*11 !第二载荷步的转速
NR2=3.141593*45 !第三载荷步的转速
UY1=60.96E-03 !轴肩沿焊接线的平动距离
/FILNAME,Friction Stir Welding,0
/PREP7
/COM, ***** 模型参数 ********
L=76.2E-03 !每一个平板的长度
W=31.75E-03 !每一个平板的宽度
T=3.18E-03 !每一个平板的厚度
R1=7.62E-03 !焊机工具的轴肩半径
/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
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
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
RMORE,,,,,,
RMORE,,,,,1000
!定义轴肩底面与焊接工件表面的接触
!*
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
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
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
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
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
! 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
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
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
!定义轴肩的刚性面约束
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
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
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
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
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
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
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
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)
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
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,
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)
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
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
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,
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
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
H=15.24E-03 !轴肩的高度H
L1=R1 !轴肩中心线在焊接线上的起始位置
L2=L-L1
NR1=3.141593*11 !第二载荷步的转速
NR2=3.141593*45 !第三载荷步的转速
UY1=60.96E-03 !轴肩沿焊接线的平动距离
/FILNAME,Friction Stir Welding,0
/PREP7
/COM, ***** 模型参数 ********
L=76.2E-03 !每一个平板的长度
W=31.75E-03 !每一个平板的宽度
T=3.18E-03 !每一个平板的厚度
R1=7.62E-03 !焊机工具的轴肩半径
/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
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
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
RMORE,,,,,,
RMORE,,,,,1000
!定义轴肩底面与焊接工件表面的接触
!*
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
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