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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » ansys热力耦合分析不收敛
121/2返回列表12下一页
查看: 3830  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 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 学教育滴
信息提示
关闭
请填处理意见
关闭