| 查看: 2954 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
ANSYS热分析时热量传递不到相邻的单元内，求教大家是什么原因 已有1人参与
|
|||
仿真条件如下图1所示，模型为长方体，单元类型为Solid 70，初始温度为172，在上表面施加移动的热源，对模型进行加热，分析其温度场，计算结果如图2所示 图1  图2 发现沿着深度方向，热量只能从一个单元的上层节点传递到本单元下层节点，但从该单元的下层节点就无法传递到与其相邻的其他单元了 对模型用不同的网格划分，其他条件不变，各次仿真温度沿深度方向的变化如下所示 0.1mm网格划分：温度升高刚好从0.1mm处开始，0.01mm之后，也就是除上表面的第一层单元以外的单元，温度均是初始温度 而在0到0.01mm之间，也就是在直接被热源作用的单元间，温度则是均匀地从最高温度降低到初始温度（172度），说明温度没有传导到下一层的单元  按其他网格大小进行划分，情况均类似 0.05mm网格划分：温度升高刚好从0.05mm处开始  0.025mm网格划分：温度升高刚好从0.025mm处开始  温度为什么会传导不到下一层单元？卡在这个问题上好几周了，各位大侠求指点迷津，感激不尽~ |
回复此楼» 猜你喜欢
求博导｜生物质基多孔碳/超级电容方向，已有相关成果，寻能源材料/碳材料方向老师
已经有3人回复
-
二苯甲酮酸类衍生物
已经有6人回复
-
接受任何调剂
已经有4人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有12人回复
-
320求调剂
已经有4人回复
-
一志愿华中农业071010，320求调剂
已经有18人回复
-
295分求调剂
已经有5人回复
-
294求调剂
已经有8人回复
-
山东双非院校考核超级无底线，领导幸灾乐祸，教师遭殃恐
已经有8人回复
-
0854求调剂
已经有25人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ansys热分析求助
已经有2人回复
- ANSYS 热分析中默认电压是多少
已经有6人回复
- ansys热分析建模求助
已经有0人回复
- ansys热分析
已经有2人回复
- ansys瞬态热分析
已经有11人回复
- ANSYS热分析求助
已经有8人回复
- ansys热分析中的问题
已经有4人回复
- ansys热应力分析
已经有4人回复
- ANSYS 热分析时热源无法施加
已经有0人回复
- ANSYS热分析
已经有4人回复
- ansys暂态热分析
已经有8人回复
- ansys瞬态热分析
已经有9人回复
- ANSYS热分析相关
已经有5人回复
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
6楼2017-04-11 07:59:50
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
2楼2017-04-07 08:35:45
|
不知道怎么上传模型  ，我把我的命令流贴上来吧，命令流文件在最后可以直接下载非常感谢~~ ********************************************************** finish /clear /filname,esize_debug /title,selective laser sintering /prep7 /unit,si et,1,solid70 !****材料属性设置**** MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,445 !粉末密度445 kg/m3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,2660 ! 比热容：2660J/（kg·K） MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,, 0.111 ! 热导率：0.111W/（m·K） !******************* !****模型几何参数设置**** length=5e-3 !长5mm width=1e-3 !宽1mm height=0.5e-3 !高0.5mm !********************** block,0, length,0, width,0, height esize,0.1e-3 vmesh,all !****在上表面建立表面效应单元，用于加载对流载荷（GUI操作）****** !* ET,2,SURF152 !* KEYOPT,2,1,0 KEYOPT,2,2,0 KEYOPT,2,3,0 KEYOPT,2,4,1 KEYOPT,2,5,0 KEYOPT,2,6,0 KEYOPT,2,7,0 KEYOPT,2,8,4 KEYOPT,2,9,0 !* !* TYPE, 2 MAT, 1 REAL, ESYS, 0 SECNUM, TSHAP,LINE !* ASEL,S, , , 2 NSLA,S,1 !* ESURF,0 !* ESEL,S,TYPE,,2 !* SFE,ALL,1,CONV, ,25 SFE,ALL,1,CONV,2,172 !* !************************* ALLSEL,ALL TUNIF,172 FINISH /SOL PI=3.1415926 !****激光参数设置**** P=12 !功率为12W w0=0.15e-3 !光斑半径为0.15mm v=3.5 !激光扫描速度为3.5m/s d=0.12e-3 !激光扫描间距为0.12mm !********************** ! ****瞬态求解参数设置**** lstep=0.1e-3 !扫描步长为0.1mm time_step=lstep/v !载荷步时间大小 ANTYPE,4 !定义新分析 TRNOPT,FULL TOFFST,273, !定义摄氏度的温度偏移量 !* TIMINT,1 !打开时间积分选项进行瞬态分析 TINTP, ,0.25250625,0.505, , , , !这条语句是GUI附送的，不知道有没有问题 !* !************************ OUTRES,ERASE OUTRES,ALL,1 my_time=0 !定义一个时间变量 XX=0.5e-3 !XX,YY用于定义高斯光束能量公式中的局部坐标的原点 YY=0 !**** 高斯光束能量公式定义（GUI方式自动生成）**** !****GAUSS=2*P/(PI*(w0)^2)*exp(-2*({X}^2+{Y}^2)/(w0)^2)**** LOCAL,11,0,XX,YY,height !定义局部坐标系11 *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNCCSYS *SET,_FNCNAME,'GAUSS' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *SET,_FNC_C1(1),P ! *SET,_FNC_C2(1),w0 ! *SET,_FNCCSYS,11 !将方程坐标系设定为自己定义的局部坐标系 ! /INPUT,GAUSS.func,,,1 *DIM,%_FNCNAME%,TABLE,6,24,1,,,,%_FNCCSYS% ! ! Begin of equation: 2*P/(PI*(w0)^2)*exp(-2*({X}^2+{Y}^2)/(w0)^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), 0.0 *SET,%_FNCNAME%(6,0,1), 0.0 *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 17 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 17 *SET,%_FNCNAME%(0,3,1), 0, -1, 0, 2, 0, 0, 18 *SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, 18, 17, -1 *SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -3 *SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 3, -3 *SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, -2, 4, -4 *SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 0, 0, -2 *SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, -2, 2, -3 *SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, -4 *SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -4, 3, -2 *SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, 2, 17, -2 *SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 2, 0, 0, 3 *SET,%_FNCNAME%(0,16,1), 0.0, -5, 0, 1, 3, 17, -2 *SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 1, -5 *SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, -3, 3, -2 *SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, 18 *SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, 18, 17, -2 *SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -4, 4, -3 *SET,%_FNCNAME%(0,22,1), 0.0, -2, 7, 1, -2, 0, 0 *SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,24,1), 0.0, 99, 0, 1, -3, 0, 0 ! End of equation: 2*P/(PI*(w0)^2)*exp(-2*({X}^2+{Y}^2)/(w0)^2) !--> !************************************************* !******高斯热源沿x方向移动并求解的程序*************** *IF,XX,EQ, 0.5e-3,THEN *DO,i, 0,10,1 !每循环一次高斯热源沿x正向移动lstep的距离 LOCAL,11,0,XX,YY,height !****每个载荷步的求解参数设置**** my_time=my_time+time_step TIME,my_time AUTOTS,1 DELTIM,time_step/100, time_step/200, time_step/10,1 KBC,1 !* TSRES,ERASE !* !******************************** !****加载**** FLST,2,1,5,ORDE,1 FITEM,2,2 /GO !* SFA,P51X,1,HFLUX,%GAUSS% !************* !****求解，完成后删除载荷**** ALLSEL,ALL !save SOLVE SFADELE,2,1,HFLUX !**************************** XX=XX+lstep *ENDDO !****显示结果**** FINISH /POST1 !* /EFACET,1 PLNSOL, TEMP,, 0 !******************** |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : APDL.doc
2017-04-07 10:30:13, 43 K
3楼2017-04-07 10:32:58
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
luoyehsw4楼2017-04-08 08:20:00
20