| 查看: 3903 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
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 |
回复此楼» 猜你喜欢
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助中 关于ANSYS接触热力耦合问题
已经有0人回复
- ansys有限元分析教程
已经有3人回复
- 中科院上海光学精密机械研究所招聘
已经有0人回复
|
大神们 可以帮我看看这个程序流吗 /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
 本帖仅楼主可见 |
2楼2016-08-19 23:51:18
sunrui19925
铁虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 38.1
- 红花: 5
- 帖子: 73
- 在线: 164.6小时
- 虫号: 3433382
- 注册: 2014-09-22
- 性别: MM
- 专业: 零件成形制造
3楼2016-08-20 20:17:55
w_kf557
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1937.5
- 散金: 66
- 帖子: 148
- 在线: 41.1小时
- 虫号: 3408975
- 注册: 2014-09-11
- 性别: GG
- 专业: 机械摩擦学与表面技术
4楼2016-08-21 20:09:50
30