24小时热门版块排行榜    

查看: 1441  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 5/250 2026-02-28 20:11 by iwuli
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[考研] 0856材料求调剂 +10 hyf hyf hyf 2026-02-28 11/550 2026-02-28 18:50 by 无际的草原
[考研] 311求调剂 +7 南迦720 2026-02-28 7/350 2026-02-28 18:28 by leonnulll
[考研] 材料学调剂 +4 提神豆沙包 2026-02-28 4/200 2026-02-28 18:26 by houyaoxu
[教师之家] 版面费该交吗 +15 苹果在哪里 2026-02-22 18/900 2026-02-28 18:20 by mibaomingg
[考研] 材料调剂 +3 爱擦汗的可乐冰 2026-02-28 3/150 2026-02-28 18:06 by houyaoxu
[考研] 化工专硕348,一志愿985求调剂 +3 弗格个 2026-02-28 5/250 2026-02-28 17:04 by sandychj
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考博] 26申博 +3 想申博! 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 295求调剂 +4 19171856320 2026-02-28 4/200 2026-02-28 13:39 by ms629
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[基金申请] 什么是人一生最重要的? +10 瞬息宇宙 2026-02-21 10/500 2026-02-27 08:46 by tfang
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见