24小时热门版块排行榜    

查看: 564  |  回复: 0

王禄友

新虫 (初入文坛)

[求助] ansys分析不收敛

/prep7
!定义线的半径、线的长度、线间距、计算线根数
Bro=0.1
Bl=20
Bs=0.5
Nbeam=bl/bs
!定义单元类型,beam188定义两种单元类型为了便于选择操作,而接触单元的keyopt不同
!定义单元实常数和材料常数、线线摩擦系数、线截面
Et,1,beam188
Et,2,beam188
Et,3,conta176
Et,4,targe170
Keyopt,3,3,1
Keyopt,3,12,2
R,1
R,2,bro,bro,0.001
Mp,ex,1,1.84e5
Mp,prxy,1,0.29
Mp,dens,1,7.85e-6
Mp,mu,2,0.1
Sectype,1,beam,csolid
Secdata,bro
!创建梁模型:此种方法可保证线方向一致,而布尔运算则不一定
!圆截面可不设方位点,但设置方位点可保证变形后的单元显示正确
K,1,,,bs/2
*do,i,2,nbeam+1
K,i,bs/2+(i-2)*bs,,bs/2
*enddo
K,nbeam+2,bl,,bs/2
*do,i,1,nbeam+1
L,i,i+1
*enddo
Lgen,nbeam,all,,,,,bs
Cm,cml1,line
Lsel,none
Kmax=kpinqr(0,14)
K,kmax+1,bs/2
*do,i,2,nbeam+1
K,kmax+i,bs/2,,bs/2+(i-2)*bs
*enddo
K,kmax+nbeam+2,bs/2,,bl
*do,i,1,nbeam+1
L,kmax+i,kmax+i+1
*enddo
Lgen,nbeam,all,,,bs
Cm,cml2,line
Lsel,all
Kmax=kpinqr(0,14)+1
K,kmax,0,1000*b
Lsel,s,tan1,x
Latt,1,1,1,,,kmax,1
Lsel,s,tan1,z
Latt,1,1,2,,,kmax,1
Lsel,all
Esize,bs/8
Lmesh,all
!创建纵线与横线之间的接触对,通过实常数2识别
Type,3
Real,2
Mat,2
Esel,s,type,,1
Nsle
Esurf
Esel,s,type,,2
Nsle
Type,4
Esurf
!施加边界条件
Nsel,s,loc,x,0
Nsel,a,loc,x,bl
Nsel,a,loc,z,0
Nsel,a,loc,z,bl
D,all,all
Cm,dnode1,node
Allsel,all
!在交叉处的重合节点施加反向力对,使其发生上下穿越变形,此力可略大些,可先脱空穿越
*do,i,1,nbeam
*do,j,1,nbeam
*if,mod(i+j,2),eq,0,then
Ftemp=125
*else
Ftemp=-125
*endif
Esel,s,type,,1
Nsle,s
F,node(bs/2+(i-1)*bs,0,bs/2+(j-1)*bs),fy,ftemp
Esel,s,type,,2
Nsle,s
F,node(bs/2+(i-1)*bs,0,bs/2+(j-1)*bs),fy,-ftemp
*enddo
*enddo
!进入求解层,设置求解选项;第一载荷步杀死全部接触对并求解
/solu
antype,static
rescontrol,define,all,1
Nlgeom,on
Outres,all,all
Nsubst,20,,5
Time,1
Esel,s,type,,3,4
Ekill,all
Allsel,all
solve
!第二载荷步定义收敛准则及载荷子步,激活纵横梁之间的接触对,将力对设置为0值,求解
time,2
cnvtol,f,,0.005
nsubst,100,,20
esel,s,type,,3,4
ealive,all
esel,s,type,,1,2
nsle,s
cmsel,u,dnode1
f,all,fy,0
allsel,all
solve
finish
如果把线的长度bl改为10就可以收敛了,但是真实情况是45呀,新人,只有这么多金币,求大神帮忙!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 王禄友 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见