24小时热门版块排行榜

返回列表

teleixi

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 3035.7
散金: 30
帖子: 252
在线: 120.4小时
虫号: 2224038
注册: 2013-01-05
性别: GG
专业: 零件成形制造

[求助] 急求ANSYS 旋转高斯体热源加载或组合体热源加载方法（APDL流程）已有1人参与

最近做激光焊接模拟，在热源在加载时，遇到问题啊，陷入困境，急求旋转高斯体热源加载或组合体热源加载方法（最好能有APDL流程）。非常感谢。
如若能帮忙解决，一定奉上金币。再次表示感谢。

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

施加高斯热源可以用temperature添加吗还是就得是heat flux函数已经有11人回复
在施加旋转高斯热源时总提示错误，但感觉程序没问题，不知错在哪？已经有5人回复
【已完结】求助：双椭球热源表格加载中其表格数据怎么来的？已经有9人回复
abaqus移动热源加载已经有4人回复
2种焊接移动高斯热源的程序，有什么区别吗已经有7人回复
fluent 流体入口的热源加载问题已经有5人回复
求助双椭球热源模型的加载问题已经有11人回复
fluent中的体热源加载问题（实际上是热流密度，不是热生成率）已经有19人回复
【求助】fluent 液体仿真的时候热源怎么加载已经有8人回复
【求助】关于高斯热源峰值功率随时间变化的问题已经有7人回复

激光辅助制造

1楼 2014-04-23 22:26:59

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

jiaka

铁虫 (初入文坛)

应助: 3 (幼儿园)
金币: 30.5
帖子: 17
在线: 2.9小时
虫号: 1807844
注册: 2012-05-10
性别: GG
专业: 计算机软件

【答案】应助回帖

/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