我的梁柱核心区网格划分怎么用映射划分，请大神指教
浏览：11829
/TITLE,TJ-5
finish
/clear
/prep7
et,1,65
et,2,8
et,3,45
!混凝土
MP,EX,1,30000
MP,PRXY,1,0.2
*SET,fc,23.4
*SET,ft,2.24
tb,concr,1
!1为材料号
tbdata,,0.5,0.95,ft,fc
TB,KINH,1,1,12,0
!Use 0 or leave blank to define stress -vs- total
TBPT,,0.00005,1.5
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.0010,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.0020,fc
tbpt,,0.0033,fc*0.85
!箍筋
MP,EX,2,2.1e5
MP,PRXY,2,0.25
TB,BKIN,2,1,2,1
TBDATA,,314.9
!纵筋
MP,EX,3,2.0e5
MP,PRXY,3,0.3
TB,BKIN,3,1,2,1
TBDATA,,378.6
!tb+tbdata
!纤维布
et,4,41
KEYOPT,4,1,2
*SET,RC,0.1
R,13,RC*2.02,RC*2.02,RC*2.02,RC*2.02
R,14,RC*3.001,RC*3.001,RC*3.001,RC*3.001
R,15,RC*3.01,RC*3.01,RC*3.01,RC*3.01
R,16,RC*2.001,RC*2.001,RC*2.001,RC*2.001
R,17,RC*2.01,RC*2.01,RC*2.01,RC*2.01
R,18,RC*2.0001,RC*2.0001,RC*2.0001,RC*2.0001
R,19,RC*2,RC*2,RC*2,RC*2
MPTEMP,,,,,,,,
!删除所有温度数据
MPTEMP,1,0
MPDATA,EX,4,,231012
!mptemp+mpdata
MPDATA,EY,4,,0.76
MPDATA,EZ,4,,0.76
MPDATA,PRXY,4,,0.34
MPDATA,PRYZ,4,,0.34
MPDATA,PRXZ,4,,0.34
MPDATA,GXY,4,,2070
MPDATA,GYZ,4,,2070
MPDATA,GXZ,4,,2070
TBDE,KINH,4,,,
!删除tb数据
TB,KINH,4,1,2,0
TBTEMP,0
!控制温度为零
TBPT,,0.0158,3650
TBPT,,0.01580001,0
!tbtemp+tbpt
!垫块
MP,EX,5,2.0e5
MP,PRXY,5,0.3
!核心区箍筋
MP,EX,6,2.1e5
MP,PRXY,6,0.25
TB,BKIN,6,1,2,1
TBDATA,,314.9
!核心区混凝土
MP,EX,7,42655
MP,PRXY,7,0.2
*SET,fc,44.9
*SET,ft,3.2
tb,concr,7
!7为材料号
tbdata,,0.5,0.95,ft,fc
!数据表中输入数据,Defines data for the data table
TB,KINH,7,1,12,0
tbpt,,0.00005,2.13275
!定义非线性数据点Defines a point on a nonlinear
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.0010,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.0020,fc
tbpt,,0.0033,fc*0.85
!材料8-角钢
mp,ex,8,206e3
mp,prxy,8,0.28
tb,biso,8,1,2,1
tbtemp,0
tbdata,1,215
tbdata,2,215
/axlab,x,strain
/axlab,y,stress
tbplot,biso,8,,,yes
r0=12
block,0,1000,0,250,0,125
!x=0~1000,y=0~250,z=0~125
wpoff,0,250,0
blc4,0,0,150,1200,100
!x=0+150=150,y=250+1200=1450,z=0~100
vadd,all
wpcsys,-1
wpoff,150,250
blc4,0,0,160,r0,100
blc4,0,r0,r0,160-r0,100
vsel,s,volu,,1,2
cm,jg,volu
!改动（斜加腋）
K,100,150+r0,410,0
K,101,150+r0,250+r0,0
!x=150,y=250,z=0
K,102,310,250+r0,0
!x=150+150+10=310,y=250,z=0
A,100,101,102
VEXT,24,,,,,125
!改动（方加腋）
wpcsys,-1
wpoff,0,250,100
blc4,0,0,150+r0,160,25
!x=0~150,y=250+160=410,z=100+25=125
blc4,150+r0,0,160-r0,r0,25
cmsel,u,jg
cm,gq,volu
vadd,all
wpcsys,0
wpoff,1000,0,0
blc4,0,0,50,250,125
!x=1000+50=1050,y=0~250,z=0~125
!划分
wpcsys,0
alls
wpoff,0,0,75
!z=75
VSBW,ALL
wpcsys,0
wpoff,0,0,100
!z=100
VSBW,ALL
wpcsys,0
wpro,,-90
wpoff,0,0,50
!y=50
VSBW,ALL
wpoff,0,0,150
!y=50+150=200
VSBW,ALL
wpoff,0,0,50
!y=200+50=250
VSBW,ALL
!改动（梁划分)
wpoff,0,0,160
!y=250+160=410(400+10)
VSBW,ALL
VSEL,S,LOC,Y,410,1450
wpoff,0,0,50
!y=410+50=460
VSBW,ALL
wpoff,0,0,150
!y=460+150=610
VSBW,ALL
wpoff,0,0,100
!y=610+100=710
VSBW,ALL
wpoff,0,0,100
!y=710+100=810
VSBW,ALL
wpoff,0,0,100
!y=810+100=910
VSBW,ALL
wpoff,0,0,100
!y=910+100=1010
VSBW,ALL
wpoff,0,0,190
!y=1010+190=1200
VSBW,ALL
!改动（柱划分）
ALLS
wpro,,,90
wpoff,0,0,100
!x=100
VSBW,ALL
wpoff,0,0,50
!x=100+50=150
VSBW,ALL
wpoff,0,0,160
!x=150+160=310
VSBW,ALL
VSEL,S,LOC,X,310,1000
wpoff,0,0,50
!x=310+50=360
VSBW,ALL
wpoff,0,0,150
!x=360+150=510
VSBW,ALL
wpoff,0,0,100
!x=510+100=610
VSBW,ALL
wpoff,0,0,100
!x=610+100=710
VSBW,ALL
wpoff,0,0,90
!x=710+90=800
VSBW,ALL
alls
!改动（角钢划分）
!wpcsys,-1
!wpro,,,90
!wpof,,,160
!x=150+10=160
!vsbw,all
!wpof,,,140
!x=160+140=300
!vsbw,all
!wpro,,-90
!wpof,,,260
!x=250+10=260
!vsbw,all
!wpof,,,140
!x=260+140=400
!vsbw,all
CSYS,0
wpcsys,-1
alls
!核心区
VSEL,S,LOC,X,-310,310
VSEL,R,LOC,Y,0,410
wpcsys,0
wpoff,310,250,125
!x=300+10=310,y=250,z=125
wpro,-135
wpro,,90
VSBW,ALL
ASBW,ALL
wpcsys,0
wpoff,150,410,125
!x=150,y=400+10=410,z=125
wpro,-135
wpro,,90
VSBW,ALL
ASBW,ALL
!wpcsys,0
!wpoff,-310,250,125
!wpro,135
!wpro,,90
!VSBW,ALL
!ASBW,ALL
wpcsys,0
wpoff,-150,410,125
!x=-150
wpro,135
wpro,,90
VSBW,ALL
ASBW,ALL
ALLSEL,ALL
*SET,pi,acos(-1)
RZ=16
RG=6
*SET,v11,0.25*pi*RZ*RZ*1.5/(125*50)!柱拉压区纵筋配筋率
*SET,v122,0.25*pi*RG*RG*4/(200*125)!柱箍筋加密区竖向
*SET,v133,0.25*pi*RG*RG*4/(200*50) !柱箍筋加密区水平向
*SET,v12,0.25*pi*RG*RG*6/(650*125)!柱箍筋非加密区竖向
*SET,v13,0.25*pi*RG*RG*6/(650*50) !柱箍筋非加密区水平向
*SET,v21,0.25*pi*RZ*RZ*1.5/(100*50)!梁拉压区纵筋配筋率
*SET,v222,0.25*pi*RG*RG*5/(250*125)!梁箍筋加密区竖向
*SET,v233,0.25*pi*RG*RG*5/(250*50) !梁箍筋加密区水平向
*SET,v22,0.25*pi*RG*RG*9/(950*125)!梁箍筋非加密区竖向
*SET,v23,0.25*pi*RG*RG*9/(950*50) !梁箍筋非加密区水平向
*SET,v0,0
r,1,3,v11,0,0,2,v133
rmore,0,90,2,v122,90,0
r,2,3,v11,0,0,2,v13
rmore,0,90,2,v12,90,0
r,3,3,v11,0,0,2,v21
rmore,90,0
r,4,2,v122,90,0
r,5,2,v12,90,0
r,6,3,V0,90,0
r,7,2,v222,0,0,2,v233
rmore,0,90,3,v21,90,0
r,8,2,v22,0,0,2,v23
rmore,0,90,3,v21,90,0
r,9,2,v222,0,0
r,10,2,v22,0,0
R,11,28.26
R,12, , , , , , ,
!赋值
vsel,s,loc,y,0,250
vsel,r,loc,x,800,1000
vsel,u,loc,y,50,200
vatt,1,1,1
!不变
!不变
vsel,s,loc,y,0,250
vsel,r,loc,x,150,800
vsel,u,loc,y,50,200
vatt,1,2,1
!不变
vsel,s,loc,y,0,250
vsel,r,loc,x,0,150
vsel,u,loc,y,50,200
vatt,1,3,1
vsel,s,loc,y,50,200
vsel,r,loc,x,800,1000
vatt,1,4,1
vsel,s,loc,y,50,200
vsel,r,loc,x,150,800
vatt,1,5,1
vsel,s,loc,y,50,200
vsel,r,loc,x,0,150
vatt,1,6,1
vsel,s,loc,y,1200,1450
vsel,r,loc,x,100,150
vatt,1,7,1
vsel,s,loc,y,250,1200
vsel,r,loc,x,100,150
vsel,r,loc,Z,0,100
vatt,1,8,1
vsel,s,loc,y,1200,1450
vsel,r,loc,x,0,100
vatt,1,9,1
vsel,s,loc,y,250,1200
vsel,r,loc,x,0,100
vsel,r,loc,Z,0,100
vatt,1,10,1
ALLSEL,ALL
VSEL,S,LOC,X,1000,1050
VATT,5,12,3
!v11
ALLSEL,ALL
VSEL,U,TYPE,,1
VSEL,U,TYPE,,3
VATT,7,12,1
!v12
ALLSEL,ALL
VSEL,U,mat,,1
VSEL,U,TYPE,,3
VSEL,r,loc,x,150,150+r0
VSEL,a,loc,y,250,250+r0
VSEL,u,loc,x,0,150
VSEL,u,loc,z,100,125
VATT,8,12,3
lSEL,S,LOC,X,0
lSEL,R,LOC,Y,50,200
lSEL,R,LOC,Z,0,100
lSEL,U,LOC,Z,0
lSEL,U,LOC,Z,75
LATT,6,11,2
!l1
LOCAL,11,0,1000,1000,1000,90,
LOCAL,12,0,1200,1200,1200,,,90
LOCAL,13,0,1400,1400,1400
LOCAL,14,0,800,800,800,-135
LOCAL,15,0,1600,1600,1600,-45
LOCAL,16,0,310,250,125,-135
LOCAL,17,0,-310,250,125,-45
LOCAL,18,0,310,250+r0,125,-135
LOCAL,19,0,-310,250+r0,125,-45
CSYS,0
ASEL,S,LOC,X,150
ASEL,R,LOC,Y,410,1010
aatt,4,13,4,11
ASEL,R,LOC,Y,910,1010
AGEN,2,ALL,,,500
ASEL,R,LOC,X,650
aatt,4,15,4,12
AGEN,2,ALL,,,-500
ASEL,S,LOC,Z,100
ASEL,R,LOC,Y,910,1010
aatt,4,15,4,13
ASEL,S,REAL,,15
ASEL,U,LOC,X,650
AGEN,2,ALL,,,,-200
ASEL,S,LOC,X,150
ASEL,R,LOC,Y,410,610
AGEN,2,ALL,,,500
ASEL,R,LOC,X,650
aatt,4,15,4,12
!a6
AGEN,2,ALL,,,-500
!a7
ASEL,S,LOC,Z,100
ASEL,R,LOC,Y,410,610
aatt,4,15,4,13
!a8
ASEL,S,REAL,,15
ASEL,R,LOC,X,650
ADELE,ALL
ASEL,S,LOC,X,310,710
ASEL,R,LOC,Y,250
ASEL,R,LOC,Z,0,100
aatt,4,13,4,13
ASEL,S,LOC,X,610,710
ASEL,R,LOC,Y,250
AGEN,2,ALL,,,,500
ASEL,R,LOC,Y,750
aatt,4,14,4,12
AGEN,2,ALL,,,,-500
ASEL,S,LOC,X,310,510
ASEL,R,LOC,Y,250
AGEN,2,ALL,,,,500
ASEL,R,LOC,Y,750
aatt,4,14,4,12
!a11
AGEN,2,ALL,,,,-500
!a12
ASEL,S,REAL,,14
ASEL,R,LOC,Y,750
ADELE,ALL
ASEL,S,LOC,X,610,710
ASEL,R,LOC,Z,125
aatt,4,14,4,11
!a13
ASEL,S,LOC,X,310,510
ASEL,R,LOC,Z,125
aatt,4,14,4,11
!a14
ASEL,S,LOC,X,610,710
ASEL,R,LOC,Y,0
aatt,4,14,4,12
!a15
ASEL,S,LOC,X,310,510
ASEL,R,LOC,Y,0
aatt,4,14,4,12
!a16
alls
csys,18
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,150,310
!常规选择面布
AGEN,2,ALL,,,600,600
ASEL,R,LOC,X,310,1100
cm,a1,area
aatt,4,16,4,12
AGEN,2,ALL,,,-600,-600
ASEL,S,LOC,X,150,310
ASEL,R,LOC,Y,0,410
ASEL,R,LOC,Z,125
aatt,4,16,4,11
!a19
ASEL,S,LOC,X,150,310
ASEL,R,LOC,Y,0
aatt,4,16,4,12
!a20
cmsel,s,a1
aatt,4,17,4,12
!a21
AGEN,2,ALL,,,-600,-600 !a22
ASEL,S,LOC,X,0,310
ASEL,R,LOC,Y,250,410
ASEL,R,LOC,Z,125
AGEN,2,ALL,,,,,500
!a23
ASEL,R,LOC,Z,625
aatt,4,17,4,13
!a24
AGEN,2,ALL,,,,,-500
!a25
ASEL,S,REAL,,17
ASEL,R,LOC,Z,625
ADELE,ALL
ASEL,S,REAL,,17
ASEL,R,LOC,X,410,1000
ADELE,ALL
LSEL,S,TYPE,,2
LSYMM,X,ALL
asel,s,mat,,4
ARSYM,X,ALL
vsel,s,all
VSYMM,X,ALL
!变动（斜面纤维布两侧）
alls
LOCAL,18,0,310,250+r0,125,-135
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,150,310
!常规选择面布
AGEN,2,ALL,,,600,600
ASEL,R,LOC,X,310,1100
cm,a1,area
aatt,4,13,4,15
AGEN,2,ALL,,,-600,-600
cmsel,s,a1
ADELE,A1
LOCAL,19,0,-310,250+r0,125,-45
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,-150,-310
!常规选择面布
AGEN,2,ALL,,,-600,600
ASEL,R,LOC,X,-310,-1100
cm,a1,area
aatt,4,13,4,14
!a33
AGEN,2,ALL,,,600,-600
!a28
cmsel,s,a1
ADELE,A1
csys,18
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,150,310
!常规选择面布
AGEN,2,ALL,,,600,600
ASEL,R,LOC,X,310,1100
cm,a1,area
aatt,4,18,4,12
!a29
AGEN,2,ALL,,,-600,-600
!a30
cmsel,s,a1
ADELE,A1
csys,19
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,-150,-310
!常规选择面布
AGEN,2,ALL,,,-600,600
ASEL,R,LOC,X,-310,-1100
cm,a1,area
aatt,4,19,4,12
!a31
AGEN,2,ALL,,,600,-600
!a32
cmsel,s,a1
ADELE,A1
ASEL,S,LOC,X,-260,60
ASEL,R,LOC,Y,0
AGEN,2,ALL,,,,-500
ASEL,R,LOC,Y,-500
aatt,4,18,4,12
!a34
AGEN,2,ALL,,,,500
!a35
ASEL,R,LOC,Y,-500
ADELE,ALL
ASEL,S,LOC,X,-60,260
ASEL,R,LOC,Y,0
AGEN,2,ALL,,,,-500
ASEL,R,LOC,Y,-500
aatt,4,19,4,12
!a36
AGEN,2,ALL,,,,500
!a37
ASEL,R,LOC,Y,-500
ADELE,ALL
ALLS
csys,16
ASEL,S,LOC,X,0,1000
ASEL,R,LOC,Y,0,-230 !在local16下的xyz坐标面
ASEL,R,LOC,Z,0
AGEN,2,ALL,,,,,500 !agen,只能用于宏观坐标系，不涉及local
ASEL,R,LOC,Z,500
aatt,4,18,4,16
!a38
AGEN,2,ALL,,,,,-500
!a39
ASEL,R,LOC,Z,500
ADELE,ALL
csys,17
ASEL,S,LOC,X,0,1000
ASEL,R,LOC,Y,0,230
ASEL,R,LOC,Z,0
AGEN,2,ALL,,,,,500
ASEL,R,LOC,Z,500
aatt,4,19,4,17
!a40
AGEN,2,ALL,,,,,-500
!a41
ASEL,R,LOC,Z,500
ADELE,ALL
CSYS,0
wpcsys,-1
ALLS
VSEL,S,LOC,Y,1200,1450
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
LSEL,S,TYPE,,2
MSHKEY,1
esize,50
LMESH,ALL
VSEL,S,LOC,Y,0,410
VSEL,R,LOC,X,-310,310
VSEL,INVE
VSEL,U,LOC,Y,1200,1450
VSEL,R,LOC,y,410,1110
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,s,LOC,x,-1010,-310
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,s,LOC,x,310,1010
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,s,LOC,x,1010,1050
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,s,LOC,x,-1050,-1010
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,S,LOC,Y,0,410
VSEL,R,LOC,X,-310,310
vsel,u,mat,,8
vsel,u,mat,,7
esize,50
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
vsel,s,mat,,7
esize,50
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
vsel,s,mat,,8
esize,50
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
ASEL,S,TYPE,,4
asel,r,loc,x,-310,310
asel,r,loc,y,0,410
cm,frp1,volu
mshape,1,2d
mshkey,0
AMESH,ALL
ASEL,S,TYPE,,4
cmsel,u,frp1,volu
esize,50
mshape,0,2d
mshkey,1
AMESH,ALL
allsel,all
nummrg,elem
nsle
nummrg,node
ksln
nummrg,kp
NUMCMP,ALL
alls
csys,0
nsel,s,loc,X,-150
nsel,r,loc,y,1300,1450
CP,1,UX,ALL
alls
!位移加载
FINISH
/SOL
nsel,s,loc,X,1050
D,ALL, , , , , ,UX,UY,UZ, , ,
nsel,s,loc,X,-1050
D,ALL, , , , , ,UY,UZ, , ,
ALLS
ASEL,S,LOC,Z,0$DA,ALL,SYMM
nsel,s,loc,X,-1050
F,ALL,FX,3704
alls
EPLOT
/SOL
anty,static
eqslv,spar,,0,
time,1
outres,all,all
nsubst,20
cnvtol,u,,0.1,0,,
lswrite,1
TIME,2
nsubst,100
d,133,ux,20
!位移加载
nropt,full,,
OUTRES,ALL,ALL
LNSRCH,on
lswrite,2
lssolve,1,2
/POST26
nsol,2,133,u,x
rforce,3,133,f,x
prod,2,2,,,,,,1
prod,3,3,,,,,,0.002
/axlab,x,ux of 133(mm)
/axlab,y,fx(N)
xvar,2
plvar,3
/POST26
lines,9999
prvar,2,3
finish
/clear
/prep7
et,1,65
et,2,8
et,3,45
!混凝土
MP,EX,1,30000
MP,PRXY,1,0.2
*SET,fc,23.4
*SET,ft,2.24
tb,concr,1
!1为材料号
tbdata,,0.5,0.95,ft,fc
TB,KINH,1,1,12,0
!Use 0 or leave blank to define stress -vs- total
TBPT,,0.00005,1.5
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.0010,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.0020,fc
tbpt,,0.0033,fc*0.85
!箍筋
MP,EX,2,2.1e5
MP,PRXY,2,0.25
TB,BKIN,2,1,2,1
TBDATA,,314.9
!纵筋
MP,EX,3,2.0e5
MP,PRXY,3,0.3
TB,BKIN,3,1,2,1
TBDATA,,378.6
!tb+tbdata
!纤维布
et,4,41
KEYOPT,4,1,2
*SET,RC,0.1
R,13,RC*2.02,RC*2.02,RC*2.02,RC*2.02
R,14,RC*3.001,RC*3.001,RC*3.001,RC*3.001
R,15,RC*3.01,RC*3.01,RC*3.01,RC*3.01
R,16,RC*2.001,RC*2.001,RC*2.001,RC*2.001
R,17,RC*2.01,RC*2.01,RC*2.01,RC*2.01
R,18,RC*2.0001,RC*2.0001,RC*2.0001,RC*2.0001
R,19,RC*2,RC*2,RC*2,RC*2
MPTEMP,,,,,,,,
!删除所有温度数据
MPTEMP,1,0
MPDATA,EX,4,,231012
!mptemp+mpdata
MPDATA,EY,4,,0.76
MPDATA,EZ,4,,0.76
MPDATA,PRXY,4,,0.34
MPDATA,PRYZ,4,,0.34
MPDATA,PRXZ,4,,0.34
MPDATA,GXY,4,,2070
MPDATA,GYZ,4,,2070
MPDATA,GXZ,4,,2070
TBDE,KINH,4,,,
!删除tb数据
TB,KINH,4,1,2,0
TBTEMP,0
!控制温度为零
TBPT,,0.0158,3650
TBPT,,0.01580001,0
!tbtemp+tbpt
!垫块
MP,EX,5,2.0e5
MP,PRXY,5,0.3
!核心区箍筋
MP,EX,6,2.1e5
MP,PRXY,6,0.25
TB,BKIN,6,1,2,1
TBDATA,,314.9
!核心区混凝土
MP,EX,7,42655
MP,PRXY,7,0.2
*SET,fc,44.9
*SET,ft,3.2
tb,concr,7
!7为材料号
tbdata,,0.5,0.95,ft,fc
!数据表中输入数据,Defines data for the data table
TB,KINH,7,1,12,0
tbpt,,0.00005,2.13275
!定义非线性数据点Defines a point on a nonlinear
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.0010,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.0020,fc
tbpt,,0.0033,fc*0.85
!材料8-角钢
mp,ex,8,206e3
mp,prxy,8,0.28
tb,biso,8,1,2,1
tbtemp,0
tbdata,1,215
tbdata,2,215
/axlab,x,strain
/axlab,y,stress
tbplot,biso,8,,,yes
r0=12
block,0,1000,0,250,0,125
!x=0~1000,y=0~250,z=0~125
wpoff,0,250,0
blc4,0,0,150,1200,100
!x=0+150=150,y=250+1200=1450,z=0~100
vadd,all
wpcsys,-1
wpoff,150,250
blc4,0,0,160,r0,100
blc4,0,r0,r0,160-r0,100
vsel,s,volu,,1,2
cm,jg,volu
!改动（斜加腋）
K,100,150+r0,410,0
K,101,150+r0,250+r0,0
!x=150,y=250,z=0
K,102,310,250+r0,0
!x=150+150+10=310,y=250,z=0
A,100,101,102
VEXT,24,,,,,125
!改动（方加腋）
wpcsys,-1
wpoff,0,250,100
blc4,0,0,150+r0,160,25
!x=0~150,y=250+160=410,z=100+25=125
blc4,150+r0,0,160-r0,r0,25
cmsel,u,jg
cm,gq,volu
vadd,all
wpcsys,0
wpoff,1000,0,0
blc4,0,0,50,250,125
!x=1000+50=1050,y=0~250,z=0~125
!划分
wpcsys,0
alls
wpoff,0,0,75
!z=75
VSBW,ALL
wpcsys,0
wpoff,0,0,100
!z=100
VSBW,ALL
wpcsys,0
wpro,,-90
wpoff,0,0,50
!y=50
VSBW,ALL
wpoff,0,0,150
!y=50+150=200
VSBW,ALL
wpoff,0,0,50
!y=200+50=250
VSBW,ALL
!改动（梁划分)
wpoff,0,0,160
!y=250+160=410(400+10)
VSBW,ALL
VSEL,S,LOC,Y,410,1450
wpoff,0,0,50
!y=410+50=460
VSBW,ALL
wpoff,0,0,150
!y=460+150=610
VSBW,ALL
wpoff,0,0,100
!y=610+100=710
VSBW,ALL
wpoff,0,0,100
!y=710+100=810
VSBW,ALL
wpoff,0,0,100
!y=810+100=910
VSBW,ALL
wpoff,0,0,100
!y=910+100=1010
VSBW,ALL
wpoff,0,0,190
!y=1010+190=1200
VSBW,ALL
!改动（柱划分）
ALLS
wpro,,,90
wpoff,0,0,100
!x=100
VSBW,ALL
wpoff,0,0,50
!x=100+50=150
VSBW,ALL
wpoff,0,0,160
!x=150+160=310
VSBW,ALL
VSEL,S,LOC,X,310,1000
wpoff,0,0,50
!x=310+50=360
VSBW,ALL
wpoff,0,0,150
!x=360+150=510
VSBW,ALL
wpoff,0,0,100
!x=510+100=610
VSBW,ALL
wpoff,0,0,100
!x=610+100=710
VSBW,ALL
wpoff,0,0,90
!x=710+90=800
VSBW,ALL
alls
!改动（角钢划分）
!wpcsys,-1
!wpro,,,90
!wpof,,,160
!x=150+10=160
!vsbw,all
!wpof,,,140
!x=160+140=300
!vsbw,all
!wpro,,-90
!wpof,,,260
!x=250+10=260
!vsbw,all
!wpof,,,140
!x=260+140=400
!vsbw,all
CSYS,0
wpcsys,-1
alls
!核心区
VSEL,S,LOC,X,-310,310
VSEL,R,LOC,Y,0,410
wpcsys,0
wpoff,310,250,125
!x=300+10=310,y=250,z=125
wpro,-135
wpro,,90
VSBW,ALL
ASBW,ALL
wpcsys,0
wpoff,150,410,125
!x=150,y=400+10=410,z=125
wpro,-135
wpro,,90
VSBW,ALL
ASBW,ALL
!wpcsys,0
!wpoff,-310,250,125
!wpro,135
!wpro,,90
!VSBW,ALL
!ASBW,ALL
wpcsys,0
wpoff,-150,410,125
!x=-150
wpro,135
wpro,,90
VSBW,ALL
ASBW,ALL
ALLSEL,ALL
*SET,pi,acos(-1)
RZ=16
RG=6
*SET,v11,0.25*pi*RZ*RZ*1.5/(125*50)!柱拉压区纵筋配筋率
*SET,v122,0.25*pi*RG*RG*4/(200*125)!柱箍筋加密区竖向
*SET,v133,0.25*pi*RG*RG*4/(200*50) !柱箍筋加密区水平向
*SET,v12,0.25*pi*RG*RG*6/(650*125)!柱箍筋非加密区竖向
*SET,v13,0.25*pi*RG*RG*6/(650*50) !柱箍筋非加密区水平向
*SET,v21,0.25*pi*RZ*RZ*1.5/(100*50)!梁拉压区纵筋配筋率
*SET,v222,0.25*pi*RG*RG*5/(250*125)!梁箍筋加密区竖向
*SET,v233,0.25*pi*RG*RG*5/(250*50) !梁箍筋加密区水平向
*SET,v22,0.25*pi*RG*RG*9/(950*125)!梁箍筋非加密区竖向
*SET,v23,0.25*pi*RG*RG*9/(950*50) !梁箍筋非加密区水平向
*SET,v0,0
r,1,3,v11,0,0,2,v133
rmore,0,90,2,v122,90,0
r,2,3,v11,0,0,2,v13
rmore,0,90,2,v12,90,0
r,3,3,v11,0,0,2,v21
rmore,90,0
r,4,2,v122,90,0
r,5,2,v12,90,0
r,6,3,V0,90,0
r,7,2,v222,0,0,2,v233
rmore,0,90,3,v21,90,0
r,8,2,v22,0,0,2,v23
rmore,0,90,3,v21,90,0
r,9,2,v222,0,0
r,10,2,v22,0,0
R,11,28.26
R,12, , , , , , ,
!赋值
vsel,s,loc,y,0,250
vsel,r,loc,x,800,1000
vsel,u,loc,y,50,200
vatt,1,1,1
!不变
!不变
vsel,s,loc,y,0,250
vsel,r,loc,x,150,800
vsel,u,loc,y,50,200
vatt,1,2,1
!不变
vsel,s,loc,y,0,250
vsel,r,loc,x,0,150
vsel,u,loc,y,50,200
vatt,1,3,1
vsel,s,loc,y,50,200
vsel,r,loc,x,800,1000
vatt,1,4,1
vsel,s,loc,y,50,200
vsel,r,loc,x,150,800
vatt,1,5,1
vsel,s,loc,y,50,200
vsel,r,loc,x,0,150
vatt,1,6,1
vsel,s,loc,y,1200,1450
vsel,r,loc,x,100,150
vatt,1,7,1
vsel,s,loc,y,250,1200
vsel,r,loc,x,100,150
vsel,r,loc,Z,0,100
vatt,1,8,1
vsel,s,loc,y,1200,1450
vsel,r,loc,x,0,100
vatt,1,9,1
vsel,s,loc,y,250,1200
vsel,r,loc,x,0,100
vsel,r,loc,Z,0,100
vatt,1,10,1
ALLSEL,ALL
VSEL,S,LOC,X,1000,1050
VATT,5,12,3
!v11
ALLSEL,ALL
VSEL,U,TYPE,,1
VSEL,U,TYPE,,3
VATT,7,12,1
!v12
ALLSEL,ALL
VSEL,U,mat,,1
VSEL,U,TYPE,,3
VSEL,r,loc,x,150,150+r0
VSEL,a,loc,y,250,250+r0
VSEL,u,loc,x,0,150
VSEL,u,loc,z,100,125
VATT,8,12,3
lSEL,S,LOC,X,0
lSEL,R,LOC,Y,50,200
lSEL,R,LOC,Z,0,100
lSEL,U,LOC,Z,0
lSEL,U,LOC,Z,75
LATT,6,11,2
!l1
LOCAL,11,0,1000,1000,1000,90,
LOCAL,12,0,1200,1200,1200,,,90
LOCAL,13,0,1400,1400,1400
LOCAL,14,0,800,800,800,-135
LOCAL,15,0,1600,1600,1600,-45
LOCAL,16,0,310,250,125,-135
LOCAL,17,0,-310,250,125,-45
LOCAL,18,0,310,250+r0,125,-135
LOCAL,19,0,-310,250+r0,125,-45
CSYS,0
ASEL,S,LOC,X,150
ASEL,R,LOC,Y,410,1010
aatt,4,13,4,11
ASEL,R,LOC,Y,910,1010
AGEN,2,ALL,,,500
ASEL,R,LOC,X,650
aatt,4,15,4,12
AGEN,2,ALL,,,-500
ASEL,S,LOC,Z,100
ASEL,R,LOC,Y,910,1010
aatt,4,15,4,13
ASEL,S,REAL,,15
ASEL,U,LOC,X,650
AGEN,2,ALL,,,,-200
ASEL,S,LOC,X,150
ASEL,R,LOC,Y,410,610
AGEN,2,ALL,,,500
ASEL,R,LOC,X,650
aatt,4,15,4,12
!a6
AGEN,2,ALL,,,-500
!a7
ASEL,S,LOC,Z,100
ASEL,R,LOC,Y,410,610
aatt,4,15,4,13
!a8
ASEL,S,REAL,,15
ASEL,R,LOC,X,650
ADELE,ALL
ASEL,S,LOC,X,310,710
ASEL,R,LOC,Y,250
ASEL,R,LOC,Z,0,100
aatt,4,13,4,13
ASEL,S,LOC,X,610,710
ASEL,R,LOC,Y,250
AGEN,2,ALL,,,,500
ASEL,R,LOC,Y,750
aatt,4,14,4,12
AGEN,2,ALL,,,,-500
ASEL,S,LOC,X,310,510
ASEL,R,LOC,Y,250
AGEN,2,ALL,,,,500
ASEL,R,LOC,Y,750
aatt,4,14,4,12
!a11
AGEN,2,ALL,,,,-500
!a12
ASEL,S,REAL,,14
ASEL,R,LOC,Y,750
ADELE,ALL
ASEL,S,LOC,X,610,710
ASEL,R,LOC,Z,125
aatt,4,14,4,11
!a13
ASEL,S,LOC,X,310,510
ASEL,R,LOC,Z,125
aatt,4,14,4,11
!a14
ASEL,S,LOC,X,610,710
ASEL,R,LOC,Y,0
aatt,4,14,4,12
!a15
ASEL,S,LOC,X,310,510
ASEL,R,LOC,Y,0
aatt,4,14,4,12
!a16
alls
csys,18
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,150,310
!常规选择面布
AGEN,2,ALL,,,600,600
ASEL,R,LOC,X,310,1100
cm,a1,area
aatt,4,16,4,12
AGEN,2,ALL,,,-600,-600
ASEL,S,LOC,X,150,310
ASEL,R,LOC,Y,0,410
ASEL,R,LOC,Z,125
aatt,4,16,4,11
!a19
ASEL,S,LOC,X,150,310
ASEL,R,LOC,Y,0
aatt,4,16,4,12
!a20
cmsel,s,a1
aatt,4,17,4,12
!a21
AGEN,2,ALL,,,-600,-600 !a22
ASEL,S,LOC,X,0,310
ASEL,R,LOC,Y,250,410
ASEL,R,LOC,Z,125
AGEN,2,ALL,,,,,500
!a23
ASEL,R,LOC,Z,625
aatt,4,17,4,13
!a24
AGEN,2,ALL,,,,,-500
!a25
ASEL,S,REAL,,17
ASEL,R,LOC,Z,625
ADELE,ALL
ASEL,S,REAL,,17
ASEL,R,LOC,X,410,1000
ADELE,ALL
LSEL,S,TYPE,,2
LSYMM,X,ALL
asel,s,mat,,4
ARSYM,X,ALL
vsel,s,all
VSYMM,X,ALL
!变动（斜面纤维布两侧）
alls
LOCAL,18,0,310,250+r0,125,-135
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,150,310
!常规选择面布
AGEN,2,ALL,,,600,600
ASEL,R,LOC,X,310,1100
cm,a1,area
aatt,4,13,4,15
AGEN,2,ALL,,,-600,-600
cmsel,s,a1
ADELE,A1
LOCAL,19,0,-310,250+r0,125,-45
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,-150,-310
!常规选择面布
AGEN,2,ALL,,,-600,600
ASEL,R,LOC,X,-310,-1100
cm,a1,area
aatt,4,13,4,14
!a33
AGEN,2,ALL,,,600,-600
!a28
cmsel,s,a1
ADELE,A1
csys,18
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,150,310
!常规选择面布
AGEN,2,ALL,,,600,600
ASEL,R,LOC,X,310,1100
cm,a1,area
aatt,4,18,4,12
!a29
AGEN,2,ALL,,,-600,-600
!a30
cmsel,s,a1
ADELE,A1
csys,19
asel,s,loc,x,0
csys,0
wpcsys,-1
asel,r,loc,x,-150,-310
!常规选择面布
AGEN,2,ALL,,,-600,600
ASEL,R,LOC,X,-310,-1100
cm,a1,area
aatt,4,19,4,12
!a31
AGEN,2,ALL,,,600,-600
!a32
cmsel,s,a1
ADELE,A1
ASEL,S,LOC,X,-260,60
ASEL,R,LOC,Y,0
AGEN,2,ALL,,,,-500
ASEL,R,LOC,Y,-500
aatt,4,18,4,12
!a34
AGEN,2,ALL,,,,500
!a35
ASEL,R,LOC,Y,-500
ADELE,ALL
ASEL,S,LOC,X,-60,260
ASEL,R,LOC,Y,0
AGEN,2,ALL,,,,-500
ASEL,R,LOC,Y,-500
aatt,4,19,4,12
!a36
AGEN,2,ALL,,,,500
!a37
ASEL,R,LOC,Y,-500
ADELE,ALL
ALLS
csys,16
ASEL,S,LOC,X,0,1000
ASEL,R,LOC,Y,0,-230 !在local16下的xyz坐标面
ASEL,R,LOC,Z,0
AGEN,2,ALL,,,,,500 !agen,只能用于宏观坐标系，不涉及local
ASEL,R,LOC,Z,500
aatt,4,18,4,16
!a38
AGEN,2,ALL,,,,,-500
!a39
ASEL,R,LOC,Z,500
ADELE,ALL
csys,17
ASEL,S,LOC,X,0,1000
ASEL,R,LOC,Y,0,230
ASEL,R,LOC,Z,0
AGEN,2,ALL,,,,,500
ASEL,R,LOC,Z,500
aatt,4,19,4,17
!a40
AGEN,2,ALL,,,,,-500
!a41
ASEL,R,LOC,Z,500
ADELE,ALL
CSYS,0
wpcsys,-1
ALLS
VSEL,S,LOC,Y,1200,1450
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
LSEL,S,TYPE,,2
MSHKEY,1
esize,50
LMESH,ALL
VSEL,S,LOC,Y,0,410
VSEL,R,LOC,X,-310,310
VSEL,INVE
VSEL,U,LOC,Y,1200,1450
VSEL,R,LOC,y,410,1110
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,s,LOC,x,-1010,-310
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,s,LOC,x,310,1010
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,s,LOC,x,1010,1050
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,s,LOC,x,-1050,-1010
MSHAPE,0,3D
MSHKEY,1
ESIZE,50
VMESH,ALL
VSEL,S,LOC,Y,0,410
VSEL,R,LOC,X,-310,310
vsel,u,mat,,8
vsel,u,mat,,7
esize,50
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
vsel,s,mat,,7
esize,50
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
vsel,s,mat,,8
esize,50
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
ASEL,S,TYPE,,4
asel,r,loc,x,-310,310
asel,r,loc,y,0,410
cm,frp1,volu
mshape,1,2d
mshkey,0
AMESH,ALL
ASEL,S,TYPE,,4
cmsel,u,frp1,volu
esize,50
mshape,0,2d
mshkey,1
AMESH,ALL
allsel,all
nummrg,elem
nsle
nummrg,node
ksln
nummrg,kp
NUMCMP,ALL
alls
csys,0
nsel,s,loc,X,-150
nsel,r,loc,y,1300,1450
CP,1,UX,ALL
alls
!位移加载
FINISH
/SOL
nsel,s,loc,X,1050
D,ALL, , , , , ,UX,UY,UZ, , ,
nsel,s,loc,X,-1050
D,ALL, , , , , ,UY,UZ, , ,
ALLS
ASEL,S,LOC,Z,0$DA,ALL,SYMM
nsel,s,loc,X,-1050
F,ALL,FX,3704
alls
EPLOT
/SOL
anty,static
eqslv,spar,,0,
time,1
outres,all,all
nsubst,20
cnvtol,u,,0.1,0,,
lswrite,1
TIME,2
nsubst,100
d,133,ux,20
!位移加载
nropt,full,,
OUTRES,ALL,ALL
LNSRCH,on
lswrite,2
lssolve,1,2
/POST26
nsol,2,133,u,x
rforce,3,133,f,x
prod,2,2,,,,,,1
prod,3,3,,,,,,0.002
/axlab,x,ux of 133(mm)
/axlab,y,fx(N)
xvar,2
plvar,3
/POST26
lines,9999
prvar,2,3
