热分析中温度控制部分的命令流怎么实现？ 300

浏览: 2147 回答: 5

其实该问题已经提出过好几次了，邹正刚老师提出了加载方式改为斜坡加载的方式及设置初始子步长、最大和最小子步长，并多次调整步长的大小，以求写入的命令可以达到与真实热电偶控温结果相同的模拟曲线，但通过几次的更改及调试，模具表面的温度仍然达不到稳定状态（即可能存在不保温的情况），下面我将最新的命令流贴出，再次求圈内大神指教（大致问题描述如下：如下图三根加热管给模具加热，顶部的筋条温度由其中一个热电偶监控，要求筋条位置的温度为设定温度的正负5℃ 例如100℃±5℃，命令流如下：

/COM,  Thermal  

/prep7

ET,1,SOLID70

TOFFST,273  

MPTEMP,,,,,,,,  

MPTEMP,1,0  

MPDATA,KXX,1,,56

MPTEMP,,,,,,,,  

MPTEMP,1,0  

MPDATA,C,1,,380

MPTEMP,,,,,,,,  

MPTEMP,1,0  

MPDATA,DENS,1,,7400


ESIZE,0.005,0,  

VMESH,all  


mytime=8000

mystep=20

/SOLU

ANTYPE,4

TRNOPT,FULL

LUMPM,0

OUTRES,ERASE

OUTRES,ALL,ALL 


AUTOTS,1

DELTIM,mystep/5,mystep/5,mystep/5,1

KBC,0

TSRES,ERASE

TUNIF,10,

SFA,all,1,CONV,10,10  


Q1=106157


*do,i,1,mytime,mystep

time,i

SFA,134,1,HFLUX,Q1

SFA,137,1,HFLUX,Q1

SFA,140,1,HFLUX,Q1

allsel

SOLVE


*GET,n4943_temp,NODE,4943,TEMP

*IF,n4943_temp,GT,105,THEN

Q1=0

*ENDIF

*IF,n4943_temp,LT,95,THEN

Q1=106157

*ENDIF

*ENDDO

FINISH

/POST1

PLNSOL,TEMP,0

FINISH

/post26

nsol,2,8219,temp,,

nsol,3,871,temp,,

nsol,4,8205,temp,,

nsol,5,872,temp,,

nsol,6,2385,temp,,

nsol,7,2384,temp,,

nsol,8,2568,temp,,

nsol,9,2569,temp,

Plvar,2,3,4,5,6,7,8,9

/post1


image.png

最终显示的监测节点位置的温度曲线如下：

image.png

理想状态的检测节点温度的曲线如下：

image.png

求大神指导！！！！！（有需要 我可以将模型发出来）

ANSYS 工程热物理及仿真 结构CAE 航空航天 有限元理论及编程

全部回答 (3)

默认 最新
邹正刚 2018年3月13日
根据你的命令流，每一次求解后，会根据 3 个节点的温度来判断是否修改 Q1、Q2、Q3，三个热流值。但是施加热流的命令在这些修改之前。而 ANSYS 中，如果载荷的修改在施加载荷的命令之后，则对于所施加的载荷是不起作用的。因此，你需要在修改热流值之后重新施加热流载荷，新的热流值才能起作用。 有多处需要修改，比如将原来的： *GET,n4943_temp,NODE,4943,TEMP *IF,n4943_temp,GT,105,THEN Q1=0 *ENDIF 需要修改为： *GET,n4943_temp,NODE,4943,TEMP *IF,n4943_temp,GT,105,THEN Q1=0 SFA,134,1,HFLUX,Q1 SFA,137,1,HFLUX,Q1 SFA,140,1,HFLUX,Q1 *ENDIF 以及所有类似情况。
2018年3月13日
评论 1 点赞
李伟 2018年3月6日
不知道是不是因为在循环里面没有确定是后处理器还是求解器；另外一个问题，会不会是因为对流系数太大，降温太快所致 Q1=106157 *do,i,1,mytime,mystep /sol time,i SFA,134,1,HFLUX,Q1 SFA,137,1,HFLUX,Q1 SFA,140,1,HFLUX,Q1 allsel SOLVE /post1 *GET,n4943_temp,NODE,4943,TEMP *IF,n4943_temp,GT,105,THEN Q1=0 *ENDIF *IF,n4943_temp,LT,95,THEN Q1=106157 *ENDIF *ENDDO
2018年3月6日
评论 点赞
技术邻问答 2018年3月6日
@小月 @周游
2018年3月6日
评论 1 点赞