在用ansys电磁感应淬火时遇到的马氏体分布求解问题

lijishan

首先计算感应加热7秒后的温度,然后水冷18秒。然后通过计算冷却速度计算马氏体分布。（生成马氏体的判定条件为马氏体转变的最低温度为800℃,冷却速度不小于200）。但是在计算淬硬层分布时出现了如下的错误。
cumulative iteration 116 may have been solved using different model or boundary condition than currently stored.post1 results may be erroneous unless you resume from a jobname.db file for this substep. 请高手帮忙,不胜感谢。
命令流如下
/BATCH
/FILNAME,2D-1-4.5E9,1
/TITLE,2D-1-4.5E9
/PREP7
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7830
MPTEMP,,,,,,,,
MPTEMP,1,25
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,600
MPTEMP,5,800
MPTEMP,6,950
MPDATA,C,1,,472
MPDATA,C,1,,521
MPDATA,C,1,,607
MPDATA,C,1,,745
MPDATA,C,1,,796
MPDATA,C,1,,677
MPTEMP,,,,,,,,
MPTEMP,1,25
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,600
MPTEMP,5,800
MPTEMP,6,950
MPDE,EX,1
MPDE,PRXY,1
MPDATA,EX,1,,2.06e11
MPDATA,EX,1,,1.92e11
MPDATA,EX,1,,1.75e11
MPDATA,EX,1,,1.53e11
MPDATA,EX,1,,1.25e11
MPDATA,EX,1,,1.11e11
MPDATA,PRXY,1,,0.28
MPDATA,PRXY,1,,.28
MPDATA,PRXY,1,,.28
MPDATA,PRXY,1,,.28
MPDATA,PRXY,1,,.28
MPDATA,PRXY,1,,.28
MPTEMP,1,25
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,600
MPTEMP,5,800
MPTEMP,6,950
UIMP,1,REFT,,,
MPDATA,ALPX,1,,2.5e-6
MPDATA,ALPX,1,,2e-5
MPDATA,ALPX,1,,4.2e-5
MPDATA,ALPX,1,,6.25e-5
MPDATA,ALPX,1,,8.5e-5
MPDATA,ALPX,1,,0.0001
MPTEMP,,,,,,,,
MPTEMP,1,25
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,600
MPTEMP,5,800
MPTEMP,6,950
MPDE,KXX,1
MPDATA,KXX,1,,49.8
MPDATA,KXX,1,,44.8
MPDATA,KXX,1,,39.4
MPDATA,KXX,1,,34
MPDATA,KXX,1,,29
MPDATA,KXX,1,,25.5
TB,BKIN,1,6,2,1
TBTEMP,25
TBDATA,,7e8,2.06e10,,,,
TBTEMP,200
TBDATA,,6.3e8,1.98e10,,,,
TBTEMP,400
TBDATA,,5.6e8,1.83e10,,,,
TBTEMP,600
TBDATA,,4.2e8,1.56e10,,,,
TBTEMP,800
TBDATA,,3.6e8,1.12e10,,,,
TBTEMP,950
TBDATA,,3.3e8,9e9,,,,
*SET,D,1
*SET,DIV,5*D
ET,1,PLANE13
KEYOPT,1,1,4
KEYOPT,1,3,1
TOFFST,273
K,1,,,,
K,2,0.045,,,
K,3,0.045,0.034,,
K,4,0.048,0.034,,
K,5,0.048,0.037,,
K,6,0.052,0.037,,
K,7,0.052,0.046,,
K,8,0,0.046,,
K,9,0.036,0,,
K,10,0.045-0.001*D,0,,
K,11,0.045-0.001*D,0.034,,
K,12,0.048,0.037+0.001*D,,
K,13,0.052,0.037+0.001*D,,
K,14,0.048,0.046,,
K,15,0.036,0.046,,
K,16,0,0.034,,
K,17,0.036,0.034,,
/PNUM,KP,1
/PNUM,LINE,1
/PNUM,AREA,1
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0
/PNUM,ELEM,0
LSTR,9,10
LSTR,10,2
LSTR,10,11
LSTR,17,11
LSTR,11,3
LARC,11,12,4,0.003+0.001*D,
LSTR,5,12
LSTR,12,14
LSTR,12,13
LSTR,6,13
LSTR,13,7
LSTR,1,9
LSTR,2,3
LARC,3,5,4,0.003,
LSTR,5,6
LSTR,7,14
LSTR,14,15
LSTR,15,8
LSTR,8,16
LSTR,16,1
LSTR,9,17
LSTR,17,15
LSTR,17,16
AL,2,3,5,13
AL,5,6,7,14
AL,7,9,10,15
AL,1,3,4,21
AL,4,6,8,17,22
AL,8,9,11,16
AL,12,20,21,23
AL,18,19,22,23
NUMMRG,ALL
NUMCMP,ALL
LSEL,S,LINE,,2
LESIZE,ALL,,,DIV,,1
LSEL,S,LINE,,5
LESIZE,ALL,,,DIV,,1
LSEL,S,LINE,,7
LESIZE,ALL,,,DIV,,1
LSEL,S,LINE,,10
LESIZE,ALL,,,DIV,,1
LSEL,S,LINE,,3
LESIZE,ALL,,,30,,1
LSEL,S,LINE,,13
LESIZE,ALL,,,30,,1
LSEL,S,LINE,,20
LESIZE,ALL,,,30,,1
LSEL,S,LINE,,21
LESIZE,ALL,,,30,,1
LSEL,S,LINE,,12
LESIZE,ALL,,,30,,1
LSEL,S,LINE,,18
LESIZE,ALL,,,30,,1
LSEL,S,LINE,,23
LESIZE,ALL,,,30,,1
LSEL,S,LINE,,17
LESIZE,ALL,,,12,,1
LSEL,S,LINE,,19
LESIZE,ALL,,,12,,1
LSEL,S,LINE,,22
LESIZE,ALL,,,12,,1
LSEL,S,LINE,,1
LESIZE,ALL,,,20,0.25,1
LSEL,S,LINE,,4
LESIZE,ALL,,,20,0.25,1
LSEL,S,LINE,,8
LESIZE,ALL,,,20,4,1
LSEL,S,LINE,,11
LESIZE,ALL,,,20,4,1
LSEL,S,LINE,,9
LESIZE,ALL,,,10,2,1
LSEL,S,LINE,,15
LESIZE,ALL,,,10,2,1
LSEL,S,LINE,,16
LESIZE,ALL,,,10,0.5,1
LSEL,S,LINE,,6
LESIZE,ALL,,,24,,1
LSEL,S,LINE,,14
LESIZE,ALL,,,24,,1
AMESH,1
AMESH,2
AMESH,3
AMESH,4
LSEL,S,LINE,,22
LSEL,A,LINE,,17
LCCAT,22,17
AMESH,5
LDELE,24,,,0
AMESH,6
AMESH,7
AMESH,8
ALLSEL,ALL
NUMMRG,ALL
NUMCMP,ALL
/AUTO,1
/REP,FAST
GPLOT
FINISH
/SOLU
ANTYPE,TRANS
TRNOPT,FULL
NLGEOM,OFF
TIME,7
AUTOTS,OFF
NSUBST,70,,,ON
OUTRES,ALL,ALL
KBC,0
ALLSEL,ALL
TUNIF,25
ASEL,S,AREA,,1!选择面1,2,3
ASEL,A,AREA,,2
ASEL,A,AREA,,3
APLOT
NSLA,S,1
BF,ALL,HGEN,4.5E9
ALLSEL,ALL
LSEL,S,LINE,,10
LSEL,A,LINE,,11
LSEL,A,LINE,,13
LSEL,A,LINE,,14
LSEL,A,LINE,,15
SFl,ALL,CONV,100,,25
lplot
ALLSEL,ALL
LSEL,S,LINE,,1
LSEL,A,LINE,,2
LSEL,A,LINE,,12
LPLOT
NSLL,S,1
D,ALL,UY,0
ALLSEL,ALL
LSEL,S,LINE,,19
LSEL,A,LINE,,20
LPLOT
NSLL,S,1
D,ALL,UX,0
ALLSEL,ALL
GPLOT
LSWRITE,1
TIME,25
AUTOTS,OFF
NSUBST,180,,,ON
OUTRES,ALL,ALL
ALLSEL,ALL
ASEL,S,AREA,,1
ASEL,A,AREA,,2
ASEL,A,AREA,,3
APLOT
NSLA,S,1
BFDELE,ALL,HGEN
ALLSEL,ALL
LSEL,S,LINE,,10
LSEL,A,LINE,,11
LSEL,A,LINE,,13
LSEL,A,LINE,,14
LSEL,A,LINE,,15
LPLOT
NSLL,S,1
SFDELE,ALL,CONV
allsel,all
LSEL,S,LINE,,10
LSEL,A,LINE,,11
LSEL,A,LINE,,13
LSEL,A,LINE,,14
LSEL,A,LINE,,15
SFl,ALL,CONV,15000,,25
ALLSEL,ALL
GPLOT
LSWRITE,2
LSSOLVE,1,2,1
FINISH
save
!以下为计算马氏体分布的程序
/PREP7

