1、/TITLE,Concrete-Filled Tube Arch Bridge/PREP7ET,1,82 !辅助单元类型CYL4,0.5,0.5,0.5,0.486 !建立一个环面SMRTSIZE,5 !网格划分AMESH,ALLSECWRITE,STEEL_TUBE,SECT,1 !截面存盘(文件名:STEEL_TUBE;文件拓展名:SECT;单元类型属性:1)SECTYPE,1,BEAM,MESH, !截面类型和ID定义SECOFFSET,CENT, !截面偏移,CENT偏移到质心SECREAD,STEEL_TUBE,SECT,MESH !读入截面ASEL,ALL !选择所有面ACLEAR
2、,ALLADELE,ALL,1 !将图形及面、线、点全部删除(只保留截面文件)CYL4,0.5,0.5,0.486 !建立一个圆面SMRTSIZE,5SECWRITE,CONCRETE,SECT,1SECTYPE,2,BEAM,MESH,SECOFFSET,CENT,CONCRETE,MESHASEL,ALLADELE,ALL,1K,1,-9K,2,-9,-0.2K,3,-7,-0.5K,4,-6.5,-1.8K,5,-5,-1.8K,6,-4.5,-0.5K,7,-4,-0.25K,8,4,-0.25K,9,4.5,-0.25K,10,5,-1.8K,11,6.5,-1.8K,12,7,-0
3、.5K,13,9,-0.2K,14,9A,1,2,3,4,5,6,7,8,9,10,11,12,13,14SECWRITE,ZL,SECT,1SECTYPE,3,BEAM,MESH,ZLET,1,BEAM44 !钢管材料特性MP,EX,1,2.1E11*2 !弹性模量,单位:PaMP,DENS,1,7800 !密度,单位:kg/立方米MP,PRXY,1,0.3 !泊松比ET,2,BEAM44 !拱轴混凝土材料特性MP,EX,2,3.45E10*2MP,DENS,2,2600MP,PRXY,2,0.1667ET,3,BEAM44 !主梁混凝土材料特性MP,EX,3,3.45E10MP,DENS,
4、3,2600MP,PRXY,3,0.1667ET,4,LINK10 !吊杆材料特性(平行钢丝束)MP,EX,4,1.95E11*2MP,DENS,4,7800MP,PRXY,4,0.3KEYOPT,4,3,0 !选择参数设定,只受拉吊杆!*建立有限元模型*定义拱轴线*DIM,x,ARRAY,179 !将顺桥向距离定义为数组*DIM,y,ARRAY,179 !将拱轴高度定义为数组f1=16 !矢高L=90 !跨径*DO,i,1,179,1x(i)=0.5*iy(i)=64/90*x(i)-64/90/90*x(i)*x(i)*enddo*建立主拱圈*DIM,GZ_NUM,ARRAY,181 !拱
5、轴节点号数组GZ_NUM(1)=1001*DO,i,2,181,1GZ_NUM(i)=GZ_NUM(i-1)+1*ENDDON,GZ_NUM(1),0,0,0 !建立拱轴节点N,GZ_NUM(181),90,0,0*DO,i,2,180,1N,GZ_NUM(i),x(i-1),y(i-1),0N,90000,0,16,0 !拱轴参考节点TYPE,1 !单元类型编号MAT,1 !材料类型编号R,1SECNUM,1 !截面编号*DO,i,1,180,1 !生成拱轴单元(钢管)E,GZ_NUM(i),GZ_NUM(i+1),90000CM,ZG_STEEL,ELEMESEL,NONETYPE,2MA
6、T,2R,2SECNUM,2生成拱轴单元(混凝土)CM,ZG_CONCRETE,ELEM*建立主梁*DIM,ZL_NUM,ARRAY,179 !主梁节点号数组ZL_NUM(1)=2002*DO,i,2,179,1ZL_NUM(i)=ZL_NUM(i-1)+1N,ZL_NUM(i),x(i),0,0TYPE,3MAT,3R,3 !定义实常数特性REAL,3 !赋予实常数特性SECNUM,3E,1001,2002,90000 !生成主梁单元*DO,i,1,178,1E,ZL_NUM(i),ZL_NUM(i+1),90000 !注意一定要添加这个截面方位参考点,要不然截面方向不对E,2180,118
7、1,90000CM,ZL,ELEM*定义吊杆参数*(每根吊杆采用54根7mm钢丝)DG_AREA=54*3.1415926*0.25*0.007*0.007 !吊杆面积TYPE,4mat,4R,4,DG_AREA !定义吊杆实常数REAL,4*DO,i,1,17,1E,1000+i*10+1,2000+i*10+1CM,DG,ELEMALLSEL*施加约束*D,1001,UXD,1001,UYD,1001,UZD,1001,ROTXD,1181,UYD,1181,UZ/ESHAPE,1$EPLOT !显示几何模型 运行以上程序,有: 主梁单元编号:361-540 钢管: 1-180 混凝土:1
8、81-360 吊杆单元编号:541-557 对桥梁桥面各节点逐点加载 在181个节点进行单位荷载加载,共181个荷载步/SOLUANTYPE,0 !定义分析类型:静力分析ACEL,0,0,0TIME,1FDELE,ALL,ALLF,1001,FY,-1000SOLVEK=2*DO,I,2,180,1 TIME,K FDELE,ALL,ALL F,2000+I,FY,-1000 SOLVE K=K+1TIME,181F,1181,FY,-1000*DIM,ZL_MZ,ARRAY,181,181!180个单元;181个荷载步*DIM,GZS_MZ,ARRAY,181,181*DIM,GZC_MZ,
9、ARRAY,181,181*DIM,GZ_MZ,ARRAY,181,181*DIM,ZL_QYI,ARRAY,180,181*DIM,ZL_QYJ,ARRAY,180,181*DIM,GZS_QYI,ARRAY,180,181*DIM,GZS_QYJ,ARRAY,180,181*DIM,GZC_QYI,ARRAY,180,181*DIM,GZC_QYJ,ARRAY,180,181*DIM,GZ_QYI,ARRAY,180,181*DIM,GZ_QYJ,ARRAY,180,181*DIM,ZLMZ_MAX,ARRAY,181!内力存储矩阵*DIM,ZLMZ_MIN,ARRAY,181*DIM,G
10、ZMZ_MAX,ARRAY,181*DIM,GZMZ_MIN,ARRAY,181*DIM,ZLQYI_MAX,ARRAY,180*DIM,ZLQYI_MIN,ARRAY,180*DIM,GZQYI_MAX,ARRAY,180*DIM,GZQYI_MIN,ARRAY,180*DIM,ZLQYJ_MAX,ARRAY,180*DIM,ZLQYJ_MIN,ARRAY,180*DIM,GZQYJ_MAX,ARRAY,180*DIM,GZQYJ_MIN,ARRAY,180*获取弯矩、剪力影响矩阵*DO,I,1,181,1 /POST1 SET,I!读取第I荷载步的数据 *DO,J,1,180,1 *GET
11、, ZL_MZ(J,I),ELEM,360+J,SMISC,5!荷载值按I,J端分别存储 *GET, GZS_MZ(J,I),ELEM, J,SMISC,5 *GET, GZC_MZ(J,I),ELEM,180+J,SMISC,5 *GET, ZL_QYI(J,I),ELEM,360+J,SMISC,3 *GET, ZL_QYJ(J,I),ELEM,360+J,SMISC,9 *GET,GZS_QYI(J,I),ELEM, J,SMISC,3 *GET,GZS_QYJ(J,I),ELEM, J,SMISC,9 *GET,GZC_QYI(J,I),ELEM,180+J,SMISC,3 *GET,
12、GZC_QYJ(J,I),ELEM,180+J,SMISC,9 GZ_MZ(J,I)=GZS_MZ(J,I)+GZC_MZ(J,I) GZ_QYI(J,I)=GZS_QYI(J,I)+GZC_QYI(J,I) GZ_QYJ(J,I)=GZS_QYJ(J,I)+GZC_QYJ(J,I) *ENDDO *GET, ZL_MZ(181,I),ELEM,540,SMISC,11 *GET, GZS_MZ(181,I),ELEM,180,SMISC,11 *GET, GZC_MZ(181,I),ELEM,360,SMISC,11 GZ_MZ(181,I)=GZS_MZ(181,I)+GZC_MZ(181,I)*公路I级汽车活载影响线加载效应计算 *车道布置:横向四车道,横向折减系数为0.67;跨径小于150m,不需要考虑纵向折减 *荷载标准值:PK=360kN,qK=10.5kN/m *不考虑冲击系数 NN=4PK=360QK=10.5*弯矩包络图
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1