| 查看: 829 | 回复: 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 |
回复此楼» 猜你喜欢
重庆科技大学化学化工学院还有8个调剂名额
已经有4人回复
-
请问有做过氯化胆碱/柠檬酸低共熔溶剂的大佬吗？
已经有1人回复
-
无机非金属材料论文润色/翻译怎么收费?
已经有255人回复
-
2025基金何时放榜？是八月多少？
已经有13人回复
-
macromolecules投稿
已经有0人回复
-
求In2S3（JCPDS card number 32-0456）的cif文件
已经有4人回复
-
成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。
已经有19人回复
-
晶圆激光切割保护液
已经有0人回复
zhukingsly
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 233.1
- 红花: 1
- 帖子: 140
- 在线: 129.5小时
- 虫号: 1031791
- 注册: 2010-05-30
- 性别: GG
- 专业: 聚合物共混与复合材料
2楼2016-04-20 21:18:16
5