| 查看: 2850 | 回复: 9 | ||
[求助]
压强大的时候不收敛，请教原因
|
| 本人最近做磨削力引起的残余应力分析，发现在压强小的时候计算收敛，在压强大的时候不收敛，请教大家原因，以及解决办法。 |
回复此楼» 猜你喜欢
085602调剂 初试总分335
已经有5人回复
-
302分 085601求调剂推荐
已经有8人回复
-
生物学308分求调剂（一志愿华东师大）
已经有10人回复
-
材料专硕(0856) 339分求调剂
已经有8人回复
-
085410人工智能 初试316分 求调剂
已经有6人回复
-
282求调剂
已经有10人回复
-
070300化学求调剂
已经有17人回复
-
318求调剂
已经有12人回复
-
308求调剂
已经有3人回复
-
315求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯计算求助！！！
已经有5人回复
- 工艺问题与收敛问题求教
已经有6人回复
- 请教，若迭代矩阵的谱半径恰好等于1，那么迭代是收敛还是发散的？
已经有8人回复
- 【已完结】abaqus分析不收敛
已经有6人回复
- ms结果出问题了，又出现错误，castep里面有个warning，如下
已经有7人回复
- 怎么计算不同温度下，100%甲醇的密度？
已经有4人回复
- 请教各位一个级数收敛问题
已经有7人回复
- 接触分析不收敛
已经有3人回复
- 为什么连续性方程一直不收敛
已经有6人回复
- 请教如何加快siesta的收敛速度
已经有9人回复
- 求助有关高斯优化和频率计算的收敛问题
已经有4人回复
- 计算graphene的问题
已经有10人回复
- 请教关于收敛精度
已经有3人回复
- Fluent连续性方程不收敛
已经有13人回复
- 想请问一下，几何优化之后的晶格角度变化的问题
已经有12人回复
- 请教一下拓扑和收敛的关系
已经有9人回复
- MS5.5 Castep 优化的时候 +U
已经有6人回复
- fluent中关于拉瓦尔喷管的问题
已经有26人回复
- 声子计算又负频的出现！！！
已经有24人回复
- 【求助】请教如何在MS，castep模块里面实现加压强这一功能？
已经有19人回复
- 【求助】计算单个Fe原子，无法收敛，请教
已经有19人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
2楼2013-10-29 07:32:53
3楼2013-10-29 08:44:08
4楼2013-10-29 09:39:19
5楼2013-10-29 17:13:54
6楼2013-10-29 18:00:23
7楼2013-10-30 07:08:58
|
命令流附上，没研究过屈曲，个人感觉，原因是压强过大，造成单元变形过大，不知如何解决？网格加密没效果，多子步，没效果。 /filename,force /prep7 /units,si /view,1,,,1 et,1,plane182 keyopt,1,1,0 keyopt,1,3,2 keyopt,1,6,0 mptemp,,,,,,,, mptemp,1,0 mpdata,ex,1,,2.19e+011 mpdata,prxy,1,,0.3 mptemp,,,,,,,, mptemp,1,0 mpdata,dens,1,,7830 sigma=1617e+6 tb,bkin,1,1,2,1 tbtemp,0 tbdata,, sigma,0.03*2.19e+11,,,, lc=3/1000 !磨削弧长是3mm stepsize=lc inlc=54 !lc分成inlc份 n=20 !工件长度是n*lc esize=stepsize/ inlc k,1, 0,10/1000,0, k,2, 0,0,0, k,3, n*lc,10/1000,0, k,4, n*lc,0,0, a,1,2,4,3 lesize, 1, , , 30,10 , , , ,1 lesize, 3 , , ,30,0.1 , , , ,1 lesize, 4 , , , inlc*n, , , , ,1 mshape,0,2d mshkey,1 amesh, 1 vw=1000/1000/60 !磨削速度 tinc=esize/vw tm=n*stepsize/vw xi=0 !热源起点 /solu antype,trans trnopt,full lsel,s, , , 2 nsll,s,1 d,all,all i=1 *do,t, tinc,tm,tinc NLGEOM,1 NSUBST,50,10000,20 AUTOTS,1 TIME,t CNVTOL,F, ,0.005,2, , outres,all,all tsres,erase xi=(i-1)*esize allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi-1e-6, xi+lc/9+1e-6 sf,all,pres,-1e+10/3 allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi+2*lc/9-1e-6, xi+3*lc/9+1e-6 sf,all,pres,-1e+10/3 allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi+4*lc/9-1e-6, xi+5*lc/9+1e-6 sf,all,pres,-1e+10/3 allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi+6*lc/9-1e-6, xi+7*lc/9+1e-6 sf,all,pres,-1e+10/3 allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi+8*lc/9-1e-6, xi+9*lc/9+1e-6 sf,all,pres,-1e+10/3 ! sfgrad,pres, ,y,xi,-60* sigma, !sfgrad,pres, ,x,xi,17/(3/1000*10/1000), !sf,all, pres,1e-6 allsel solve allsel lsel,s,,,4 nsll,s,1 sfdele,all, pres i=i+1 *enddo |
8楼2013-10-30 08:57:27
【答案】应助回帖
|
试了一下你的命令流，似乎是求解设置有问题。对命令流做了一些修改，然后可以计算了 (用你的命令流，一执行 Solve 命令就马上报错并停止计算)，不过需要的时间步长太小，估计需要算很长时间，所以算了一段时间就强制退出了。下面是修改后的命令流，修改的地方用 !*** 做了标记，你试试看，然后自己再修改一下吧： fini /clear /filename,force !* /prep7 /units,si /view,1,,,1 et,1,plane182 keyopt,1,1,0 keyopt,1,3,2 keyopt,1,6,0 !* mptemp,,,,,,,, mptemp,1,0 mpdata,ex,1,,2.19e+011 mpdata,prxy,1,,0.3 mptemp,,,,,,,, mptemp,1,0 mpdata,dens,1,,7830 sigma=1617e+6 tb,bkin,1,1,2,1 tbtemp,0 tbdata,, sigma,0.03*2.19e+11,,,, lc=3/1000 ! 磨削弧长是3mm = 0.003 stepsize=lc ! stepsize = 0.003 inlc=54 ! lc分成inlc份 n=20 ! 工件长度是 n*lc = 0.06 esize=stepsize/ inlc ! esize = 0.003 / 54 = 0.00005556 k,1, 0,10/1000,0, k,2, 0,0,0, k,3, n*lc,10/1000,0, k,4, n*lc,0,0, a,1,2,4,3 lesize, 1, , , 30,10 , , , ,1 lesize, 3 , , ,30,0.1 , , , ,1 lesize, 4 , , , inlc*n, , , , ,1 ! Line 4 - 单元分段数 54*20 mshape,0,2d mshkey,1 amesh, 1 vw=1000/1000/60 ! 磨削速度 = 1/60 = 0.01667 tinc=esize/vw ! tinc = 0.003333 tm=n*stepsize/vw ! tm = 3.6 xi=0 ! 热源起点 /solu antype,trans trnopt,full lsel,s, , , 2 nsll,s,1 d,all,all iin = nint(tm/tinc) ! 总计算步数 i=1 ! *do, t, tinc, tm,tinc *do, ii, 1, 3, 1 ! iin, 1 !*** 采用整数计数比较好； 先试 3 个循环 t = ii*tinc NLGEOM,1 NSUBST,10,10000,5 !*** 子步数开始可小一点，最大值可以大一点 AUTOTS,1 TIME,t CNVTOL,F, ,0.005,2, -1, !*** 未设置最小参考自由度值，第 5 个参数应设为负值 kbc,0 !*** 设置为 ramped 方式加载，steped 方式载荷一下子到载荷子步的最大值，不易收敛 /GST, 1 !*** 显示收敛追踪曲线 outres,all,all tsres,erase xi=(i-1)*esize allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi-1e-6, xi+lc/9+1e-6 sf,all,pres,-1e+10/3 allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi+2*lc/9-1e-6, xi+3*lc/9+1e-6 sf,all,pres,-1e+10/3 allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi+4*lc/9-1e-6, xi+5*lc/9+1e-6 sf,all,pres,-1e+10/3 allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi+6*lc/9-1e-6, xi+7*lc/9+1e-6 sf,all,pres,-1e+10/3 allsel lsel,s,,,4 nsll,s,1 nsel,r,loc,x, xi+8*lc/9-1e-6, xi+9*lc/9+1e-6 sf,all,pres,-1e+10/3 ! sfgrad,pres, ,y,xi,-60* sigma, ! sfgrad,pres, ,x,xi,17/(3/1000*10/1000), ! sf,all, pres,1e-6 allsel solve allsel lsel,s,,,4 nsll,s,1 sfdele,all, pres i=i+1 *enddo |
» 本帖已获得的红花（最新10朵）
我是东成西就9楼2013-10-31 10:36:19
10楼2013-10-31 14:26:54
