24小时热门版块排行榜    

查看: 716  |  回复: 0

安仁兵

新虫 (初入文坛)

[求助] ansys初学者自己编写的命令流,出现大变形结构不收敛,求大神帮我找找问题

附命令流,请大神解答
FINISH
/CLEAR                                                                  
/PREP7
!定义单元类型
ET,1,beam188                            !定义钢筋的单元类型
sectype,1,beam,csolid
secdata,6,,
MP,EX,1,2.06E5
MP,PRXY,1,0.25
MP,DENS,1,7850
ET,2,beam188  
sectype,2,beam,csolid
secdata,4,,
MP,EX,2,2.06E5
MP,PRXY,2,0.25
MP,DENS,2,7850
TB,BISO,1,1,2,1                                             !双线性各向同性强化模型
TBDATA,,235,0
TB,BISO,2,1,2,1
TBDATA,,300,0
ET,3,SOLID65                            !定义混凝土的单元类型
MP,EX,3,13585
MP,PRXY,3,0.2
fc=14.3$ft=1.43
tb,concr,3$tbdata,,0.5,0.95,ft,-1
tb,miso,3,,11
tbpt,,0.0002,fc*0.19$tbpt,,0.0004,fc*0.36$tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64$tbpt,,0.0010,fc*0.75$tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91$tbpt,,0.0016,fc*0.96$tbpt,,0.0018,fc*0.99
ET,4,solid45                            !定义钢管的单元类型
MP,EX,4,2.06E5
MP,PRXY,4,0.25
MP,DENS,4,7850
TB,BISO,4,1,2,1                                             !双线性各向同性强化模型
TBDATA,,235,0

!建立几何模型(钢筋)
k,1,20,20,820
k,2,20,20,1280
l,1,2
k,3,20,80,820
k,4,20,80,1280
l,3,4
k,5,80,20,820
k,6,80,20,1020
l,5,6
k,7,80,20,1080
k,8,80,20,1280
l,7,8
k,9,80,80,820
k,10,80,80,1020
l,9,10
k,11,80,80,1080
k,12,80,80,1280
l,11,12
k,13,280,20,1020
l,6,13
k,14,280,80,1020
l,10,14
k,15,280,20,1080
l,7,15
k,16,280,80,1080
l,11,16
k,17,20,20,860
k,18,20,80,860
k,19,80,80,860
k,20,80,20,860
l,17,18
l,18,19
l,19,20
l,20,17
lglue,11,12,13,14
lgen,3,11,14,1,,,60
k,29,20,20,1120
k,30,20,80,1120
k,31,80,80,1120
k,32,80,20,1120
l,29,30
l,30,31
l,31,32
l,32,29
lglue,23,24,25,26
lgen,3,23,26,1,,,60
k,41,120,20,1020
k,42,120,80,1020
k,43,120,80,1080
k,44,120,20,1080
l,41,42
l,42,43
l,43,44
l,44,41
lglue,35,36,37,38
lgen,3,35,38,1,60
!建立几何模型(混凝土)
k,53,,100
k,54,100,100
k,55,100
k,56,,,
a,53,54,55,56
vext,all,,,,,2100
k,100,100,100,1100
k,101,100,100,1000
k,102,100,,1000
k,103,100,,1100
a,100,101,102,103
vext,7,,,1000
vglue,all

!钢筋划分
elemsiz=20
lsel,s,line,,1,10
latt,1,,1$lesize,all,elemsiz
lsel,s,line,,11,46
latt,2,,1$lesize,all,elemsiz
lsel,s,type,,1
lmesh,all

!建立外包钢管
k,201,-2,102
k,202,102,102
k,203,102,-2
k,204,-2,-2
k,205,,100
k,206,100,100
k,207,100,
k,208,,
a,201,202,203,204
a,205,206,207,208
asba,3,4
vext,5,,,,,1000
vgen,2,1,,,,,1100
k,209,102,102,1102
k,210,102,102,998
k,211,102,-2,998
k,212,102,-2,1102
k,213,102,100,1100
k,214,102,100,1000
k,215,102,,1000
k,216,102,,1100
a,209,210,211,212
a,213,214,215,216
asba,34,35
vext,36,,,998

wpcsys,-1$wprota,,-90$vsbw,1$vsbw,4$vsbw,5
wpcsys,-1$wpoff,,100$wprota,,-90$vsbw,7$vsbw,8$vsbw,9
wpcsys,-1$wprota,,,-90$vsbw,11$vsbw,13$vsbw,1$vsbw,6
wpcsys,-1$wpoff,100$wprota,,,-90$vsbw,16$vsbw,17$vsbw,18$vsbw,19
wpcsys,-1$wpoff,,,1000$vsbw,4$vsbw,15
wpcsys,-1$wpoff,,,1100$vsbw,24$vsbw,25
wpcsys,-1$wpoff,1000$wprota,,,-90$asbw,all
wpcsys,-1

esize,20
vsel,s,loc,x,,100$vsel,r,loc,y,,100$vsel,r,loc,z,0,2100
vatt,3,,3$vsweep,all
vsel,all
vsel,s,loc,x,100,1100$vsel,r,loc,y,,100$vsel,r,loc,z,1000,1100
vatt,3,,3$vsweep,all
vsel,all
vsel,s,loc,x,-2,0$vsel,r,loc,y,-2,102$vsel,r,loc,z,0,2100
vatt,4,,4$vsweep,all
vsel,all
vsel,s,loc,x,60,102$vsel,r,loc,y,-2,102$vsel,r,loc,z,0,2100
vatt,4,,4$vsweep,all
vsel,all
vsel,s,loc,x,-2,2100$vsel,r,loc,y,100,102$vsel,r,loc,z,0,2100
vatt,4,,4$vsweep,all
vsel,all
vsel,s,loc,x,-2,2100$vsel,r,loc,y,-2,0$vsel,r,loc,z,0,2100
vatt,4,,4$vsweep,all
vsel,all
vsel,s,loc,x,100,1100$vsel,r,loc,y,-2,102$vsel,r,loc,z,998,1000
vatt,4,,4$vsweep,all
vsel,all
vsel,s,loc,x,100,1100$vsel,r,loc,y,-2,102$vsel,r,loc,z,1100,1102
vatt,4,,4$vsweep,all

nummrg,node


/solu
asel,s,loc,z,0
da,all,all
asel,s,loc,z,2100
da,all,all
allsel,all
sfa,169,1,pres,10
allsel,all
antype,0
nsubst,100
outres,all,all
autos,on
neqit,100
cnvtol,u,,0.05
solve
回复此楼

» 收录本帖的淘帖专辑推荐

ANSYS命令流

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 安仁兵 的主题更新
信息提示
请填处理意见