24小时热门版块排行榜    

查看: 1398  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见