应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » ansys热力耦合分析不收敛
121/2返回列表12下一页
查看: 3507  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

sunrui19925

铁虫 (小有名气)

[求助] 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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-08-19 22:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

sunrui19925

铁虫 (小有名气)
引用回帖:
5楼: Originally posted by shgao20 at 2016-08-22 08:03:45
首先检查一下温度计算结果是否正确，特别是看一下温度梯度结果，是否有过大的地方。在温度梯度过大的部位会造成应力集中，使结构分析不收敛。

温度梯度是正常的啊 第一道加热冷却后再进行的第二道加热 两次的温度场是相似的 第一道温度场对应的应力场可以算 等第二道加热开始时 应力计算就不收敛了 这是怎么回事呢
赞一下(1人) 回复此楼
6楼2016-08-22 09:44:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

匿名

本帖仅楼主可见
2楼2016-08-19 23:51:18
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

sunrui19925

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by whxn at 2016-08-19 23:51:18
正想学习这个

多多交流啊
回复此楼
3楼2016-08-20 20:17:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w_kf557

金虫 (小有名气)
这命令流是自己写的吗？  大神啊
赞一下 回复此楼
4楼2016-08-21 20:09:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
sunrui19925: 金币+5, ★★★很有帮助, 谢谢 希望能私下交流O(∩_∩)O 2016-08-22 10:29:27
首先检查一下温度计算结果是否正确，特别是看一下温度梯度结果，是否有过大的地方。在温度梯度过大的部位会造成应力集中，使结构分析不收敛。
赞一下 回复此楼
5楼2016-08-22 08:03:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunrui19925

铁虫 (小有名气)
引用回帖:
4楼: Originally posted by w_kf557 at 2016-08-21 20:09:50
这命令流是自己写的吗？  大神啊

就是不收敛呢
回复此楼
7楼2016-08-22 09:46:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

曹珍123

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

我也在学习热力耦合，可以咨询一下你吗？

发自小木虫IOS客户端
赞一下 回复此楼
8楼2017-11-24 09:42:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
9楼2017-11-24 09:54:23
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

奥利奥麦旋风

新虫 (初入文坛)
引用回帖:
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
赞一下 回复此楼
10楼2018-02-06 18:07:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sunrui19925 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