声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1462|回复: 1

[综合讨论] 求高手帮忙,老不收敛

[复制链接]
发表于 2012-11-22 13:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
这是我建的模型 老是不收敛 求高手给看看
/FILNAME,shitanxi001,1
/TITLE, moxing
/REPLOT
/config,fsplit,750
/units,si !定义公制,米、千克、秒
/PREP7  !进入预处理

!******定义桁架参数********
hjps=12       !桁架品数
zgwj=0.245   !主管外径
zgbh=0.007   !主管壁厚
plgwj=0.180  !平联杆外径
plgbh=0.005  !平联杆壁厚
fgwj=0.127    !腹杆外径
fgbh=0.004   !腹杆壁厚
wbgbbh=0.012 !外包钢板壁厚
wbgbw=0.776  !外包钢板宽度
wbgbh=1.320 !外包钢板高度
lbt=0.504       !螺栓杆长度
dbt=0.027      !螺栓杆直径或有效直径
dbth=0.054     !螺栓头和螺母直径
lbth=0.025     !螺栓头厚度
preten=230000   !螺栓施工预拉力 (查要施加多大的预应力)
dh0=dbt+0.001     !螺栓孔直径  
d=0.016  !植筋直径
jjbbh=0.010  !加劲板壁厚

!建立外包钢板wbgb11
wpcsys,wn,kcn
wprota,,,90
blc4,-0.188,-0.22,wbgbw,wbgbh,-wbgbbh
!!在外包钢板上切除圆柱面
cyl4,,,zgwj/2,360,,,-0.1
vgen,2,2,,,,,-0.4
vgen,2,2,3,1,,0.88
vptn,all
vdele,6,8,2,1
vdele,10,12,2,1
cm,wbgb1,volu
!!外包钢板打螺栓孔
wpoff,-0.08,-0.1768,
cylind,dh0/2,,-wbgbbh-0.01,wbgbbh+0.01,0,360
vgen,2,1,,,,,-0.16
vgen,2,1,2,1,,,-0.4
vgen,2,1,4,1,,1.2336,
vsel,s,,,1,6,1,
vsel,a,,,8,10,2
cm,bolt,volu
cmsel,all
vsbv,wbgb1,bolt,keep
!!在外包钢板上切出焊接缝
wpcsys,wn,kcn
vsbw,all
wprota,,90
vsbw,all
wpoff,,,-0.223
vsbw,all
wpoff,,,-0.434
vsbw,all
wpcsys,wn,kcn
wpoff,,0.88
wprota,,90
vsbw,all
wpoff,,-0.4,
wprota,,90
vsbw,all
!!切植筋缝
wpcsys,wn,kcn
vsel,s,,,9
vsel,a,,,31,32,1
wpoff,,,0.08
vsbw,all
wpoff,,,-0.16
vsbw,all
wpoff,,,-0.24
vsbw,all
wpoff,,,-0.16
vsbw,all
vsel,all
cm,wbgb11,volu
!镜像边跨对称面的外包钢板wbgb12
allsel,
wpoff,-0.252,,
csys,4
vsymm,x,all,,,,0,0
cmsel,u,wbgb11
cm,wbgb12,volu

!建立边跨螺栓ls11
csys,0
wpcsys,kn,wcn
wpoff,0,-0.1768,0.08
wprota,,,90
wpoff,,,-lbt
cyl4,,,dbt/2,,,,lbt+4*lbth
vgen,2,71,,,-lbth*2,,,,,1
rpr4,6,,,dbth/2,,-lbth
wpoff,,,lbt
rpr4,6,,,dbth/2,,lbth
allsel
vsel,s,,,71,73,1
vptn,all
vgen,2,all, , , , ,-0.16 ,0
vgen,2,all, , , , ,-0.4 ,0
vgen,2,all, , , ,1.2336 , ,0
cm,ls11,volu

