24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

bluecmy

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 189.3
帖子: 69
在线: 14小时
虫号: 840812
注册: 2009-09-05

[交流] 【求助】非线性分析得到的结果不收敛，怎么改进已有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

回复此楼