24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sunrui19925

铁虫 (小有名气)

应助: 2 (幼儿园)
金币: 38.1
红花: 5
帖子: 73
在线: 164.6小时
虫号: 3433382
注册: 2014-09-22
性别: MM
专业: 零件成形制造

[求助] ansys热力耦合分析不收敛已有1人参与

最近在用ansys命令流仿真激光加热板材发生弯曲过程，用的间接耦合法先计算出温度场，在进行应力场计算，温度场能收敛，但是应力场不收敛
之前用的求解器是稀疏矩阵法应力开始计算时就不收敛
换了预共轭梯度求解器第一道加热收敛了等第二道加热开始又不收敛了不知道是因为什么原因有没有大神帮忙分析一下是什么原因吗
应力部分命令流如下
!----------应力部分-------------
/prep7
etchg,tts
mp,dens,1,7850
mptemp,1,20,100,300,450,550,600
mptemp,7,720,1000,1500
!mpdata,dens,1,1,7820,7740,7710,7650,7620,7580
!mpdata,dens,1,7,7500
mpdata,ex,1,1,2.07e11,2.03e11,2.0e11,1.62e11,0.98e11,0.60e11
mpdata,ex,1,7,0.40e11,0.10e11,0.10e11
mpdata,nuxy,1,1,0.28,0.31,0.33,0.34,0.36,0.37
mpdata,nuxy,1,7,0.37,0.47,0.49
mpdata,alpx,1,1,1.1e-5,1.1e-5,1.2e-5,1.3e-5,1.4e-5,1.4e-5
mpdata,alpx,1,7,1.4e-5,1.5e-5,1.5e-5

tb,bkin,1,6
tbtemp,20,1
tbdata,1,345.0e6,1.03e10
tbtemp,300,2
tbdata,1,278.0e6,1.00e10
tbtemp,600,3
tbdata,1,127.0e6,0.3e10
tbtemp,800,4
tbdata,1,64e6,0.20e10
tbtemp,1000,5
tbdata,1,11e6,0.05e10
tbtemp,1500,6
tbdata,1,5e6,0.05e10
fini

!----------求解计算-------------
/solu
antype,static
nsel,s,loc,z,0
d,all,all
allsel
tref,20
nlgeom,on
outres,all,all
nsubst,4,10,2
kbc,0
eqslv,pcg,1e-5

*do,tm2,tst,tt+tst,tinc
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

*do,tm2,tt+3*tst,90,1
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

*do,tm2,90,100,2
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

*do,tm2,100+tst,300,10
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

outres,all,all
ldread,temp,,,300,,,rth
time,300
allsel,all
solve

*do,tm2,tst+300,tt+tst+300,tinc
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

*do,tm2,300+tt+2*tst,390,1
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

*do,tm2,390,400,2
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

*do,tm2,400+tst,1200,10
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

*do,tm2,1200,1800,25
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

*do,tm2,1800+tst,4000,100
outres,all,all
ldread,temp,,,tm2,,,rth
time,tm2
allsel,all
solve
*enddo

回复此楼

» 猜你喜欢

青椒八年已不青，大家都被折磨成啥样了？已经有11人回复
限项规定已经有9人回复
免疫学博士有名额，速联系已经有4人回复
交叉科学部支持青年基金，对三无青椒是个机会吗？已经有5人回复
国家基金申请书模板内插入图片不可调整大小? 已经有6人回复
国家级人才课题组招收2026年入学博士已经有5人回复
Fe3O4@SiO2合成已经有6人回复
青年基金C终止已经有4人回复
26申博求博导推荐-遥感图像处理方向已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2016-08-19 22:25:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

奥利奥麦旋风

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 18
帖子: 3
在线: 1.6小时
虫号: 7890988
注册: 2018-02-04
专业: 机械结构强度学

引用回帖:

6楼: Originally posted by sunrui19925 at 2016-08-22 09:44:23
温度梯度是正常的啊第一道加热冷却后再进行的第二道加热两次的温度场是相似的第一道温度场对应的应力场可以算等第二道加热开始时应力计算就不收敛了这是怎么回事呢...

