钢管混凝土网格划分问题

leozhao

我做的是钢管混凝土柱，下面是我的模型和网格划分的命令流，高手帮我看看我这样划分网格行不行？谢谢啦！
看看我标红的那部分这样划分行不行，谢谢
附我的命令流：finish
/clear
/prep7
et,1,solid65                    !混凝土
et,2,solid45                    !上垫板及薄壁钢管
MP,EX,1,3.35E4          !混凝土
MP,PRXY,1,0.2            
TB,CONCR,1                  
TBDATA,,0.25,0.99,2.64,-1              
TB,KINH,1,1,13
TBTEMP,0
TBPT,,0,0   
TBPT,,0.00001,0.335  
TBPT,,0.001,14.4
TBPT,,0.002,23.082
TBPT,,0.003,26.047
TBPT,,0.004,25.486
TBPT,,0.005,23.006  
TBPT,,0.007,17.306
TBPT,,0.008,15.302  
TBPT,,0.009,13.792   
TBPT,,0.015,9.501   
TBPT,,0.02,8.076
TBPT,,0.025,7.221

mp,ex,2,2.06e5    ! 薄壁钢管
mp,prxy,2,0.3   
TB,BKIN,2,1,2    !BKIN 经典双线性随动强化              
TBTEMP,0        
TBDATA,,345
mp,ex,3,2.06e9    ! 垫板
mp,prxy,3,0.3   
TB,BKIN,3,1,2    !BKIN 经典双线性随动强化              
TBTEMP,0        
TBDATA,,345


l=200    ! 肢长
b=100    !肢宽
h=1580   !总高度
t=3      !钢管厚度
p=7
!!!建立关键点               
k,1,0,0,0
k,2,b,0,0
k,3,l,0,0
k,4,l,b,0
k,5,b,b,0
k,6,b,l,0
k,7,0,l,0
k,8,t,t,0
k,9,b-t,t,0
k,10,b-t,l-t,0
k,11,t,l-t,0
k,12,b,t,0
k,13,l-t,t,0
k,14,l-t,b-t,0
k,15,b,b-t,0
a,1,2,3,4,5,6,7
a,8,9,10,11
a,12,13,14,15
asel,s,,,1
asel,a,,,2,3,1
asba,1,all
a,8,9,10,11
a,12,13,14,15
wprota,,,90
wpoffs,,,t
asbw,all
wpoffs,,,b-2*t
asbw,all
wpoffs,,,t
asbw,all
wpoffs,,,b-t
asbw,all
wprota,,-90
wpoffs,,,t
asbw,all
wpoffs,,,b-2*t
asbw,all
wpoffs,,,t
asbw,all
wpoffs,,,b-t
asbw,all
wpoffs,,,-l+t
wprota,,90
wpoffs,,,-l+t
wprota,,,-90
vext,all,,,,,h+20
wpoffs,,,h
vsbw,all
wpoffs,,,-h
!!!垫块划分
vsel,s,loc,z,h,h+20
vatt,3,,2
esize,20
vmesh,all
ALLSEL APLOT
!!!混凝土划分
vsel,s,loc,z,0,h
vsel,r,loc,x,t,b-t
vsel,r,loc,y,t,l-t
vatt,1,,1
esize,25
vmesh,all
vsel,s,loc,z,0,h
vsel,r,loc,x,b,l-t
vsel,r,loc,y,t,b-t
vatt,1,,1
esize,25
vmesh,all
ALLSEL APLOT
!!!钢管划分
vsel,s,loc,z,0,h
vsel,r,loc,x,0,t
vsel,r,loc,y,t,l-t
vatt,2,,2
esize,25
vmesh,all
vsel,s,loc,z,0,h
vsel,r,loc,x,b-t,b
vsel,r,loc,y,t,l-t
vatt,2,,2
esize,25
vmesh,all
vsel,s,loc,z,0,h
vsel,r,loc,x,l-t,l
vsel,r,loc,y,t,b-t
vatt,2,,2
esize,25
vmesh,all
vsel,s,loc,z,0,h
vsel,r,loc,x,0,l
vsel,r,loc,y,0,t
vatt,2,,2
esize,25
vmesh,all
vsel,s,loc,z,0,h
vsel,r,loc,x,0,b
vsel,r,loc,y,l-t,l
vatt,2,,2
esize,25
vmesh,all
vsel,s,loc,z,0,h
vsel,r,loc,x,b,l
vsel,r,loc,y,b-t,b
vatt,2,,2
esize,25
vmesh,all
NUMMRG,ALL

zengfanyong

不错,做得挺好

lizhi17

想请问一下，混凝土的本构考虑了约束效应了么

jackinsjtu2012

对呀，我也是想问下，这个是不是做了？

leozhao

jackinsjtu2012 发表于 2013-5-9 20:41
对呀，我也是想问下，这个是不是做了？

做的时间很长了，这个混凝土的本构关系考虑约束效应的，

流酸衲

楼主以为你做的是ABAQUS呢 还是感谢