24小时热门版块排行榜    

查看: 1419  |  回复: 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的回帖

rrainbow

铜虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
adu886886(金币+1):谢谢提供意见 2010-04-19 14:43
增加子步数,增加迭代次数,改用弧长法
资源有限,交流无限
2楼2010-04-19 11:47:38
已阅   回复此楼   关注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的回帖

rrainbow

铜虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
adu886886(金币+1):谢谢应助 2010-05-04 08:10
nsubst,200,0,20
资源有限,交流无限
4楼2010-05-03 13:37:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
5楼2015-05-15 17:14:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bluecmy 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见