24小时热门版块排行榜    

查看: 828  |  回复: 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhukingsly

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
请问下,你这个问题解决了吗?
2楼2016-04-20 21:18:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 国王灬咸鱼 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 295求调剂 +3 一志愿京区211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 9/450 2026-03-18 15:05 by stone_128
[考研] 314求调剂 +8 无懈可击的巨人 2026-03-12 8/400 2026-03-18 14:50 by haxia
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 材料专硕274一志愿陕西师范大学求调剂 +6 薛云鹏 2026-03-13 6/300 2026-03-18 14:14 by 脱颖而出
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 265求调剂 +3 梁梁校校 2026-03-17 3/150 2026-03-18 09:12 by zhukairuo
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考博] 26申博 +4 八旬速览 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
信息提示
请填处理意见