| 查看: 905 | 回复: 1 | |||
[交流]
我仿真一个打印过程的温度变化，但是一直出错，请各位大神帮忙指点一下 已有1人参与
|
|
/BATCH FINISH /CLEAR,START /CWD,G:\Exercise /FILNAME,FDMTEMPERTURE,0 /TITLE,ANSYS OF FDM TEMPERTURE *SET,L1,1E-2 *SET,W1,1.2E-3 *SET,H1,0.3E-3 !设置模型尺寸为10mm*1.2mm*0.3mm !设置网格尺寸 *SET,L2,L1/25 *SET,W2,0.4E-3 *SET,H2,H1/3 !设置单元尺寸为0.4mm*0.4mm*0.1mm（喷嘴直径为0.4，分层厚度0.1mm） ARG1=L1/L2*W1/W2*H1/H2 !计算单元总数量 /CONFIG,NRES,10000 !************************************* !进入前处理 !************************************* /PREP7 ET,1,SOLID70 MP,DENS,1,1150 !设置密度 MPTEMP,1,50,100,150,200 !指定四个温度 MPDATA,C,1,1,1470,1490,1710,2020 !比热容 MPDATA,KXX,1,1,3.0E-2,2.8E-2,2.9E-2,3.3E-2 !导热系数 MPDATA,EX,1,1,3.5E5,2.48E9,1.68E9,0.5E9 !弹性模量 MPDATA,NUXY,1,1,0.38,0.39,0.40,0.41 !泊松比 MPDATA,ALPX,1,1,8.51E-5,8.42E-5,8.40E-5,8.38E-5 !热膨胀系数 BLOCK,0,L1,0,W1,0,H1 !建立成型件的20mm*1mm*0.3mm的模型 !网格划分 /PNUM,AREA,1 !显示面编号 VATT,1,,0 !体单元类型solid70 /PNUM,LINE,1 MSHAPE,0,3D !单元形状六面体 MSHKEY,1 !映射网格划分 LESIZE,1,W2 LESIZE,4,L2 LESIZE,9,H2 !单元尺寸大小 VMESH,1 EPLOT /VIEW,,-1,-2,-3 /IMAGE,SAVE,VMESH,BMP FINISH !************************************** !进入载荷加载和求解 !************************************* /SOLU ANTYPE,TRANSIENT,NEW AR21=0 !每一个载荷步加载时间是0.02s，单元激活时间间隔也是0.02s，扫描速度20mm/s AR22=0.02 TIME,AR21 AUTOTS,ON !自动时间步选项（非线性分析打开加速收敛） NROPT,FULL,,ON ! 牛顿-拉普森选项 TRNOPT,FULL !瞬态分析选项，指定最大分析模态数 LUMPM,OFF !使用一致质量矩阵 KBC,1 DELTIM,AR21/4,AR21/4,AR21/2,0 !载荷子步设定 TIMINT,ON MAX_TEMP=220 !设定温度值220摄氏度 NSEL,S,LOC,X,0,L1,L2 SF,ALL,CONV,72,25 NSEL,S,LOC,Y,0,W1,W2 SF,ALL,CONV,72,25 NSEL,S,LOC,Z,0,H1,H2 SF,ALL,CONV,72,25 !初始化所有节点，加载对流 /REPLOT,FAST OUTRES,ALL,ALL OUTPR,NSOL,ALL !生死单元循环命令 ESEL,S,ELEM,,1,ARG1,1,0 EKILL,ALL ALLSEL,ALL *DO,AR25,1,H1/H2,1 *DO,AR23,1,L1/L2,1 ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0 !依次选择下一个单元激活 EALIVE,ALL *DO,J,1,8,1 D,NELEM((AR25-1)*W1/W2*L1/L2+AR23,J),TEMP,MAX_TEMP !给每个单元的八个节点施加恒定温度 *ENDDO ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0 NSLE,S,ALL !选择单元中的所有节点 SF,ALL,CONV,72,25 ALLSEL,ALL SOLVE ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0 *DO,J,1,8,1 DDELE,NELEM((AR25-1)*W1/W2*L1/L2+AR23,J),TEMP !删除每个单元的八个节点施加恒定温度 *ENDDO ALLSEL,ALL AR22=AR22+AR21 !*IF,AR22,GT,AR21*ARG1,*EXIT FINISH ANTYPE,TRANSIENT,RESTART ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0 *DO,J,1,8,1 DDELE,NELEM((AR25-1)*W1/W2*L1/L2+AR23,J),TEMP !删除每个单元的八个节点施加恒定温度 *ENDDO TIME,AR22 *ENDDO *DO,AR26,L1/L2,1,-1 ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+L1/L2+AR26,,,0 EALIVE,ALL *ENDDO *DO,AR27,L1/L2,1,-1 ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0 EALIVE,ALL *DO,J,1,8,1 D,NELEM((AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,J),TEMP,MAX_TEMP !给每个单元的八个节点施加恒定温度 *ENDDO ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0 NSLE,S,ALL !选择单元中的所有节点 SF,ALL,CONV,72,25 ALLSEL,ALL SOLVE ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0 *DO,J,1,8,1 DDELE,NELEM((AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,J),TEMP !删除每个单元的八个节点施加恒定温度 *ENDDO ALLSEL,ALL AR22=AR22+AR21 !*IF,AR22,GT,AR21*ARG1,*EXIT FINISH ANTYPE,TRANSIENT,RESTART ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0 *DO,J,1,8,1 DDELE,NELEM((AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,J),TEMP !删除每个单元的八个节点施加恒定温度 *ENDDO TIME,AR22 *ENDDO *ENDDO SAVE /POST1 PLNSOL,TEMP,,1 ANDATA,0.02,1,2,0,0,1,0,0 !生成某一范围内的结果数据的顺序梯度线动画 FINISH |
回复此楼» 猜你喜欢
材料/力学/岩土方向，分享三套顶刊复现实战（PINN+机器学习+生成式AI）
已经有0人回复
-
如何用原子层沉积（ALD)高温固体源
已经有0人回复
-
金属材料论文润色/翻译怎么收费?
已经有253人回复
-
ai+cfd智能流体工程实战应用
已经有10人回复
-
树脂浇铸体出现纹理问题是什么原因？
已经有5人回复
-
ai+cfd智能流体工程前沿技术
已经有5人回复
-
ai+cfd：fluent、pinn、python、pyfluent智能流体技术前沿
已经有1人回复
-
工作党想考在职研究生
已经有21人回复
zhukingsly
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 234.6
- 红花: 1
- 帖子: 140
- 在线: 129.5小时
- 虫号: 1031791
- 注册: 2010-05-30
- 性别: GG
- 专业: 聚合物共混与复合材料
2楼2016-04-20 21:18:16
100