！本文为天大读研期间编写文件，有问题可以交流，可供科研使用
/COM,Structural
/config,NRES,90000
/UNITS,SI
/prep7
gybh=0.010 !钢梁翼板厚度10mm
gybk=0.150 !钢梁翼板宽度150mm
gfbh=0.007 !钢梁腹板厚度7mm
gfbg=0.130 !钢梁腹板高度130mm
hntbh=0.075 !混凝土板厚度75mm
hntbk=0.600 !混凝土板宽度600mm
hntbch=6.000 !每节混凝土板长度6000mm
glch=6.200 !每节钢梁总长度600mm
xhntbh=0.060 !混凝土板厚度60mm
xhntbk=0.220 !混凝土板宽度220mm
xhntbch=1.200 !每节混凝土板长度1200mm
et,1,solid65,,,,,,3 !定义混凝土板单元
et,2,shell43 !定义钢梁单元
et,3,combin39,,,0 !定义弹簧单元combin39
et,4,solid45 !定义垫板单元
et,5,link8 !定义纵向钢筋单元
mp,ex,1,1.9095E10 !定义混凝土板的弹性模量N/mm2
mp,prxy,1,0.20 !定义混凝土板的泊松比
mp,dens,1,2.6e3 !定义混凝土板的密度
mp,ex,2,2.10E11 !定义钢梁及一级钢筋的弹性模量
mp,prxy,2,0.30 !定义钢梁及一级钢筋的泊松比
mp,dens,2,7.85e3 !定义钢梁及一级钢筋的密度
mp,ex,3,2.00E11 !定义二级钢筋的弹性模量
mp,prxy,3,0.3 !定义二级钢筋的泊松比
mp,dens,3,7.85e3 !定义二级钢筋的密度
mp,ex,4,2.02E11 !定义垫板材料特性
mp,prxy,4,0.3
tb,miso,1 !混凝土屈服准则
*do,e,0.0002,0.002,0.0001
tbpt,,e,20.1e6*(1000*e-500*500*e*e)
*enddo
tbpt,,0.0035,20.1e6 !该处将0.0033改为了0.0035
tbpl,miso,1
tb,concr,1 !混凝土破坏准则
tbdata,,0.6,0.8,2.01e6,-1
tb,bkin,2,1,2 !钢梁及一级钢筋的本构关系
tbdata,,235e6,10.3e6
tbpl,bkin,2
tb,bkin,3,1,2 !二级钢筋的本构关系
tbdata,,335e6,10.3e6 !335屈服应力，10代表屈服后的切线模量
tbpl,bkin,3
r,1,2,0.00745,90,90 !定义上混凝土板实常数
r,2,2,0.00673,90,90 !定义下混凝土板实常数
r,3,gybh,gybh,gybh,gybh !定义钢梁上下翼缘的厚度
r,4,gfbh,gfbh,gfbh,gfbh !定义钢梁腹板厚度
r,5,0,0,0.00002,12175,0.00004,16021 !定义弹簧的荷载变形曲线
rmore,0.00005,17492,0.00006,18789,0.00007,19956
rmore,0.00008,21022,0.00009,22006,0.0001,22921
rmore,0.0002,29831,0.0003,34610,0.0004,38313
rmore,0.0005,41338,0.0006,43887,0.0007,46078
rmore,0.0008,47989,0.0009,49673,0.001,51169
rmore,0.002,61124,0.003,63521
r,6 !垫板实常数
r,7,0.00004398 !上混凝土板中纵筋的横截面积（折合而来）
r,8,0.000050265 !下混凝土板中纵筋的横截面积
!***********************************************************************
!注：钢梁总长6200mm，跨度为3000+3000，在两端支座外各留出了100mm
blc4,0,0,glch,gfbg !钢梁腹板
csys,4
wprota,,-90
blc4,0,-gybk*0.5,glch,gybk*0.5 !钢梁下翼板之半(一次建成整个板时无法用aglue命令)
blc4,0,0,glch,gybk*0.5 !钢梁下翼板之半(一次建成整个板时无法用aglue命令)
wpoffs,0,0,gfbg
blc4,0,-gybk*0.5,glch,gybk*0.5 !钢梁上翼板之半(一次建成整个板时无法用aglue命令)
blc4,0,0,glch,gybk*0.5 !钢梁上翼板之半(一次建成整个板时无法用aglue命令)
wpoffs,0.100 !用工作平面切割梁体然后添加加劲肋
wprota,,,-90
asbw,all
wpoffs,,,-hntbch*0.5
asbw,all
wpoffs,,,-0.5*hntbch
asbw,all
nummrg,all
numcmp,all
*do,i,0,2,1
*do,j,0,1,1 !连线，为增加加劲肋作准备
l,15+j+6*i,17+j+6*i
*enddo
*enddo
al,28,50,21,16
al,31,16,25,51
al,33,52,7,3
al,34,3,11,53
al,44,54,27,19
al,47,19,41,55
asel,all
aglue,all !粘结各板形成钢梁
nummrg,all
numcmp,all
asel,s,loc,x,-gfbg !选中上下翼缘板，准备对钢梁翼板进行单元划分
asel,a,loc,x,0
aplot
aatt,2,3,2 !赋予钢梁上下翼缘板单元特性
lesize,31,0.0375$lesize,28,0.0375 !定义上下翼缘有限元网格尺寸
lesize,25,0.0375$lesize,21,0.0375
lesize,18,0.100$lesize,17,0.100
lesize,36,0.100$lesize,35,0.100
lesize,23,0.100$lesize,20,0.100
lesize,15,0.100$lesize,14,0.100
mshkey,1 !声明使用映射网格划分
mshape,0 !声明使用平面四边形网格划分
amesh,all
allsel,all
asel,s,loc,y,0 !选中腹板，准备对钢梁腹板进行单元划分
aplot
aatt,2,4,2
lesize,16,,,4 !定义腹板有限元网格尺寸
mshkey,1 !声明使用映射网格划分
mshape,0 !声明使用平面四边形网格划分
amesh,all
allsel,all
asel,s,loc,z,0
asel,a,loc,z,hntbch*0.5
asel,a,loc,z,hntbch
aplot
aatt,2,4,2
mshkey,1 !声明使用自由网格划分
mshape,0 !声明使用平面三角形网格划分
amesh,all
allsel,all
!*************************************************************************
wpoffs,gybh*0.5 !准备建立上侧混凝土板几何模型
blc4,0,-hntbk*0.5,hntbh,hntbk,hntbch
!**********************************************************************
!对上混凝土板进行映射划分
wpoffs,,0.0375
wprota,,90
vsbw,all
wpoffs,,,0.075
vsbw,all
vsel,all
vatt,1,1,1
lesize,78,,,3$lesize,87,,,1
lesize,80,,,3$lesize,60,,,2
lesize,66,0.100
mshkey,1 !声明使用映射网格划分
mshape,0,3D !声明使用平面四边形网格划分
vmesh,all
!esel,s,real,,1
!/DEVICE,VECTOR,1
!/TYPE,,BASIC !查看钢筋的方向
!/eshape,1
!allsel,all
!****************************************************
!加弹簧、耦合
!上混泥土板
!***************************
allsel,all
wpoffs,,,-0.0375
nsel,s,loc,z,-0.0375,0.0375
nsel,r,loc,x,-0.5*gybh,0
type,3 !设定弹簧单元
real,5
e,24,917
*do,i,0,13,1
e,330+2*i,977-2*i
*enddo
e,299,949
*do,i,0,13,1
e,545+2*i,947-2*i
*enddo
e,423,919
cp,1,uy,24,917
*do,i,-1,27,1
cp,3+i,uy,330+i,977-i
*enddo
cp,31,uy,299,949
*do,i,-1,27,1
cp,33+i,uy,545+i,947-i
*enddo
cp,61,uy,423,919
cp,62,uz,24,917
*do,i,-1,27,1
cp,64+i,uz,330+i,977-i
*enddo
cp,92,uz,299,949
*do,i,-1,27,1
cp,94+i,uz,545+i,947-i
*enddo
cp,122,uz,423,919
e,118,1600
*do,i,0,13,1
e,390+2*i,1660-2*i
*enddo
e,359,1632
*do,i,0,13,1
e,603+2*i,1630-2*i
*enddo
e,485,1589
!**********************
cp,123,uy,118,1600
*do,i,-1,27,1
cp,125+i,uy,390+i,1660-i
*enddo
cp,153,uy,359,1632
*do,i,-1,27,1
cp,155+i,uy,603+i,1630-i
*enddo
cp,183,uy,485,1589
cp,184,uz,118,1600
*do,i,-1,27,1
cp,186+i,uz,390+i,1660-i
*enddo
cp,214,uz,359,1632
*do,i,-1,27,1
cp,216+i,uz,603+i,1630-i
*enddo
cp,244,uz,485,1589
allsel,all
!搞定
!**********************************
!为上混凝土板增设纵向钢筋
nsel,s,loc,x,0.5*hntbh
type,5 !设定弹簧单元
real,7
*do,j,0,1,1
e,1232+j,1587-j
*do,i,0,57,1
e,1587-2*i-j,1585-2*i-j
*enddo
e,1471-j,1047-j
*enddo
e,918,1038
*do,i,0,57,1
e,1038-i,1037-i
*enddo
e,980,979
e,1602,1721
*do,i,0,57,1
e,1721+i,1722+i
*enddo
e,1779,1597
e,1787,2260
*do,i,0,57,1
e,2260-i,2259-i
*enddo
e,2202,1598
e,1788,2319
*do,i,0,57,1
e,2319-i,2318-i
*enddo
e,2261,1599
allsel,all
!**********************************
!为荷载和支座处添加弹性垫块
wpoffs,-0.5*gybh-gfbg,0.5*hntbch,-0.5*gybk
blc4,-0.5*gybh-0.010,-0.100,0.010,0.200,gybk
vsel,s,,,2
vatt,4,6,4
lesize,93,,,4$lesize,89,,,4$lesize,90,,,1
mshkey,1 !声明使用映射网格划分
mshape,0,3D !声明使用平面四边形网格划分
vmesh,all
nsel,s,loc,x,-0.5*gybh,0
nsel,r,loc,y,-0.100,0.100
cp,1000,uy,238,2325$cp,1001,uy,267,2346
cp,1002,uy,178,2347$cp,1003,uy,207,2348
cp,1004,uy,149,2332$cp,1005,uy,208,2328
cp,1006,uy,209,2368$cp,1007,uy,57,2365
cp,1008,uy,87,2362$cp,1009,uy,27,2334
cp,1010,uy,424,2326$cp,1011,uy,453,2343
cp,1012,uy,58,2344$cp,1013,uy,88,2345
cp,1014,uy,28,2331
cp,1005+35,uz,208,2328
cp,1006+35,uz,209,2368$cp,1007+35,uz,57,2365
cp,1008+35,uz,87,2362$cp,1009+35,uz,27,2334
cp,1005+70,ux,208,2328
cp,1006+70,ux,209,2368$cp,1007+70,ux,57,2365
cp,1008+70,ux,87,2362$cp,1009+70,ux,27,2334
allsel,all
nsel,s,loc,x,-0.5*gybh-0.010
nsel,r,loc,y,0
d,all,uy,0,,,,ux,uz
nsel,s,loc,x,0
nsel,r,loc,y,0.5*hntbch
d,all,uy,0,,,,uz
nsel,s,loc,x,0
nsel,r,loc,y,-0.5*hntbch
d,all,uy,0,,,,uz
allsel,all
/solu
antype,0
solcontrol,on
outres,all,5
nlgeom,off
kbc,0
nropt,full
autots,on
time,1
acel,,9.8
allsel,all
nsubst,5,5,5
solve
save
time,350
nsel,s,loc,x,0.5*gybh+gfbg+hntbh
nsel,r,loc,y,0.25*hntbch
f,902,fy,-175000,,1706,804
nsel,s,loc,x,0.5*gybh+gfbg+hntbh
nsel,r,loc,y,-0.25*hntbch
f,872,fy,-175000,,1676,804
allsel,all
nsubst,800,12800,800
solve
save
閱讀更多 肖河南6868 的文章
