24小时热门版块排行榜    

查看: 904  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 国王灬咸鱼 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 帮忙看看fileCode +4 wwncly 2026-08-10 7/350 2026-08-10 21:24 by Tide man
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 15/750 2026-08-10 19:29 by 要乐观耀哥
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见