| 查看: 1167 | 回复: 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 |
回复此楼» 猜你喜欢
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有9人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有6人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有4人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有3人回复
-
退学或坚持读
已经有28人回复
-
面上项目申报
已经有3人回复
-
酰胺脱乙酰基
已经有9人回复
-
博士延得我，科研能力直往上蹿
已经有7人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
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 |
哀莫子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