*CSET,1,2,1,3,


!这个可能有问题
MPCOPY, ,1,3

TBCOPY,ALL,1,3

!此处拟定3号材料参数为马氏体材料参数
FINISH
/POST1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!


说明


!!!!!!!!
!!!!!!


在使用本程序前需对下列参数进行设置

!!!!!!!!
!!!!!!
ETEMP、EDERI、NOFLS、TOFFST
!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*SET, ETEMP, 800
!马氏体转变温度底限
*SET, EDERI, 200
!马氏体转变冷却速度底限
*SET, NOFLS, 180
!冷却载荷步的子步数
*SET, TOFFST, 7

!加热载荷步停止时间
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL, ALL
*GET, A, ELEM, , COUNT

!计算出所选单元的总个数
*SET, ENUMBER, A

!计算模型单元总数
*DIM, TANDERI1, , ENUMBER, 3

!临时马氏体转变判断结果矩阵
*DIM, EDERIV,,ENUMBER
!各单元某时刻冷速矩阵
SET, 1, LAST, 1

!读入载荷步1最后一个子步

ETABLE, ELETEMP, TEMP,
!定义单元温度表

*VGET,T1,ELEM,,ETAB,ELETEMP,,2
!创建单元温度参数ELETEMP
*DO,J,1,NOFLS-1,1


