| 查看: 4481 | 回复: 3 | ||
teleixi木虫 (小有名气)
|
[求助]
急求ANSYS 旋转高斯体热源加载 或 组合体热源加载方法（APDL流程） 已有1人参与
|
|
最近做激光焊接模拟，在热源在加载时，遇到问题啊，陷入困境，急求 旋转高斯体热源加载 或 组合体热源加载方法（最好能有APDL流程）。非常感谢。 如若能帮忙解决，一定奉上金币。再次表示感谢。  |
回复此楼» 猜你喜欢
智能制造和机器人国家科技重大专项研究团队招收调剂研究生10人
已经有4人回复
-
燕山大学先进结构材料团队招收金属材料相关方向 2026年硕士研究生
已经有0人回复
-
机械工程论文润色/翻译怎么收费?
已经有287人回复
-
青岛大学纺织服装学院-极端环境纺织材料与制品团队-招2026年考研调剂
已经有0人回复
-
求调剂建议
已经有0人回复
-
国产沪模心肺复苏：从厂家到品牌，深度剖析
已经有0人回复
-
高温合金厂家
已经有0人回复
-
085500机械 335求调剂
已经有0人回复
-
动力工程求调剂
已经有0人回复
-
工科370求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 施加高斯热源可以用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
