| 查看: 4609 | 回复: 3 | ||
teleixi木虫 (小有名气)
|
[求助]
急求ANSYS 旋转高斯体热源加载 或 组合体热源加载方法（APDL流程） 已有1人参与
|
|
最近做激光焊接模拟，在热源在加载时，遇到问题啊，陷入困境，急求 旋转高斯体热源加载 或 组合体热源加载方法（最好能有APDL流程）。非常感谢。 如若能帮忙解决，一定奉上金币。再次表示感谢。  |
回复此楼» 猜你喜欢
TC4钛合金激光焊接金相组织分析
已经有2人回复
-
定向工程博士如何申请
已经有2人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有168人回复
-
紧急招收2026年秋季入学博士生1名（湘潭大学 固体废弃物低碳利用湖南省工程研究中心）
已经有20人回复
-
【急招】“双一流”高校-新能源材料课题组招收2026年秋季入学博士生1名(湘潭大学)
已经有22人回复
-
合金高温性能
已经有0人回复
-
祈祷面上基金能中
已经有129人回复
-
力学研究所某专项团队支撑岗位和特别研究助理岗位招聘启事(长期有效)
已经有10人回复
-
力学研究所某专项团队支撑岗位和特别研究助理岗位招聘启事
已经有10人回复
-
自用分享｜桌面深层文件夹分组映射工具 HyperDesk 免费体验版
已经有2人回复
-
求一篇论文，出自Molecular Simulation，
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 施加高斯热源可以用temperature添加吗 还是就得是heat flux函数
已经有11人回复
- 在施加旋转高斯热源时总提示错误，但感觉程序没问题，不知错在哪？
已经有5人回复
- 【已完结】求助：双椭球热源表格加载中其表格数据怎么来的？
已经有9人回复
- abaqus移动热源加载
已经有4人回复
- 2种焊接移动高斯热源的程序，有什么区别吗
已经有7人回复
- fluent 流体入口的热源加载问题
已经有5人回复
- 求助双椭球热源模型的加载问题
已经有11人回复
- fluent中的体热源加载问题（实际上是热流密度，不是热生成率）
已经有19人回复
- 【求助】fluent 液体仿真的时候热源怎么加载
已经有8人回复
- 【求助】关于高斯热源峰值功率随时间变化的问题
已经有7人回复
【答案】应助回帖
|
/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
3楼2014-10-24 12:06:38
|
本帖内容被屏蔽 |
4楼2020-05-07 11:06:46