!读入载荷步2内的所有子步

SET,2,J,1

ETABLE,ELETEMP,TEMP,
!定义单元温度表

*VGET,T2,ELEM,,ETAB,ELETEMP,,2
!创建单元温度参数ELETEMP

SET,2,J+1,1


ETABLE,ELETEMP,TEMP,

!定义单元温度表

*VGET,T3,ELEM,,ETAB,ELETEMP,,2 !创建单元温度参数ELETEMP

*IF,J,EQ,1,THEN



*DO,K,1,ENUMBER,1

!计算每个单元温度变化的微分

EDERIV(K)=10*(T2(K)-T1(K))


*ENDDO

!结束计算每个单元温度的微分


*DO,K,1,ENUMBER,1

!存储分析结果到TANDERI1矩阵


TANDERI1(K,1)=TOFFST+0.1*J
!马氏体组织转变时间


TANDERI1(K,2)=T2(K)

TANDERI1(K,3)=EDERIV(K)



*ENDDO

!结束归纳存储分析结果

*ELSEIF,J,GT,1,AND,J,LT,NOFLS-1


*DO,K,1,ENUMBER,1

!计算每个单元温度的微分

EDERIV(K)=5*(T3(K)-T1(K))


*ENDDO


!结束计算每个单元温度的微分

*DO,K,1,ENUMBER,1


!存储分析结果到TANDERI1


TANDERI1(K,1)=TOFFST+0.1*J
!马氏体组织转变时间

TANDERI1(K,2)=T2(K)

TANDERI1(K,3)=EDERIV(K)



*ENDDO

!结束存储分析结果到TANDERI1

*ELSEIF,J,EQ,NOFLS-1


*DO,K,1,ENUMBER,1

!计算每个单元温度的微分

EDERIV(K)=10*(T3(K)-T2(K))


*ENDDO

!结束计算每个单元温度变化的微分


*DO,K,1,ENUMBER,1
!存储分析结果到TANDERI1


TANDERI1(K,1)=TOFFST+0.1*J
!马氏体组织转变时间


TANDERI1(K,2)=T2(K)

TANDERI1(K,3)=EDERIV(K)



*ENDDO

!结束存储分析结果到TANDERI1

*ENDIF

*DO,K,1,ENUMBER,1


T1(K)=T2(K)
!转移温度矩阵准备下次计算

T2(K)=T3(K)
!进入单元温度结果数据处理阶段

*IF,TANDERI1(K,2),ABGT,ETEMP,AND,TANDERI1(K,3),ABGT,EDERI,THEN

/PREP7

ESEL,S,,,K

MPCHG,3,ALL

!材料属性转变为马氏体材料

ESEL,ALL

FINISH

/POST1

*ENDIF

*ENDDO

!结束结果数据处理阶段
*ENDDO
!结束读入载荷步2内的所有子步
FINISH
/PREP7

/PNUM,MAT,1
/NUM,1
ESEL,ALL
EPLOT
/AUTO,1
/REP,FAST

FINISH

xiaolaohu01

很强大啊