24小时热门版块排行榜    

查看: 827  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 国王灬咸鱼 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 324分 +4 llllkkkhh 2026-03-18 4/200 2026-03-18 14:46 by haxia
[考研] 化工学硕306求调剂 +10 42838695 2026-03-12 10/500 2026-03-18 14:42 by haxia
[考研] 271材料工程求调剂 +3 .6lL 2026-03-18 3/150 2026-03-18 14:36 by laoshidan
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +6 Liwangman 2026-03-15 6/300 2026-03-18 13:21 by 尽舜尧1
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 285化工学硕求调剂(081700) +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
信息提示
请填处理意见