!建立边跨混凝土块hnt11
wpcsys,wn,kcn
wpoff,0,-0.1768,0.08
wprota,,,90
wpoff,-0.8+0.28,0.015,-wbgbbh
wprota,,,90
blc4,,,0.48,1.2036,1.6
cmsel,u,ls11
!!在混凝土htn11上切出植筋的位置
wpcsys,wn,kcn
wpoff,,,0.188
vsbw,all
wpoff,,,-0.188+0.08
vsbw,all
wpoff,,,-0.16
vsbw,all
wpoff,,,-0.24
vsbw,all
wpoff,,,-0.16
vsbw,all
wpoff,,,-0.108
vsbw,all
wprota,,90,
wpoff,,,-0.223
vsbw,all
wpoff,,,-0.434
vsbw,all
cm,hnt11,volu

!mesh外包钢板
allsel,
cmsel,s,wbgb11
cmsel,a,wbgb12

aslv,s
lsla,s
cm,wbgb0line,line

lsel,r,loc,x,0
cm,wbgb11line,line
cmsel,s,wbgb0line
lsel,r,loc,x,-0.504
cm,wbgb12line,line
cmsel,s,wbgb0line
lsel,r,loc,x,-0.504+3.98-0.001,-0.504+3.98+0.001
cm,wbgb21line,line
cmsel,a,wbgb11line
cmsel,a,wbgb12line
cm,wbgbline,line

lsel,r,radius,,dh0/2
lesize,all,,,10
cmsel,s,wbgbline
lsel,r,radius,,zgwj/2
lesize,all,,,15
cmsel,s,wbgbline
lsel,r,length,,0.188-zgwj/2
lesize,all,,,4

cmsel,s,wbgbline
lsel,r,length,,0.22-zgwj/2
lesize,all,,,6
cmsel,s,wbgbline
lsel,r,length,,wbgbh-0.88-0.22-zgwj/2
lesize,all,,,6
cmsel,s,wbgbline
lsel,r,length,,0.4-zgwj
lesize,all,,,10
cmsel,s,wbgbline
lsel,r,length,,0.223-zgwj/2
lesize,all,,,6
cmsel,s,wbgbline
lsel,r,length,,0.434
lesize,all,,,10
cmsel,s,wbgbline
lsel,r,length,,0.88-0.223-0.434-zgwj/2
lesize,all,,,6

cmsel,s,wbgbline
lsel,r,length,,0.4
lesize,all,,,20
cmsel,s,wbgbline
lsel,r,length,,0.22
lesize,all,,,15
cmsel,s,wbgbline
lsel,r,length,,0.223
lesize,all,,,10
cmsel,s,wbgbline
lsel,r,length,,0.88-0.223-0.434
lesize,all,,,10
cmsel,s,wbgbline
lsel,r,length,,wbgbh-0.88-0.22
lesize,all,,,15
cmsel,s,wbgbline
lsel,r,length,,0.188
lesize,all,,,15
cmsel,s,wbgbline
lsel,r,length,,0.24
lesize,all,,,12
cmsel,s,wbgbline
lsel,r,length,,zgwj/2
lesize,all,,,10
cmsel,s,wbgbline
lsel,r,length,,0.08
lesize,all,,,4
cmsel,s,wbgbline
lsel,r,length,,0.108
lesize,all,,,5
cmsel,s,wbgb0line
lsel,r,length,,wbgbbh
lesize,all,,,2
!定义外包钢板单元和材料Q345钢材
et,5,solid45
mp,ex,5,2.06e11
mp,prxy,5,0.3
tb,biso,5,1,2
tbtemp,0
tbdata,,345000000,0,,,,
r,5,
!给外包钢板赋予材料和属性,并划分网格
type,5
mat,5
real,5
vsweep,wbgb11
vsweep,wbgb12
cm,gbe,elem
!定义材料和单元,用于高强螺栓!!!!
et,6,solid45
mp,ex,6,2.06e11
mp,prxy,6,0.3
r,6,
tb,biso,6,1,2
tbtemp,0
tbdata,,940000000,0,,,,

