找回密码
 注册
Simdroid-非首页
查看: 55|回复: 4

[问题讨论] 压强大的时候不收敛,请教原因

[复制链接]
发表于 2013-10-28 20:58:30 | 显示全部楼层 |阅读模式 来自 山东济南
本人最近做磨削力引起的残余应力分析,发现在压强小的时候计算收敛,在压强大的时候不收敛,请教大家原因,以及解决办法。
发表于 2013-10-29 20:31:53 | 显示全部楼层 来自 湖北武汉
Simdroid开发平台
首先自己应该去查看计算过程不收敛的因素是什么,比如警告信息,错误文件等等,再判断
回复 不支持

使用道具 举报

 楼主| 发表于 2013-10-29 20:42:38 | 显示全部楼层 来自 山东济南
senja 发表于 2013-10-29 20:31
首先自己应该去查看计算过程不收敛的因素是什么,比如警告信息,错误文件等等,再判断 ...

错误信息为Solution not converged at time 3.333333333E-07 (load step 1 substep 1).  Run terminated.

请教版主什么原因,个人认为是压强过大,导致的网格变形过大,不知如何解决。
附上命令流
/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


回复 不支持

使用道具 举报

发表于 2013-10-29 21:13:16 | 显示全部楼层 来自 湖北武汉
ifly7 发表于 2013-10-29 20:42
错误信息为Solution not converged at time 3.333333333E-07 (load step 1 substep 1).  Run terminated. ...

如果是单元变形过大,那么就找到相应的位置,网格需要细化
回复 不支持

使用道具 举报

发表于 2013-12-23 11:03:38 | 显示全部楼层 来自 江苏淮安
不收敛往往和网格有关系
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 20:40 , Processed in 0.032140 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表