！本文為天大讀研期間編寫文件，有問題可以交流，可供科研使用
/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 的文章
