请大侠帮忙看看，我的模型没法划分单元
我建了一个关于钢管混凝土节点的模型，钢牛腿与钢管都是用SOLID45单元，将钢牛腿与钢管建成体后，将两者粘结（VGLUE），模拟其焊接后共同受力的效果，但在划分网格的时候出现了一个问题，不管是用MAP、SWEEP还是FREE，都无法划分网格，请大侠帮我看看，急啊。/PREP7
ET,1,SOLID45 !钢板
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.0e11 !纵筋弹性模量
MPDATA,PRXY,1,,0.3 !纵筋泊松比
TB,BISO,1,1,2, !纵筋的本够关系为理想弹塑性
TBDATA,,335e6,,,,, !钢筋的屈服应力
ET,2,SOLID65 !混凝土
KEYOPT,2,1,0
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,34.5E9 !混凝土弹性模量
MPDATA,PRXY,2,,0.2 !混凝土泊松比
MPDATA,DENS,2,,2500 !混凝土密度
PressYield=32.4e6 !混凝土单轴抗压应力
TensionYield=2.64e6 !混凝土单轴抗拉应力
c2=0.5 !张开裂缝的剪切传递系数
c3=0.95 !笔合裂缝的剪切传递系数
!TB,CONCR,2
!TBDATA,,c2,c3,TensionYield,-1!W-W破坏准则
R,1,2,0.01, , , , , !混凝土实常数
R,2,,, !牛腿实常数
R,3,,,, !内侧加强板实常数
R,4,,,, !钢管实常数
K,1,0,-4.5
K,2,0.71,-4.5
K,3,0.75,-4.5
K,4,0.71,0.072
K,5,0.71,0.112
K,6,0.35,0.072
K,7,0.35,0.112
K,8,0,0.072
K,9,0,0.112
K,10,0.67,1.138
K,11,0.71,1.138
K,12,0,1.138
K,13,0.71,1.888
K,14,0.71,1.928
K,15,0.35,1.888
K,16,0.35,1.928
K,17,0.67,1.888
K,18,0,1.888
K,19,0,1.928
K,20,0.71,6.5
K,21,0.75,6.5
K,22,0,6.5
A,1,2,4,8
A,6,7,9,8
A,4,5,7,6
A,5,9,12,11
A,10,12,18,17
A,10,11,13,17
A,18,19,16,15
A,13,14,16,15
A,19,14,20,22
A,2,3,21,20
VROTAT,ALL,,,,,,1,22,, !旋转成体
WPOFFS,,0.112 !坐标轴沿竖向（Y轴）移动0.112
WPROTA,,-90 !沿X轴旋转，即将X-Y平面旋转为水平面
CYL4,,,0.75,,1,360
WPROTA,-1.146 !沿Z轴旋转1.146°
CSYS,4 !激活当前坐标系
K,1000,0.75,0,0
WPROTA,0.286 !沿Z轴旋转0.86°
CSYS,4 !激活当前坐标系
K,1001,1,0,0
L,1000,1001
L,68,73
LSBL,195,199,,DELETE,KEEP
LSBL,194,199,,DELETE,KEEP
AL,199,195,202,200
WPROTA,-22.718 !沿Z轴旋转23.578°
CSYS,4 !激活当前坐标系
K,1002,1,0,0
WPROTA,-8.653 !沿Z轴旋转32.231°
CSYS,4 !激活当前坐标系
K,1003,0.75,0,0
L,1002,1003
LSBL,201,194,,DELETE,KEEP
LSBL,203,194,,DELETE,KEEP
WPROTA,-2.752 !沿Z轴旋转34.983°
CSYS,4 !激活当前坐标系
K,1004,0.75,0,0
WPROTA,9.515 !沿Z轴旋转25.468°
CSYS,4 !激活当前坐标系
K,1005,1,0,0
L,1004,1005
LSBL,204,203,,DELETE,KEEP
LSBL,201,203,,DELETE,KEEP
AL,194,203,208,209
WPROTA,25.468 !沿Z轴旋转到零度角
ASEL,S,,,166
ARSYM,Y,ALL,,,,0,0
AADD,166,168
ALLSEL,ALL
AGEN,2,165,167,2,,,1.776
VEXT,165,,,,,-0.04
VEXT,167,,,,,0.43
VEXT,166,,,,,0.04
VEXT,168,,,,,-0.71
LGEN,2,205,206,,,,0.43,,1
LGEN,2,199,,,,,0.43,,1
!下面将压缩单元节点号和单元号
ALLSEL,ALL
NUMMRG,ALL
NUMCMP,ALL
AL,257,258,300,301,299,256
VEXT,210,,,,,0.03
VGEN,2,45,,,,,0.636
VEXT,169,,,,,1.776
VADD,44,46 !牛腿体相加
VADD,42,45
VSEL,S,,,44,48,4
VSYM,Y,ALL,,,,0,0
!复制其他三个方向的牛腿
VSEL,S,,,42
VSEL,A,,,44,45
VSEL,A,,,47,48
VSYM,X,ALL,,,,0,0
WPSTYL,,,,,,,,0 !返回总体坐标系
CSYS,1 !选用柱坐标系
VSEL,S,,,42
VSEL,A,,,44,45
VSEL,A,,,47,48
LOCAL,12,1,,,,,90 !定义局部坐标系，将XY平面旋转为水平面
VGEN,2,ALL,,,,90
CSDELE,12 !删除局部坐标系
WPSTYL,,,,,,,,0 !返回笛卡尔总体坐标系
VSEL,S,,,53,57
VSYM,Z,ALL,,,,0,0
!粘结牛腿
ALLSEL,ALL
VADD,42,44,45,47,48
VADD,53,54,55,56,57
VADD,46,49,50,51,52
VADD,58,59,60,61,62
!粘结钢管
VADD,10,20,30,40
!粘结1#环板与钢管
VGLUE,41,43,46
!粘结牛腿与钢管
VGLUE,10,42,44,45,63
!粘结牛腿与1#环板
VGLUE,20,41
VGLUE,20,43
VGLUE,10,30
VGLUE,41,43
VGLUE,40,42
VGLUE,40,10
VGLUE,41,46
VGLUE,40,42
!!!建混凝土梁，并镜像其他3个梁
ALLSEL,ALL
WPOFFS,,,-0.072 !坐标轴移动到原点
CYL4,,,0.75,,1.1,360 !建混凝土梁
RECTNG,0.8,4.0,0.5,-0.5
ASEL,S,,,234,
ARSYM,X,ALL,,,,0,0
CSYS,1
ASEL,S,,,234,
LOCAL,12,1,,,,,90 !定义局部坐标系，将XY平面旋转为水平面
AGEN,2,ALL,,,,90
CSDELE,12 !删除局部坐标系
WPSTYL,,,,,,,,0 !返回笛卡尔总体坐标系
ASEL,S,,,271,
ARSYM,Z,ALL,,,,0,0
ALLSEL,ALL
!把梁与环形梁相接
AADD,234,231
LFILLT,868,869,0.2
LFILLT,867,870,0.2
AL,485,278,738
AL,741,742,871
AADD,305,231,234
!把梁与环形梁相接
AADD,306,244
LFILLT,485,871,0.2
LFILLT,738,742,0.2
AL,788,872,858
AL,283,480,747
AADD,244,234,231
!把梁与环形梁相接
AADD,305,271
LFILLT,747,858,0.2
LFILLT,480,872,0.2
AL,859,738,861
AL,862,868,873
AADD,234,244,231
!把梁与环形梁相接
AADD,271,273
LFILLT,868,859,0.2
LFILLT,861,873,0.2
AL,863,865,485
AL,866,867,874
AADD,234,244,231
WPSTYL,,,,,,,,0 !返回总体坐标系
ALLSEL,ALL
VEXT,271,,,,2.0 !混凝土梁
ALLSEL,ALL !混凝土梁属性
VSEL,S,,,1,2
*DO,I,1,3,1
VSEL,A,,,1+10*I,2+10*I
*ENDDO
*DO,I,0,3,1
VSEL,A,,,4+10*I,5+10*I
VSEL,A,,,7+10*I,9+10*I
*ENDDO
VSEL,A,,,42
VATT,2,1,2
ALLSEL,ALL
VSEL,S,,,47 !钢管属性（实常数为4）
VATT,2,4,2
ALLSEL,ALL
VSEL,S,,,10,40,10 !牛腿属性（实常数为2）
VSEL,A,,,41,43,2
VATT,2,2,2
ALLSEL,ALL
VSEL,S,,,3,8,5 !钢管内2#内隔板和内测加劲环属性（实常数为3）
VSEL,A,,,13,18,5
VSEL,A,,,23,28,5
VSEL,A,,,33,38,5
VSEL,A,,,6,16,10
VSEL,A,,,26,36,10
VATT,2,3,2
!粘结梁与钢管
ALLSEL,ALL
VGLUE,42,47
!划分单元，体扫略网格
!无法划分单元，急啊
!下面将压缩单元节点号和单元号
ALLSEL,ALL
NUMMRG,ALL
NUMCMP,ALL
/SOLU
ANTYPE,0 !静力分析
PSTRES,ON
R,1,2,0.01, , , , , !混凝土实常数
R,2,,, !牛腿实常数
R,3,,,, !内侧加强板实常数
R,4,,,, !钢管实常数
ALLSEL,ALL
ASEL,S,LOC,Y,-4.5,-4.5
DA,ALL,ALL
FCUM,ADD !集中荷载叠加（默认加载方式为替代）
ALLSEL,ALL
NSEL,S,LOC,Y,6.5
*get,countnode,NODE,0,count
*IF,countnode,GT,0,THEN
F,ALL,FY,1*(-4.0805E7)/countnode !施加柱轴向压力
*ENDIF
/PBC,F,1 !用箭头显示集中荷载方向
ALLSEL,ALL
NSEL,S,LOC,X,4.0
*get,countnode,NODE,0,count
*IF,countnode,GT,0,THEN
F,ALL,FY,1*(-8.284E6)/countnode !施加第1根梁端剪力
*ENDIF
FCUM,ADD !集中荷载叠加（默认加载方式为替代）
/PBC,F,1 !用箭头显示集中荷载方向
ALLSEL,ALL
NSEL,S,LOC,X,-4.0
*get,countnode,NODE,0,count
*IF,countnode,GT,0,THEN
F,ALL,FY,1*(-7.7E6)/countnode !施加第2根梁端剪力
*ENDIF
FCUM,ADD !集中荷载叠加（默认加载方式为替代）
/PBC,F,1 !用箭头显示集中荷载方向
ALLSEL,ALL
NSEL,S,LOC,Z,4.0
*get,countnode,NODE,0,count
*IF,countnode,GT,0,THEN
F,ALL,FY,1*(-6.318E6)/countnode !施加第3根梁端剪力
*ENDIF
FCUM,ADD !集中荷载叠加（默认加载方式为替代）
/PBC,F,1 !用箭头显示集中荷载方向
ALLSEL,ALL
NSEL,S,LOC,Z,-4.0
*get,countnode,NODE,0,count
*IF,countnode,GT,0,THEN
F,ALL,FY,1*(-2.522E6)/countnode !施加第4根梁端剪力
*ENDIF
FCUM,ADD !集中荷载叠加（默认加载方式为替代）
/PBC,F,1 !用箭头显示集中荷载方向
ALLSEL,ALL
OUTPR,BASIC,ALL,
CNVTOL,U, ,0.05,2, ,
NSUBST,20,0,0
NLGEOM,1 !打开几何非线性
SOLVE
FINISH
页:
[1]