应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 专业学科区 » 机械 » CAD/CAM » 急求ANSYS 旋转高斯体热源加载 或 组合体热源加载方法（APDL流程）
51/1返回列表
查看: 4503  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

teleixi

木虫 (小有名气)

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
激光辅助制造
1楼 2014-04-23 22:26:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
3楼2014-10-24 12:06:38
已阅   申请机械EPI   回复此楼   编辑   查看我的主页
查看全部 4 个回答

jiaka

铁虫 (初入文坛)

【答案】应助回帖

/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
赞一下 回复此楼
不屈不挠，一往无前
2楼2014-05-04 19:34:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wfdyel

禁虫
本帖内容被屏蔽
4楼2020-05-07 11:06:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上和省基金B类撒花 +17 花田半亩～白 2026-04-21 17/850 2026-04-23 07:30 by john198197
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +5 yexuqing 2026-04-19 5/250 2026-04-22 10:01 by easeheart
[论文投稿] 急需审稿人！！！ +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
关闭
请填处理意见
关闭