| 查看: 1396 | 回复: 4 | ||||
[求助]
三铰拱非线性屈曲分析的收敛问题
|
|
需要计算一个三铰拱的非线性屈曲荷载，可以正常进行特征值屈曲分析，但之后的非线性部分总是不收敛，希望有人能解答一下。具体命令流如下： finish /clear /filname,gong /prep7 et,1,beam3 mp,ex,1,2e11 mp,prxy,1,0.3 r,1,0.011309734,2.89812E-06,0.2 w=0.1 k,2,0,0$k,3, k,1,-1.663,-0.96$k,4,1.663,-0.96 l,1,2 l,3,4 allsel lesize,all,,,10 lmesh,all cp,next,ux,2,12 cp,next,uy,2,12 allsel dk,1,ux dk,1,uy dk,4,ux dk,4,uy /solu fk,2,fy,-1 antype,0 !!!!!静力分析 pstres,on solve finish /solu antype,1 !!!!屈曲分析 bucopt,lanb,1 mxpand,1 outres,all,all solve finish /prep7 upgeom,0.01,1,1,gong,rst !!!施加初始缺陷 finish !!!!!求解 /solu antype,0 nlgeom,on !!!!!大变形效应 pstres,on !!!!!预应力效应 nsubst,100 !!!子步数 arclen,on !!!激活弧长法 arctrm,u,1.7,,uy outpr,all outres,all,all allsel,all fk,3,fy,-2E6 solve finish !!!!后处理 /post26 nsol,2,2,u,y prod,3,2,,,uy,,,-1000 prod,4,1,,,p,,,2e6 xvar,3 plvar,4 /axlab,x,uy(mm) /axlab,y,q(kN) lines,9999 prvar,3,4 /post1 set,list |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
问题处在两个 cp 命令上，将 cp 命令去掉，改为合并节点就可以收敛了。修改后的命令流如下，修改部分见 !*** 作为标记的地方： finish /clear /filname,gong /prep7 et,1,beam3 mp,ex,1,2e11 mp,prxy,1,0.3 r,1,0.011309734,2.89812E-06,0.2 w=0.1 k,2,0,0$k,3, k,1,-1.663,-0.96$k,4,1.663,-0.96 l,1,2 l,3,4 allsel lesize,all,,,10 lmesh,all !*** cp,next,ux,2,12 !*** cp,next,uy,2,12 NUMMRG,NODE, , , ,LOW !*** allsel dk,1,ux dk,1,uy dk,4,ux dk,4,uy /solu fk,2,fy,-1 antype,0 !!!!!静力分析 pstres,on solve finish /solu antype,1 !!!!屈曲分析 bucopt,lanb,1 mxpand,1 outres,all,all solve finish /prep7 upgeom,0.01,1,1,gong,rst !!!施加初始缺陷 finish !!!!!求解 /solu antype,0 nlgeom,on !!!!!大变形效应 pstres,on !!!!!预应力效应 nsubst,100 !!!子步数 arclen,on !!!激活弧长法 arctrm,u,1.7,,uy outpr,all outres,all,all allsel,all fk,3,fy,-2E6 solve finish !!!!后处理 /post26 nsol,2,2,u,y prod,3,2,,,uy,,,-1000 prod,4,1,,,p,,,2e6 xvar,3 plvar,4 /axlab,x,uy(mm) /axlab,y,q(kN) lines,9999 prvar,3,4 /post1 set,list |
» 本帖已获得的红花（最新10朵）
哀莫子2楼2016-01-01 08:06:40
3楼2016-01-01 13:13:50
zyh20150930
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 925.1
- 散金: 236
- 红花: 3
- 帖子: 259
- 在线: 37.6小时
- 虫号: 4113529
- 注册: 2015-09-30
- 专业: 结构工程
4楼2016-01-17 20:24:33
5楼2016-01-18 23:46:58