!给高强螺栓赋予属性
allsel,
cmsel,s,ls11
cm,bolt,volu
vatt,6,6,6,0
!!!mesh螺栓杆和螺母!!!
aslv,s,1
lsla,s,1
cm,lsline,line
lsel,r,loc,x,lbth*2
lesize,all,,,3

type,6
mat,6
real,6
!smrtsize,0.1
vsweep,all
cmsel,u,gbe
cm,lse,elem
!!!定义材料,用于混凝土块
et,7,solid65
KEYOPT,7,1,0
KEYOPT,7,5,0
KEYOPT,7,6,0
KEYOPT,7,7,1
mp,ex,7,3e10
mp,prxy,7,0.2
mp,dens,7,2500
mp,ex,8,2e11
mp,prxy,8,0.25
!屈服准则
TB,BISO,8,1,2,
TBTEMP,0
TBDATA,,310000000,20000,,,,
r,7,8,0.01,90,0,8,0.01,0,0    !设置配筋率

tb,conc,7,1,9
tbdata,,0.35,1,3.1125e6,-1  !张开传递系数,裂缝闭合传递系数,单轴抗拉强度,单轴抗压
!屈服准则
TB,MISO,7,1,5,
TBTEMP,0
TBPT,,0.0005,15e6
TBPT,,0.001,21e6
TBPT,,0.0015,24e6
TBPT,,0.002,27e6
TBPT,,0.003,24e6
!!赋予混凝土块材料和属性
allsel,
cmsel,s,hnt11
vatt,7,7,7
!!网格混凝土块
aslv,s
lsla,s
cm,hntline,line
lsel,r,length,,0.48
lesize,all,,,12,,,,,1
cmsel,s,hntline
lsel,r,length,,0.412
lesize,all,,,10,,,,,1
cmsel,s,hntline
lsel,r,length,,0.434
lesize,all,,,10,,,,,1
cmsel,s,hntline
lsel,r,length,,0.3848
lesize,all,,,10,,,,,1
cmsel,s,hntline
lsel,r,length,,0.108
lesize,all,,,5,,,,,1

cmsel,s,hntline
lsel,r,length,,0.16
lesize,all,,,8,,,,,1

!!!mesh混凝土块
type,7
mat,7
real,7
mshkey,1
vsweep,all
cmsel,u,gbe
cmsel,u,lse
cm,hnte,elem

!!!!对高强螺栓psmesh,建立预紧单元!!!!!!!!
cmsel,s,bolt
cm,boltshank1,volu
wpoff,-0.252,,
PSMESH,10,bolt10, ,v,75,4,x,0,,,,bolt_1
PSMESH,11,bolt11, ,V,72,4,x,0,,,,bolt_2
PSMESH,12,bolt12, ,V,89,4,x,0,,,,bolt_3
PSMESH,13,bolt13, ,V,86,4,x,0,,,,bolt_4
PSMESH,14,bolt14, ,v,103,4,x,0,,,,bolt_5
PSMESH,15,bolt15, ,V,100,4,x,0,,,,bolt_6
PSMESH,16,bolt16, ,V,114,4,x,0,,,,bolt_7
PSMESH,17,bolt17, ,V,117,4,x,0,,,,bolt_8
!!1
allsel
cmsel,s,wbgb11
aslv,s
asel,r,loc,x,-wbgbbh
cm,_contact1,area
cmsel,s,hnt11
aslv,s
asel,r,loc,x,-wbgbbh
cm,_target1,area
et,10,targe170
et,11,conta174
keyopt,11,9,0
KEYOPT,11,5,1
!KEYOPT,9,10,1
!KEYOPT,9,11,1
!KEYOPT,9,12,2
mp,mu,9,0.2
r,9,,,0.1,0.5,,    !计算中提示要么增大fkn或者ftoln,ftoln从0.1变成0.5,fkn也改成0.5
real,9
mat,9
cmsel,s,_target1
type,10
nsla,s,1
esln,s,0
esurf,all
allsel
cmsel,s,_contact1
type,11
nsla,s,1
esln,s,0
esurf,all
cmsel,s,wbgb11
aslv,s
asel,r,loc,x,0
cm,_target2,area
allsel
cmsel,s,bolt
vsel,r,,,80
vsel,a,,,84
vsel,a,,,94
vsel,a,,,98
vsel,a,,,108
vsel,a,,,112
vsel,a,,,122
vsel,a,,,126
aslv,s
asel,r,loc,x,0
cm,_contact2,area
et,12,targe170
et,13,conta174
keyopt,13,9,0
KEYOPT,13,5,1
mp,mu,10,0.35
r,10,,,0.1,0.5,,
real,10
mat,10
cmsel,s,_target2
type,12
nsla,s,1
esln,s,0
esurf,all
allsel
cmsel,s,_contact2
type,13
nsla,s,1
esln,s,0
esurf,all

