请问在端板连接的梁柱分析中加了梁段荷载,计算的两个荷载步能收敛,怎么不见梁端变化呢?第二个荷载步在梁端加集中力,计算收敛后内力变化也很小,是单元shell181没设置好还是接触或是别的问题,接触应该是起作用了,只定义了梁端和柱的接触,请大家看看。弄了好多天了就是不知道怎么回事,对shell181的keyopt反复改0,2了都不行,命令流如下:
*SET,hc,200
*SET,bc,200
*SET,tc,9
*SET,lc,1500
*SET,hb,194
*SET,bb,150
*SET,tbf,9
*SET,tbw,6
*SET,lb,1050
*SET,tep,9 !端板厚度
*SET,bep,200 !端板宽度
*SET,hep,400 !端板高度
*SET,dbt,22
*SET,dbth,33.4
*SET,lbth,12.5
*SET,preten,155000
*SET,miu,0.45
*SET,dho,dbt+2
*SET,aa,70 !螺栓到端板顶的距离
*SET,ab,40 !螺栓到端板边的距离
/prep7
/pnum,area,1
!梁建模
rect,0,lb,-hb/2+tbf/2,hb/2-tbf/2
wprota,,90
wpoffs,,,hb/2-tbf/2
rect,0,lb,0,bb/2
wpoffs,,,-hb+tbf
rect,0,lb,0,bb/2
aglue,all
cm,beam,area
!端板建模
wprota,,,-90
wpave,0,0,0
rect,-hep/2,hep/2,0,bep/2
asel,s,loc,x,0
wprota,,,90
wpoffs,,,hb/2-tbf/2
a***w,all
wpoffs,,,-hb+tbf
a***w,all
allsel,all
aglue,all
!柱建模
wprota,,,-90
wpave,0,0,0
rect,-hep/2-20,hep/2+20,0,bc/2
rect,-lc/2,-hep/2-20,0,bc/2
rect,hep/2+20,lc/2,0,bc/2
wpoffs,,,hc-tc
rect,-lc/2,lc/2,0,bc/2
wpoffs,,bc/2
wprota,,-90
rect,-lc/2,lc/2,0,hc-tc
allsel,all
!建端板螺栓
wpave,0,0,0
wprota,,90
cyl4,hep/2-aa,bep/2-ab,0,,dho/2
asel,s,loc,y,-hep/2+aa
agen,2,all,,,,hep-2*aa
cm,bolto1,area
allsel,all
lsel,s,length,,2*(hep/2+20)
asll,s,0 !包含target的面
asel,u,,,bolto1
cm,target,area
cmsel,a,bolto1
a***a,target,bolto1
cm,target,area
allsel,all
cyl4,hep/2-aa,bep/2-ab,0,,dho/2
asel,s,loc,y,-hep/2+aa
agen,2,all,,,,hep-2*aa
cm,bolto2,area
allsel,all
asel,s,loc,x,0
asel,r,loc,y,-hep/2-20,hep/2+20
asel,u,,,bolto2
asel,u,,,target
cm,plate,area
allsel,all
cmsel,s,plate
cmsel,a,bolto2
a***a,plate,bolto2
cm,endplate,area
!建螺栓
allsel,all
wpoffs,hep/2-aa,bep/2-ab
cylind,dbt/2,,-lbth,0,0,360
cylind,dbth/2,dbt/2,-lbth,0,0,360
vglue,all
cylind,dbt/2,,0,lbth,0,360
cylind,dbth/2,dbt/2,0,lbth,0,360
vsel,s,loc,y,-hep/2+aa
vsel,r,loc,x,-lbth,0
vglue,all
allsel,all
vsel,s,loc,y,-hep/2+aa
vgen,2,all,,,,hep-2*aa
allsel,all
cm,bolt,volu
aslv,s
cm,bolta,area
save
!cm梁柱螺栓
asel,s,loc,x,-hc+tc,0
asel,u,,,endplate
asel,u,,,bolta
aglue,all
cm,column,area
lsel,s,length,,2*(hep/2+20)
asll,s,0 !包含target的面
asel,r,loc,x,0
cm,target,area
allsel,all
asel,s,loc,x,lbth+10,lb
cm,beam,area
allsel,all
!单元属性
et,1,solid45
mp,ex,1,206e3
mp,prxy,1,0.3
tb,biso,1,1
tbtemp,0
tbdata,,940
et,2,shell181
mp,ex,2,206e3
mp,prxy,2,0.3
tb,biso,2,1
tbtemp,0
tbdata,,235
et,3,mesh200
!体积面积性质
allsel,all
vatt,1,,1,0
allsel,all
asel,s,,,beam
asel,r,loc,y,0
sect,1,shell
secdata,tbw/2,2,0,3
r,3,tbw/2
aatt,2,3,2,0,1
allsel,all
asel,s,,,beam
asel,u,loc,y,0
r,4,tbf
sect,2,shell
secdata,tbf,2,0,3
aatt,2,4,2,0,2
allsel,all
asel,s,,,column
asel,r,loc,x,0
r,5,tc
sect,3,shell
secdata,tc,2,0,3
secoffset,bot
aatt,2,5,2,0,3
allsel,all
asel,s,,,column
asel,u,loc,x,0
r,6,tc
sect,4,shell
secdata,tc,2,0,3
aatt,2,6,2,0,4
allsel,all
asel,s,,,endplate
r,7,tep
sect,5,shell
secdata,tep,2,0,3
aatt,2,6,2,0,5
!mesh螺栓梁柱
allsel,all
lsel,s,radius,,dbt/2
asll,s,1
vsla,s,0
cm,boltlink,volu
vsel,r,loc,x,0,lbth
aslv,s
asel,r,loc,x,lbth
lsel,s,radius,,dbt/2
lesize,all,,,4
lsel,s,length,,lbth
lesize,all,,,3
type,3
keyopt,3,1,4
mshkey,0
mshape,1,2d
amesh,all
extopt,aclear,1
vsweep,all
allsel,all
vsel,s,,,boltlink
vsel,r,loc,x,-lbth,0
aslv,s
asel,r,loc,x,0
lsel,s,radius,,dbt/2
lesize,all,,,4
lsel,s,length,,lbth
lesize,all,,,3
type,3
keyopt,3,1,4
mshkey,0
mshape,1,2d
amesh,all
extopt,aclear,1
vsweep,all
allsel,all
vsel,s,,,bolt
vsel,u,,,boltlink
vsel,r,loc,x,0,lbth
lsel,s,length,,lbth
lesize,all,,,3
lsel,s,radius,,dbth/2
lesize,all,,,4
mshkey,1
mshape,0,3d
vsweep,all
allsel,all
vsel,s,,,bolt
vsel,u,,,boltlink
vsel,r,loc,x,-lbth,0
lsel,s,length,,lbth
lesize,all,,,3
lsel,s,radius,,dbth/2
lesize,all,,,4
mshkey,1
mshape,0,3d
vsweep,all
!mesh柱
allsel,all
asel,s,,,column
lsel,s,radius,,dho/2
lesize,all,,,4
lsel,s,length,,bc/2
lesize,all,,,10
lsel,s,length,,hc-tc
lesize,all,,,15
lsel,s,length,,lc
lesize,all,,,120
lsel,s,length,,2*(hep/2+20)
lesize,all,,,36
lsel,s,length,,lc/2-hep/2-20
lesize,all,,,42
mshkey,0
mshape,0,2d
amesh,all
!mesh梁
allsel,all
asel,s,,,beam
lsel,s,length,,lb
lesize,all,,,60
lsel,s,length,,bb/2
lesize,all,,,8
lsel,s,length,,hb-tbf
lesize,all,,,20
amesh,all
!mesh端板
allsel,all
asel,s,,,endplate
lsel,s,length,,bep/2
lesize,all,,,10
lsel,s,length,,(hep-hb+tbf)/2
lesize,all,,,8
amesh,all
!ce自由度
allsel,all
vsel,s,,,bolt
vsel,s,loc,x,0,lbth
vsel,u,,,boltlink
aslv,s
asel,r,loc,x,0
cm,boltc1,area
allsel,all
vsel,s,,,bolt
vsel,s,loc,x,-lbth,0
vsel,u,,,boltlink
aslv,s
asel,r,loc,x,0
cm,boltc2,area
allsel,all
asel,s,,,endplate
esla,s
asel,s,,,boltc1
nsla,s,1
ceintf,1,ux,uy,uz
allsel,all
asel,s,,,target
esla,s
asel,s,,,boltc2
nsla,s,1
ceintf,1,ux,uy,uz
!接触建立
allsel,all
cmsel,s,target
cmsel,a,endplate
r,7,,,1,0.2
real,7
et,4,170
et,5,174 !contact
keyopt,5,9,0
keyopt,5,12,0
asel,s,,,target
type,4
nsla,s,1
esln,s,1
/psymb,esys,1
esurf,,bottom
mp,mu,3,miu
mat,3
asel,s,,,endplate
real,7
type,5
nsla,s,1
esln,s,1
esurf,,top
allsel,all
!对高强螺栓psmesh
vsel,s,,,boltlink
vsel,r,loc,y,hep/2-aa
aslv,s
asel,r,loc,x,0
nsla,s,0
nummrg,node
esln,s,0
r,8
real,8
psmesh,6,,,all,,0,x,0
allsel,all
vsel,s,,,boltlink
vsel,r,loc,y,-hep/2+aa
aslv,s
asel,r,loc,x,0
nsla,s,0
nummrg,node
esln,s,0
real,8
psmesh,7,,,all,,0,x,0
!下面施加荷载和位移条件求解
allsel,all
/solu
antype,0
nlgeom,1
nsubst,10
outres,basic,1
neqit,200
cnvtol,f,,0.01,2
lsel,s,loc,z,0
dl,all,,symm
allsel,all
nsel,s,loc,y,-lc/2
nsel,a,loc,y,lc/2
d,all,ux,0,,,,uy,uz
allsel,all
sload,all,9,lock,forc,preten,1,2
time,1
solve
!下面施加 loadstep2
allsel,all
nlgeom,1
nsubst,20,5000,1
outres,basic,1
cnvtol,f,,0.01,2
lsel,s,loc,z,0
dl,all,,symm
allsel,all
nsel,s,loc,y,-lc/2
nsel,a,loc,y,lc/2
d,all,ux,0,,,,uy,uz
allsel,all
num=node(lb-50,hb/2-tbf/2,0)
d,num,uy,-200
allsel,all
time,2
solve
save
模型的图贴上,请大家看看怎么回事呢? |