试了一下你的命令流,发现如果进一步切割会有很大困难,不如直接对方管部分划分四面体网格。由于方管部分的刚度远大于混凝土,厚度方向划分一层单元也是可以的。我采用单元尺寸 40 (与其它部分相同),可以成功划分四面体网格,不过有很多警告:由 solid185 退化得到的四面体单元 (4 节点四面体) 计算精度很差,不推荐。后来改用 solid186 单元,同样可以划分四面体网格 (退化成 10 节点四面体单元),且只有一个警告信息,是有少量单元质量超过警告线,但可以忽略。
相应的命令流如下:
!1.定义单元与材料属性
FINISH
/CLEAR
/filn, how_meshing-186, 1
/PREP7
ET,1,SOLID65
ET,2,LINK180
ET,3,SOLID186
keyopt,1,1,1 !不考虑形函数的附加项
keyopt,1,3,2 !无筋单元的压碎后行为
keyopt,1,7,1 !考虑开裂后的应力释放
!混凝土材料性质
fck=20.1 !C30混凝土抗压强度标准值
ftk=2.01 !C30混凝土抗拉强度标准值
mp,ex,1,fck*0.19/0.0002
mp,prxy,1,0.2
tb,concr,1
tbdata,,0.5,0.95,ftk,-1
tb,miso,1,,11
tbpt,,0.0002,fck*0.19
tbpt,,0.0004,fck*0.36
tbpt,,0.0006,fck*0.51
tbpt,,0.0008,fck*0.64
tbpt,,0.0010,fck*0.75
tbpt,,0.0012,fck*0.84
tbpt,,0.0014,fck*0.91
tbpt,,0.0016,fck*0.96
tbpt,,0.0018,fck*0.99
tbpt,,0.002,fck
tbpt,,0.0033,fck
!纵筋材料性质
mp,ex,2,2e5
mp,prxy,2,0.3
mp,dens,2,7.85e-9
tb,bkin,2
tbdata,,400,0
!箍筋材料性质
mp,ex,3,2e5
mp,prxy,3,0.3
mp,dens,3,7.85e-9
tb,bkin,3
tbdata,,300,0
!型钢材料性质
mp,ex,4,2.06e5
mp,prxy,4,0.3
mp,dens,4,7.85e-9
tb,bkin,4
tbdata,,235,0
r,1
pi=acos(-1)
r,2,18*18*pi/4
r,3,8*8*pi/4
!2.创建几何模型
BLC5,0,0,360,360,20
wpoff,0,0,20
BLC5,0,0,100,200,1120
BLC5,0,0,84,184,1120
wpoff,0,0,1120
BLC5,0,0,360,360,20
VSBV,2,3
wpoff,0,0,20
BLC5,0,0,360,360,1040
wpoff,0,0,200
wpro,,,-90.000000
wpoff,0,0,180
BLC5,0,0,400,280,1000
wpro,,,-180.000000
wpoff,0,0,360
BLC5,0,0,400,280,1000
VADD,ALL
!切分上部型钢柱
WPCSYS,-1
wpoff,0,0,20
VSBW,ALL
wpoff,0,0,80
VSBW,ALL
wpoff,0,0,200
VSBW,ALL
wpoff,0,0,840
VSBW,ALL
wpoff,0,0,20
VSBW,ALL
VSEL,S,LOC,Z,20,1140
wpro,,,-90.000000
wpoff,0,0,42
VSBW,ALL
wpoff,0,0,-84
VSBW,ALL
wpro,,-90.000000,
wpoff,0,0,92
VSBW,ALL
wpoff,0,0,-184
VSBW,ALL
!切分下部钢筋混凝土
ALLSEL,ALL
/REPLOT
WPCSYS,-1
VSEL,S,LOC,Z,1160,2200
!切出纵筋位置
wpoff,0,0,1200
VSBW,ALL
wpoff,0,0,160
VSBW,ALL
wpoff,0,0,160
VSBW,ALL
wpro,,-90.000000,
VSBW,ALL
wpoff,0,0,100
VSBW,ALL
wpoff,0,0,40
VSBW,ALL
wpoff,0,0,-240
VSBW,ALL
wpoff,0,0,-40
VSBW,ALL
WPCSYS,-1
wpro,,,-90.000000
VSBW,ALL
wpoff,0,0,140
VSBW,ALL
wpoff,0,0,-280
VSBW,ALL
!切出箍筋位置
wpoff,0,0,-1040
*DO,I,1,12
wpoff,0,0,80
VSBW,ALL
*ENDDO
wpoff,0,0,40
VSBW,ALL
wpoff,0,0,360
VSBW,ALL
wpoff,0,0,40
VSBW,ALL
*DO,I,1,11
wpoff,0,0,80
VSBW,ALL
*ENDDO
WPCSYS,-1
wpoff,0,0,1280
VSBW,ALL
wpoff,0,0,160
VSBW,ALL
wpoff,0,0,120
VSBW,ALL
wpoff,0,0,40
VSBW,ALL
*DO,I,1,7
wpoff,0,0,80
VSBW,ALL
*ENDDO
!3.划分网格
!划分下部钢筋混凝土网格
!划分纵筋
WPCSYS,-1
ALLSEL,ALL
LPLOT
LSEL,S,LOC,Y,100
LSEL,A,LOC,Y,-100
LSEL,R,LOC,Z,1200
LATT,2,2,2
LESIZE,ALL,40
LMESH,ALL
ALLSEL,ALL
LPLOT
LSEL,S,LOC,Y,100
LSEL,A,LOC,Y,-100
LSEL,R,LOC,Z,1360
LATT,2,2,2
LESIZE,ALL,40
LMESH,ALL
ALLSEL,ALL
LPLOT
LSEL,S,LOC,Y,100
LSEL,A,LOC,Y,-100
LSEL,R,LOC,Z,1520
LATT,2,2,2
LESIZE,ALL,40
LMESH,ALL
ALLSEL,ALL
LPLOT
LSEL,S,LOC,X,140
LSEL,A,LOC,X,0
LSEL,A,LOC,X,-140
LSEL,R,LOC,Y,140
LATT,2,2,2
LESIZE,ALL,40
LMESH,ALL
ALLSEL,ALL
LPLOT
LSEL,S,LOC,X,140
LSEL,A,LOC,X,-140
LSEL,R,LOC,Y,0
LATT,2,2,2
LESIZE,ALL,40
LMESH,ALL
ALLSEL,ALL
LPLOT
LSEL,S,LOC,X,140
LSEL,A,LOC,X,0
LSEL,A,LOC,X,-140
LSEL,R,LOC,Y,-140
LATT,2,2,2
LESIZE,ALL,40
LMESH,ALL
!划分箍筋
*DO,I,-1100,-220,80
ALLSEL,ALL
LPLOT
LSEL,S,LOC,Z,1200
LSEL,A,LOC,Z,1520
LSEL,R,LOC,Y,-100,100
LSEL,R,LOC,X,I
LATT,3,3,2
LESIZE,ALL,40
LMESH,ALL
ALLSEL,ALL
LPLOT
LSEL,S,LOC,Y,100
LSEL,A,LOC,Y,-100
LSEL,R,LOC,Z,1200,1520
LSEL,R,LOC,X,I
LATT,3,3,2
LESIZE,ALL,40
LMESH,ALL
*ENDDO
*DO,I,220,1100,80
ALLSEL,ALL
LPLOT
LSEL,S,LOC,Z,1200
LSEL,A,LOC,Z,1520
LSEL,R,LOC,Y,-100,100
LSEL,R,LOC,X,I
LATT,3,3,2
LESIZE,ALL,40
LMESH,ALL
ALLSEL,ALL
LPLOT
LSEL,S,LOC,Y,100
LSEL,A,LOC,Y,-100
LSEL,R,LOC,Z,1200,1520
LSEL,R,LOC,X,I
LATT,3,3,2
LESIZE,ALL,40
LMESH,ALL
*ENDDO
*DO,I,1200,2160,80
ALLSEL,ALL
LPLOT
LSEL,S,LOC,X,-140
LSEL,A,LOC,X,140
LSEL,R,LOC,Y,-140,140
LSEL,R,LOC,Z,I
LATT,3,3,2
LESIZE,ALL,40
LMESH,ALL
ALLSEL,ALL
LPLOT
LSEL,S,LOC,Y,-140
LSEL,A,LOC,Y,140
LSEL,R,LOC,X,-140,140
LSEL,R,LOC,Z,I
LATT,3,3,2
LESIZE,ALL,40
LMESH,ALL
*ENDDO
!*
/VIEW,1,,-1
/ANG,1
/REP,FAST
VPLOT
!划分混凝土
ALLSEL,ALL
MSHKEY,1
MSHAPE,0,3d
VPLOT
VSEL,S,LOC,Z,1160,2200
VATT,1,1,1
ESIZE,40
VMESH,ALL
!*
!划分上部型钢网格
!上部型钢都材料性质都为VATT,4,1,3
ALLSEL,ALL
/REPLOT
WPCSYS,-1
VSEL,S,LOC,Z,0,1160
!* 设置型钢部分的单元属性
TYPE, 3
MAT, 4
REAL,
ESYS, 0
SECNUM,
! 对上部型钢划分 FREE 网格
MSHKEY,0
MSHAPE,1,3d
FLST,5,26,6,ORDE,6
FITEM,5,1
FITEM,5,-6
FITEM,5,8
FITEM,5,-21
FITEM,5,28
FITEM,5,-33
VSEL, , , ,P51X
CHKMSH,'VOLU'
!*
VMESH,all
!*
allsel,all
eplot
save
!* |