!!2
cmsel,s,wbgb12
aslv,s
asel,r,loc,x,-lbt+wbgbbh
cm,_contact3,area
cmsel,s,hnt11
aslv,s
asel,r,loc,x,-lbt+wbgbbh
cm,_target3,area
et,14,targe170
et,15,conta174
keyopt,15,9,0
KEYOPT,15,5,1
mp,mu,11,0.2
r,11,,,0.1,0.5,,
real,11
mat,11
cmsel,s,_target3
type,14
nsla,s,1
esln,s,0
esurf,all
allsel
cmsel,s,_contact3
type,15
nsla,s,1
esln,s,0
esurf,all
cmsel,s,wbgb12
aslv,s
asel,r,loc,x,-lbt
cm,_target4,area
allsel
cmsel,s,bolt
vsel,r,,,79
vsel,a,,,83
vsel,a,,,93
vsel,a,,,97
vsel,a,,,107
vsel,a,,,111
vsel,a,,,121
vsel,a,,,125
aslv,s
asel,r,loc,x,-lbt
cm,_contact4,area
et,16,targe170
et,17,conta174
keyopt,17,9,0
KEYOPT,17,5,1
mp,mu,12,0.35
r,12,,,0.1,0.5,,
real,12
mat,12
cmsel,s,_target4
type,16
nsla,s,1
esln,s,0
esurf,all
allsel
cmsel,s,_contact4
type,17
nsla,s,1
esln,s,0
esurf,all

!下面施加loadstep1,边界约束条件、螺栓预应力
csys,0
allsel,
!esel,s,type,,10,21,1
!esurf,all,rev        !!反转接触单元和目标单元的法向
!allsel,
finish
/solu
lsclear,all
!约束边界
allsel,

asel,s,loc,z,-1
asel,u,loc,x,1,5
da,all,all
asel,s,loc,z,0.6
asel,u,loc,x,1,5
da,all,all
allsel,
sload,all,9,lock,forc,preten,1,2  !对螺栓施加预应力
nlgeom,on
NROP,UNSYM  !指定牛顿拉夫逊法求解的选项
pred,on
lnsrch,on
antype,0
autots,on
nsubst,100,1000,1
NEQIT,100       !平衡迭代次数  
outres,all,1   !1表示不写入所有子步内容
eqslv,sparse,1e-8
cnvtol,f, ,0.05,2,0.1
allsel
lswrite,1,
solve
save
finish
回复
分享到:

使用道具 举报

发表于 2012-12-17 15:28 | 显示全部楼层
请问楼主解决这个不收敛问题了吗。我现在也遇到这种问题了。
另外我查过资料,说可以增加迭代次数或者放松收敛准则。但是我不会怎么设置。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 19:54 , Processed in 0.071055 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表