| 查看: 4653 | 回复: 3 | |||
teleixi木虫 (小有名气)
|
[求助]
急求ANSYS 旋转高斯体热源加载 或 组合体热源加载方法（APDL流程） 已有1人参与
|
|
最近做激光焊接模拟，在热源在加载时，遇到问题啊，陷入困境，急求 旋转高斯体热源加载 或 组合体热源加载方法（最好能有APDL流程）。非常感谢。 如若能帮忙解决，一定奉上金币。再次表示感谢。  |
回复此楼» 猜你喜欢
航天医学动物模型构建，鼠尾吊失重模型实验箱如何评估心血管系统响应？
已经有0人回复
-
焦虑动物模型如何评估？多通道大小鼠饮水冲突分析系统量化冲突下行为决策
已经有0人回复
-
金属材料论文润色/翻译怎么收费?
已经有115人回复
-
ZL-AI-BT 大小鼠步态精细行为分析系统怎样提取多维步态参数？
已经有0人回复
-
ZL-AI-NORAI 新物体识别精细行为分析系统如何自动判别探索行为？
已经有0人回复
-
ZL-AI-Y AI 大小鼠 Y 迷宫精细行为分析系统怎样自动计算自发交替率、臂区停留时长？
已经有0人回复
-
小动物主动跑轮系统如何捕捉跑轮方向、转数等运动数据？
已经有0人回复
-
英国埃克塞特大学（QS 136）招CSC博士生，3D打印增材制造方向
已经有0人回复
-
曼彻斯特大学2027 中英CSC博士全额奖学金申报-智能制造/金属增材/焊接
已经有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
200