24小时热门版块排行榜    

查看: 1842  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

bluecmy

新虫 (小有名气)

[交流] 【求助】非线性分析得到的结果不收敛,怎么改进 已有3人参与

这里有两个模型,(1)梁,(2)圆柱壳。
用梁得到的结果很好,采用相同的求解设置,圆柱壳的结果不收敛,最后得到的变形也不好不知道什么原因。以下是命令流。
(1)
/prep7           !模型一
et,1,beam3

mp,ex,1,200000
mp,prxy,1,.3

r,1,100,833.333,10

k,1,0,0,0
k,2,0,100,0
l,1,2
esize,1
lmesh,all
finish

/solu
antype,0
nlgeom,1
outres,all,all
nsubst,20
neqit,1000
autots,on
lnsrch,on


dk,1,all,0
fk,2,fy,-50000
fk,2,fx,-250

solve
finish


(2)

t=0.016
L=8
R=0.408
s=0.06
d=0.2
opint=1500000
*set,max_defect,0.0002
/PREP7
ET,1,SHELL93
ET,2,mass21
R,2,1e-6
R,1,t
MP,EX,1,2e11
MP,PRXY,1,0.3

K,1,0,0,0
K,2,0,0,L
LSTR, 1, 2
circle,1,R
lsel,s,loc,z,0
ADRAG,all,0,0,0,,,1
ldele,1
aesize,all,s
MSHAPE,0,2D
MSHKEY,0
amesh,all


et,3,184                    !弯矩
keyopt,3,1,1
nsel,s,loc,z,L,L+1
n,opint ,0,0,L+d
type,3
*get,nnum,node,0,count
*get,ND,node,0,num,min
*do,i,2,nnum
E, opint,ND
ND=NDNEXT(ND)
*enddo
allsel



/solu                    !特征值
antype,0
nsel,s,,,opint
f,all,mx,1

lsel,s,loc,z,0
dl,all,,ux
dl,all,,uy
dl,all,,uz

allsel
pstres,on
solve
finish

/solu
antype,1
bucopt,lanb,1
mxpand,1
solve
finish

/post1
!set,list
set,first
pldisp,0

*get,pressure1,mode,1,freq
nsort,u,sum
*get,max_disp_sum,sort,,max
*set,scale_factor,max_defect/max_disp_sum
finish

/prep7
upgeom,scale_factor,last,last,'shell_mpc','rst',''

/solu
antype,0
nsel,s,,,opint
f,all,mx,pressure1

!nlgeom,1
!nsubst,10,0,0
!arclen,1,4,0.001
!outres,all,all
!arctrm,l
!allsel
!solve
!finish

allsel
nlgeom,1
outres,all,all
nsubst,20
neqit,1000
autots,on
!lnsrch,on
solve
finish
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluecmy

新虫 (小有名气)

引用回帖:
Originally posted by rrainbow at 2010-04-19 11:47:38:
增加子步数,增加迭代次数,改用弧长法

这些方法我都试过了,但结果都不好。是不是具体的参数(比如载荷步加大到多少)很不容易确定
3楼2010-04-29 10:15:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

rrainbow

铜虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
adu886886(金币+1):谢谢提供意见 2010-04-19 14:43
增加子步数,增加迭代次数,改用弧长法
资源有限,交流无限
2楼2010-04-19 11:47:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见