大神们可以帮我看看这个程序流吗
/PREP7
ET, 1, SOLID70
MP, DENS, 1, 7.839e3
MP, EMIS, 1, 0.8
MPTEMP, 1, 20,100,200,300,400,500
MPTEMP, 7, 600, 700, 765, 800, 900, 1000
MPDATA, KXX, 1, 1, 41.84, 40.88, 39.58, 37.69, 35.35, 32.15
MPDATA, KXX, 1, 7, 28.44, 24.61, 26.77, 17.38, 22.07, 24.88
MPTEMP, 1, 20,100,200,300,400,500
MPTEMP,7,600,700,765,800,900,1000
MPDATA,C,1,1,4.67e2,4.91e2,5.27e2,5.65e2,6.11e2,6.78e2
MPDATA,C,1,7,7.87e2,10.50e2,14.23e2,6.72e2,6.34e2,6.49e2
MPTEMP,1,40,60,80,100,120,140
MPTEMP,7,160,180,220,260,300,340
MPTEMP,13,380,480,580,680,780,980
MPTEMP,19,1180,1380,1580
MPDATA,HF,1,1,4.85,6.0,6.71,7.21,7.57,7.92
MPDATA,HF,1,7,8.21,8.45,8.73,8.89,9.91,9.34
MPDATA,HF,1,13,9.49,9.76,9.68,9.7,9.65,9.46
MPDATA,HF,1,19,9.23,9.05,8.87
TOFFST,273
/VIEW,1,1,1,1
BLOCK,0,0.03875,0,0.004,0,0.015
SAVE
Wpoff,0,0,0.012
VSBW,    1
KWPAVE,    2
/DIST, 1, 1.08222638492,1
/REP, FAST
FLST, 5, 6, 4, ORDE, 4
FITEM, 5, 1
FITEM,5,-4
FITEM,5,13
FITEM,5,-14
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,0.000775/2, , , , , , ,1
!*
FLST, 5, 8, 4, ORDE, 4
FITEM, 5, 5
FITEM, 5,-8
FITEM, 5, 21
FITEM, 5,-24
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,0.003875/2, , , , , , ,1
!*
FLST, 5, 2, 6, ORDE, 2
FITEM, 5, 2
FITEM, 5,-3
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1
MSHKEY,0
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
/AUTO, 1
/REP, FAST
SAVE
/SOLU
Q=12984120
L=0.00775
LSIZE=0.1*L
V=1/150
TINC=LSIZE/V
MAX_X=50
MAX_Y=6
MAX_TIME=50
*DIM, FLUX2, TABLE, MAX_X, MAX_Y, MAX_TIME, X, Y, TIME
*DO, K, 1, MAX_TIME, 1
  *DO, I, 1, MAX_X, 1
  FLUX2 (I, 0, K) = (I-1)*LSIZE
  *ENDDO

*DO, J, 1, MAX_Y, 1
  FLUX2 (0, J, K) = (J-1)*LSIZE
  *ENDDO

*DO, I, 1, MAX_X, 1
*DO, J, 1, MAX_Y, 1
   XCENTER=V*(K-1)*TINC
   YCENTER=0.002
   DISTANCE=(I-1)*LSIZE-XCENTER
   *IF,ABS(DISTANCE),GT,L/2,THEN
   FLUX2(I,J,K)=0
   *ELSE
   FLUX2(I,J,K)=Q*(1+2*DISTANCE/L)
   *ENDIF
*ENDDO
*ENDDO

  FLUX2(0,0,K)=(K-1)*TINC
*ENDDO
!*
ANTYPE,TRANS
TUNIF,20
KBC,1
*SET,DELT,1
TIMINT, ON
AUTOTS, OFF
TT=4.650001
*DO, TM, 1E-6, TT, TINC
  FLST, 2, 1, 5, ORDE, 1
  FITEM, 2, 1
  /GO
  !*
  TIME, TM
  AUTOTS, 0
  KBC,0
  !*
  TSRES,ERASE
  FLST,2,1,5,ORDE,1
  FITEM, 2, 1
  /GO
  !*
  !*
  SFA, P51X, 1, HFLUX, %FLUX2%
  /STATUS, SOLU
  SOLVE
  OUTRES, ALL, ALL,
  /PSF, HFLUX, 2
  /REPLOT
*ENDDO
OUTRES, ALL, ALL,
FINISH
/PREP7
ETCHG,TTS
ET, 1, 45
!*
MPTEMP,,,,,,,,
MPTEMP, 1,293,373,473,573,673,773
MPDATA, EX, 1, 1,209E9, 205E9, 200E9, 190E9, 183E9, 174E9
MPTEMP,,,,,,,,
MPTEMP, 1,293,373,473,573,673,773
MPDATA,PRXY,1,1,0.295,0.291,0.295,0.277,0.319,0.336
MPTEMP,,,,,,,,
MPTEMP, 1,373,473,573,673,773,873
MPTEMP, 7,973
MPDATA,ALpX,1,1,1.28E-5,1.34E-5,1.46E-5,1.46E-5,1.47E-5,1.47E-5
MPDATA, ALpX, 1, 7, 1.47E-5
TB, BKIN, 1, 6,,,
TBTEMP, 293
TBDATA, 1, 9.5E8, 8.1E10
TEMP, 373
TBDATA, 1, 7.95E8, 8E10
TBTEMP, 473
TBDATA, 1, 8.35E8, 7.7E10
TBTEMP, 573
TBDATA, 1, 8.15E8, 7.5E10
TBTEMP,673
TBDATA,1,7.65E8,7E10
TEMP,773
TBDATA,1,6.75E8,6.5E10
FINI
/SOLU
LSCLEAR,SOLID
ANTYPE,STATIC,NEW
KBC,1
NSUBST,5,20,3
DELTIM, , , ,3
TIMINT，ON
AUTOTS, ON
L=0.00775
LSIZE=0.1*L
V=1/150
TINC=LSIZE/V
TT=2.90625
*DO, TM, 1E-6, TT, TINC
ASEL, S, LOC, Z, 0.015
LSEL, S, EXT
NSLL, S, L
D, ALL, ALL
ALLSEL
TREF, 293
LDREAD, TEMP,,, TM,,TEMP2,RTH
FINISH
/SOLU
ANTYPE,, REST
TIME, TM
SOLVE
SAVE
*ENDDO
OUTRES, ALL, ALL,
FINISH