登录  
 加关注
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

为梦飞翔

我希望我身边的亲人和朋友天天快乐!

 
 
 

日志

 
 

钢管混凝土拱桥  

2010-04-27 19:46:34|  分类: Ansys |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

/prep7
/title,the analyse of steel arch bridge
!**********************************************************
et,1,82                     !建立辅助单元作出截面特性
cyl4,0.6,0.6,0.6,,0.586     !以下命令是用于建立一个空心钢管截面
cyl4,0.6,2.4,0.6,,0.586
rectng,0.35,0.364,0.8,2.2
rectng,0.85,0.836,0.8,2.2
asel,all
aptn,all
adele,5,8,1                  !删除多余的面积
allsel
aadd,all
smrtsize,5
amesh,all
secwrite,gg,sect,,1
sectype,1,beam,mesh          !截面代号1
secoffset,cent,,,
secread,'gg','sect',,mesh    !读入空心钢管截面
asel,all
aclear,all
adele,all,,,1                !将图形及面、线、点全部删除(只保留截面文件)


/replot
cyl4,0.6,0.6,0.586           !以下命令是用于建立一个空心钢管内的混凝土截面
cyl4,0.6,2.4,0.586
rectng,0.364,0.836,0.8,2.2
allsel
aadd,all
smrtsize,5
amesh,all
secwrite,hnt,sect,,1
sectype,2,beam,mesh           !截面代号2
secoffset,cent,,,
secread,'hnt','sect',,mesh    !读入混凝土(空心钢管内的)截面
asel,all
aclear,all
adele,all,,,1                 !将图形及面、线、点全部删除(只保留截面文件)

/replot

k,1,-0.5,0                    !以下命令是用于建立一个横梁标准截面
k,2,-0.75,0.25
k,3,-0.75,0.5
k,4,-0.55,0.5
k,5,-0.55,0.85
k,6,0.55,0.85
k,7,0.55,0.5
k,8,0.75,0.5
k,9,0.75,0.25
k,10,0.5,0
a,1,2,3,4,5,6,7,8,9,10        !生成横梁标准截面1 A-A
smrtsize,5
amesh,all
secwrite,hl1,sect,,1
sectype,4,beam,mesh           !截面代号4
secoffset,cent,,,
secread,'hl1','sect',,mesh    !读入横梁标准截面1 A-A
asel,all
aclear,all
adele,all,,,1

/replot

k,1,-0.5,0                     !以下命令是用于建立一个横梁标准截面b-b
k,2,-0.5,1
k,3,-0.75,1.25
k,4,-0.75,1.5
k,5,0.75,1.5
k,6,0.75,1.25
k,7,0.5,1
k,8,0.5
a,1,2,3,4,5,6,7,8                !生成横梁标准截面2 b-b
smrtsize,5
amesh,all
secwrite,hl2,sect,,1             !读入横梁标准截面2
sectype,5,beam,mesh              !截面代号5
secoffset,cent,,,
secread,'hl2','sect',,mesh
asel,all
aclear,all
adele,all,,,1


/replot

k,1,-0.5,0                        !以下命令是用于建立一个横梁标准截面c-c
k,2,-0.5,1.04
k,3,-0.75,1.04+0.25
k,4,-0.75,1.04+0.5
k,5,-0.55,1.04+0.5
k,6,-0.55,1.04+0.5+0.35
k,7,0.55,1.04+0.5+0.35
k,8,0.55,1.04+0.5
k,9,0.75,1.04+0.5
k,10,0.75,1.04+0.25
k,11,0.5,1.04
k,12,0.5,0
k,13,-0.2,0.58
k,14,-0.2,1.04+0.5-0.3
k,15,0.2,1.04+0.5-0.3
k,16,0.2,0.58
a,1,2,3,4,5,6,7,8,9,10,11,12
a,13,14,15,16
asba,1,2
smrtsize,5
amesh,all
secwrite,hl3,sect,,1
sectype,6,beam,mesh                    ! 截面代号6
secoffset,cent,,,
secread,'hl3','sect',,mesh             !横梁标准截面3
asel,all
aclear,all
adele,all,,,1

/replot

ldele,all
kdele,all


k,1,-0.5,0                              !以下命令是用于建立一个横梁标准截面d-d
k,2,-0.5,1.23
k,3,-0.75,1.23+0.25
k,4,-0.75,1.23+0.5
k,5,-0.55,1.23+0.5
k,6,-0.55,1.23+0.5+0.35
k,7,0.55,1.23+0.5+0.35
k,8,0.55,1.23+0.5
k,9,0.75,1.23+0.5
k,10,0.75,1.23+0.25
k,11,0.5,1.23
k,12,0.5,0
k,13,-0.2,0.38
k,14,-0.3,0.38+0.2
k,15,-0.3,1.23+0.5-0.1-0.2
k,16,-0.2,1.23+0.5-0.1
k,17,0.2,1.23+0.5-0.1
k,18,0.3,1.23+0.5-0.1-0.2
k,19,0.3,0.38+0.2
k,20,0.2,0.38
a,1,2,3,4,5,6,7,8,9,10,11,12
a,13,14,15,16,17,18,19,20
asba,1,2                                     !横梁标准截面4 D-D
smrtsize,5
amesh,all
secwrite,hl4,sect,,1
sectype,7,beam,mesh                          !截面代号7
secoffset,cent,,,
secread,'hl4','sect',,mesh                   !读入横梁标准截面4
asel,all
aclear,all
adele,all,,,1

 

/replot

k,1,-0.5,0                                   !以下命令是用于建立一个横梁标准截面e-e
k,2,-0.5,1.23
k,3,-0.75,1.23+0.25
k,4,-0.75,1.23+0.5
k,5,-0.55,1.23+0.5
k,6,-0.55,1.23+0.5+0.35
k,7,0.55,1.23+0.5+0.35
k,8,0.55,1.23+0.5
k,9,0.75,1.23+0.5
k,10,0.75,1.23+0.25
k,11,0.5,1.23
k,12,0.5,0
a,1,2,3,4,5,6,7,8,9,10,11,12                  !横梁标准截面5 E-E
smrtsize,5
amesh,all
secwrite,hl5,sect,,1
sectype,8,beam,mesh                           ! 截面代号8
secoffset,cent,,,
secread,'hl5','sect',,mesh                    !横梁标准截面5
asel,all
aclear,all
adele,all,,,1

/replot
etdele,1              !将辅助单元类型删除
/replot
!*********************************************************************
et,1,beam44          !钢管特性
mp,ex,1,2.1e11
mp,dens,1,7800
mp,prxy,1,0.3
n,90000,0,0,30       !参考点

et,2,beam44         !拱肋内混凝土特性(50号混凝土)
mp,ex,2,3.5e10
mp,dens,2,2600
mp,prxy,2,0.1667

et,3,beam44         !纵梁特性(30号混凝土)
mp,ex,3,3e10
mp,dens,3,2600
mp,prxy,3,0.1667

et,4,beam44         !横梁特性(30号混凝土)
mp,ex,4,3e10
mp,dens,4,2600
mp,prxy,4,0.1667

et,5,beam44          !风撑特性
mp,ex,5,2.1e11
mp,dens,5,7800
mp,prxy,5,0.3

et,6,link10          !吊杆特性(钢绞线)
mp,ex,6,1.9e11
mp,dens,6,7800
mp,prxy,6,0.3
keyopt,6,3,0        !只拉吊杆

et,7,beam44         !盖梁特性(30号混凝土)
mp,ex,7,3e10
mp,dens,7,2600
mp,prxy,7,0.1667

et,8,beam44         !墩柱特性(40号混凝土)
mp,ex,8,3.3e10
mp,dens,8,2600
mp,prxy,8,0.1667

et,9,beam44         !横系梁特性(30号混凝土)
mp,ex,9,3e10
mp,dens,9,2600
mp,prxy,9,0.1667

et,10,link10          !系杆特性(钢绞线)
mp,ex,10,2e11
mp,dens,10,7800
mp,prxy,10,0.3
keyopt,10,3,0        !只拉吊杆

et,11,beam44         !承台桩基础特性(25号混凝土)
mp,ex,11,2.85e10
mp,dens,11,2500
mp,prxy,11,0.1667

n,90000,0,0,30
!*****************************************************************
!开始建模
n,1,0,2.025      !0号横梁(端横梁)
n,2,6.25,2.025
n,3,6.25+5.1,2.025
n,4,6.25+5.1+0.2+1.9/2,2.025        !12.5米位置
n,5,6.25+5.1+2.3,2.025
n,6,6.25+5.1+2.3+2.23,2.025
nsym,x,10,2,6,1

type,4                !横梁
r,1
real,1
mat,4
secnum,7
                   !单元1到10为0号横梁
e,1,2      !截面7 D-D

secnum,6
e,2,3      !截面6 D-D

secnum,5
e,3,4      !截面5 D-D
e,4,5      !截面5 D-D

secnum,4
e,5,6      !截面4 a-a

secnum,7
e,1,12     !截面7 D-D

secnum,6
e,12,13    !截面6 D-D

secnum,5
e,13,14    !截面5 D-D
e,14,15    !截面5 D-D

secnum,4
e,15,16    !截面4 a-a

!1号端横梁
egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2

!其余横梁(2到17号)
egen,17,100,11,20,1,,,,,,0,5.3+1.5

!18号横梁(端横梁)
egen,2,100,171,180,1,,,,,,0,4.7+1.5/2+1.25/2

cm,hl,elem
esel,none

!******************************************************
!纵梁单元
type,3
r,2
real,2
mat,3
sectype,10,beam,rect,zl            !纵梁截面代号10
secdata,0.4,1.2
secnum,10
                    !0号纵梁(端横梁)
e,1,101
e,2,102
e,3,103
e,4,104
e,5,105
e,12,112
e,13,113
e,14,114
e,15,115

                    !中间的纵梁(1到16号)
*do,i,0,15,1
 e,101+i*100,201+i*100
 e,102+i*100,202+i*100
 e,104+i*100,204+i*100
 e,112+i*100,212+i*100
 e,114+i*100,214+i*100
*enddo

!17号纵梁(端横梁)
e,1701,1801
e,1702,1802
e,1703,1803
e,1705,1805

e,1712,1812
e,1713,1813
e,1715,1815

cm,zl,elem
esel,none

!定义拱轴线(参考姚玲森《桥梁工程》p299页相关公式)
*dim,y,array,16       !纵桥向距离
*dim,cs,array,16      !定义分值1-2x/L
*dim,z,array,16       !拱轴高度

!吊杆上的点
                              !吊杆的y向坐标(顺桥向)

y(1)=8.1,14.9,21.7,28.5,35.3,42.1,48.9,55.7,62.5

!风撑上的点
                              !风撑的y向坐标(顺桥向)
y(11)=21.7,28.5,42.1,48.9,59.1

                              !0号横梁的y向的坐标(顺桥向)
y(16)=2.025
f1=25        !矢高
m1=1.1       !悬链线拱轴系数
k=log(m1+sqrt(m1*m1-1))     !k=ch-1m

*do,i,1,16,1
  cs(i)=1-y(i)*2/125
  z(i)=25-f1/(m1-1)*(cosh(k*cs(i))-1)    !拱轴高度
*enddo

*dim,ygd,array,23       !定义预拱度(是在后续工作中反复调算的结果),初始运行时可将预拱度全部设置为0
ygd(1)=0.40004E-02+1.29E-02
ygd(2)=0.10685E-01+2.22E-02
ygd(3)=0.19483E-01+3.03E-02
ygd(4)=0.29219E-01+3.69E-02
ygd(5)=0.39022E-01+4.23E-02
ygd(6)=0.48243E-01+4.65E-02
ygd(7)=0.56390E-01+4.95E-02
ygd(8)=0.63070E-01+5.15E-02
ygd(9)=0.68014E-01+5.21E-02
ygd(10)=0.71003E-01+5.15E-02
ygd(11)=0.71739E-01+4.95E-02
ygd(12)=0.69962E-01+4.65E-02
ygd(13)=0.65514E-01+4.23E-02
ygd(14)=0.58289E-01+3.69E-02
ygd(15)=0.48248E-01+3.04E-02
ygd(16)=0.35507E-01+2.24E-02
ygd(17)=0.20415E-01+1.23E-02
ygd(18)=0.65770E-01 +0.40234E-01
ygd(19)=0.69769E-01 +0.41215E-01


!生成各个吊杆的z向节点(竖直方向)
ngen,2,10000,104,114,10,0,0,z(1)+ygd(1)     !1到17号吊杆节点
ngen,2,10000,204,214,10,0,0,z(2)+ygd(2)
ngen,2,10000,304,314,10,0,0,z(3)+ygd(3)
ngen,2,10000,404,414,10,0,0,z(4)+ygd(4)
ngen,2,10000,504,514,10,0,0,z(5)+ygd(5)
ngen,2,10000,604,614,10,0,0,z(6)+ygd(6)
ngen,2,10000,704,714,10,0,0,z(7)+ygd(7)
ngen,2,10000,804,814,10,0,0,z(8)+ygd(8)
ngen,2,10000,904,914,10,0,0,z(9)+ygd(9)
ngen,2,10000,1004,1014,10,0,0,z(8)+ygd(10)
ngen,2,10000,1104,1114,10,0,0,z(7)+ygd(11)
ngen,2,10000,1204,1214,10,0,0,z(6)+ygd(12)
ngen,2,10000,1304,1314,10,0,0,z(5)+ygd(13)
ngen,2,10000,1404,1414,10,0,0,z(4)+ygd(14)
ngen,2,10000,1504,1514,10,0,0,z(3)+ygd(15)
ngen,2,10000,1604,1614,10,0,0,z(2)+ygd(16)
ngen,2,10000,1704,1714,10,0,0,z(1)+ygd(17)

!重复定义吊杆节点是因为以后吊杆转动要进行放松

ngen,2,3000,104,114,10,0,0,0        !1到17号吊杆节点
ngen,2,3000,204,214,10,0,0,0
ngen,2,3000,304,314,10,0,0,0
ngen,2,3000,404,414,10,0,0,0
ngen,2,3000,504,514,10,0,0,0
ngen,2,3000,604,614,10,0,0,0
ngen,2,3000,704,714,10,0,0,0
ngen,2,3000,804,814,10,0,0,0
ngen,2,3000,904,914,10,0,0,0
ngen,2,3000,1004,1014,10,0,0,0
ngen,2,3000,1104,1114,10,0,0,0
ngen,2,3000,1204,1214,10,0,0,0
ngen,2,3000,1304,1314,10,0,0,0
ngen,2,3000,1404,1414,10,0,0,0
ngen,2,3000,1504,1514,10,0,0,0
ngen,2,3000,1604,1614,10,0,0,0
ngen,2,3000,1704,1714,10,0,0,0

ngen,2,6000,104,114,10,0,0,z(1)+ygd(1)     !1到17号吊杆节点
ngen,2,6000,204,214,10,0,0,z(2)+ygd(2)
ngen,2,6000,304,314,10,0,0,z(3)+ygd(3)
ngen,2,6000,404,414,10,0,0,z(4)+ygd(4)
ngen,2,6000,504,514,10,0,0,z(5)+ygd(5)
ngen,2,6000,604,614,10,0,0,z(6)+ygd(6)
ngen,2,6000,704,714,10,0,0,z(7)+ygd(7)
ngen,2,6000,804,814,10,0,0,z(8)+ygd(8)
ngen,2,6000,904,914,10,0,0,z(9)+ygd(9)
ngen,2,6000,1004,1014,10,0,0,z(8)+ygd(10)
ngen,2,6000,1104,1114,10,0,0,z(7)+ygd(11)
ngen,2,6000,1204,1214,10,0,0,z(6)+ygd(12)
ngen,2,6000,1304,1314,10,0,0,z(5)+ygd(13)
ngen,2,6000,1404,1414,10,0,0,z(4)+ygd(14)
ngen,2,6000,1504,1514,10,0,0,z(3)+ygd(15)
ngen,2,6000,1604,1614,10,0,0,z(2)+ygd(16)
ngen,2,6000,1704,1714,10,0,0,z(1)+ygd(17)

!定义吊杆的参数(每根吊杆都采用109根直径为7mm的高强钢丝)
dgarea=109*3.1415926*0.25*0.007*0.007    !吊杆面积
dgyl=0.3*1670*1e6                        !吊杆初应力取用0.3倍的抗拉强度(N/m2)
dgl=dgyl*dgarea                          !吊杆力
edg=1.95e11

*dim,L0,array,17                        !吊杆初应力为0时的长度
*dim,L,array,17                         !吊杆后的长度
*dim,deltL,array,17                     !吊杆受力后的伸长量
*dim,istrn,array,17                     !吊杆应变
*dim,fdg,array,17                       !吊杆上施加力

*do,i,1,9,1
  L0(i)=z(i)          !应力为0时的单元长度
  L0(18-i)=z(i)       !沿中轴对称的单元
*enddo
fdg(1)=dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl,dgl     !1到17号吊杆施加的预应力

type,6
mat,6
*do,i,1,17,1
deltL(i)=fdg(i)*L0(i)/(edg*dgarea)   !吊杆受力后的增量
L(i)=L0(i)+deltL(i)
istrn(i)=deltL(i)/L(i)               !吊杆受力后的应变
r,i+1000,dgarea,istrn(i)             !定义吊杆的实常数
*enddo

*do,i,1,17,1                  !生成吊杆单元
real,i+1000
e,i*100+4+3000,i*100+4+6000
e,i*100+10+4+3000,i*100+10+4+6000
*enddo

cm,dg,elem
esel,none
allsel

!***************************************
r,1018
real,1018     !定义一个全为0的实常数便于利用截面实常数时不产生冲突
!***********************************************************************
!生成风撑上的节点
                            !风撑上节点坐标
n,20001,0,21.7,z(11)
n,20002,0,28.5,z(12)
n,20003,0,42.1,z(13)
n,20004,0,48.9,z(14)

n,20005,0,59.1,z(15)
n,20006,0,65.9,z(15)
n,20007,0,76.1,z(14)
n,20008,0,82.9,z(13)

n,20009,0,96.5,z(12)
n,20010,0,103.3,z(11)
n,20011,12.5,21.7,z(11)
n,20012,12.5,28.5,z(12)

n,20013,12.5,42.1,z(13)
n,20014,12.5,48.9,z(14)
n,20015,12.5,59.1,z(15)+ygd(18)
n,20016,12.5,65.9,z(15)+ygd(19)

n,20017,12.5,76.1,z(14)
n,20018,12.5,82.9,z(13)
n,20019,12.5,96.5,z(12)
n,20020,12.5,103.3,z(11)

n,20021,-12.5,21.7,z(11)
n,20022,-12.5,28.5,z(12)
n,20023,-12.5,42.1,z(13)
n,20024,-12.5,48.9,z(14)

n,20025,-12.5,59.1,z(15)+ygd(18)
n,20026,-12.5,65.9,z(15)+ygd(18)
n,20027,-12.5,76.1,z(14)
n,20028,-12.5,82.9,z(13)

n,20029,-12.5,96.5,z(12)
n,20030,-12.5,103.3,z(11)
n,20000,0,62.5,z(15)              !中撑中点

!风撑1(倒K撑)
type,5
mat,5
sectype,500,beam,ctube,windc          !定义中撑的形状为空心圆
secdata,0.486,0.50                    !赋予截面形状参数,直径1m,厚度1.4cm
secnum,500                            !调用截面

!生成风撑单元
                        !风撑1(倒K撑)
e,20002,10404,90000
e,20002,10414,90000
e,20002,10304,90000
e,20002,10314,90000
                        !风撑2(倒K撑)
e,20004,10704,90000
e,20004,10714,90000
e,20004,10604,90000
e,20004,10614,90000
                        !中撑3(中撑X撑)
e,20000,20015,90000
e,20000,20016,90000
e,20000,20025,90000
e,20000,20026,90000
                        !风撑4(K撑)
e,20007,11104,90000
e,20007,11114,90000
e,20007,11204,90000
e,20007,11214,90000
                        !风撑5(K撑)
e,20009,11404,90000
e,20009,11414,90000
e,20009,11504,90000
e,20009,11514,90000

cm,fc,elem
esel,none

!******************************************
!连接拱轴线
!定义拱轴上的坐标点
n,30000,12.5,0,0                    !左拱轴起点
n,30001,12.5,125,0                  !左拱轴端点
n,30002,12.5,2.205,z(16)            !左横梁起点
n,30003,12.5,125-2.205,z(16)        !左横梁端点

n,40000,-12.5,0,0                    !右拱轴起点
n,40001,-12.5,125,0                  !右拱轴端点
n,40002,-12.5,2.205,z(16)            !右横梁起点
n,40003,-12.5,125-2.205,z(16)        !右横梁端点

*dim,zgzx,array,23             !定义拱轴线,左轴线
*dim,ygzx,array,23             !定义拱轴线,右轴线
                      !拱轴线上的节点号
zgzx(1)=30000,30002,10104,10204,10304,10404,10504,10604,10704,10804,20015,10904,20016,11004,11104
zgzx(16)=11204,11304,11404,11504,11604,11704,30003,30001

ygzx(1)=40000,40002,10114,10214,10314,10414,10514,10614,10714,10814,20025,10914,20026,11014,11114
ygzx(16)=11214,11314,11414,11514,11614,11714,40003,40001

n,90001,12.5,0,30       !钢管左轴线参考点
n,90002,-12.5,0,30      !钢管右轴线参考点

type,1
mat,1
r,2000
secnum,1
*do,i,1,22,1                    !生成拱轴单元(钢管)
 e,zgzx(i),zgzx(i+1),90001
 e,ygzx(i),ygzx(i+1),90002
*enddo

cm,gzgg,elem
esel,none

type,2
mat,2
r,2001
secnum,2
*do,i,1,22,1                    !生成拱轴单元(混凝土)
 e,zgzx(i),zgzx(i+1),90001
 e,ygzx(i),ygzx(i+1),90002
*enddo

cm,gzhnt,elem
esel,none

!******************************
!建立拱轴底座
type,9
r,6
real,6
mat,9
sectype,70,beam,rect,dizuo   !定义截面形状
secdata,1.4,4                !赋予截面参数
secnum,70
e,30001,30003,90001
e,40001,40003,90002
e,30000,30002,90001
e,40000,40002,90002

cm,gzdz,elem
esel,none

!******************************
!建立盖梁
ngen,2,2000,1,16,1,0,-2.025,0      !复制节点到拱轴线起点位置
ngen,2,100,1801,1816,1,0,2.025,0   !复制节点到拱轴线终点位置

type,7
mat,7
r,7
real,7
sectype,60,beam,rect,gail
secdata,4,1.6

secnum,60
n,90005,0,0,20
n,90006,0,125,20
                       !盖梁单元125m处
e,1901,1902,90006
e,1902,1903,90006
e,1903,1904,90006
e,1904,1905,90006
e,1905,1906,90006

e,1901,1912,90006
e,1912,1913,90006
e,1913,1914,90006
e,1914,1915,90006
e,1915,1916,90006

e,2001,2002,90006
e,2002,2003,90006
e,2003,2004,90006
e,2004,2005,90006
e,2005,2006,90006

e,2001,2012,90006
e,2012,2013,90006
e,2013,2014,90006
e,2014,2015,90006
e,2015,2016,90006

cm,gailiang,elem
esel,none

!******************************
!建立墩柱
n,90003,14.5,0,0
n,90004,-14.5,0,0
type,8
mat,8
r,8
real,8
sectype,50,beam,hrec,dun
secdata,3,4,0.8,0.8,0.8,0.8
ngen,2,10000,2004,2014,10,0,0,-10   !墩节点(距离桥面-10m)
ngen,2,10000,1904,1914,10,0,0,-10

secnum,50
e,2004,12004,90003
e,2014,12014,90004
e,1904,11904,90003
e,1914,11914,90004

cm,dunzhu,elem
esel,none

!******************************
!承台单元
type,11
mat,11
r,1
sectype,40,beam,rect,chengtai
secdata,12.7,9
r,11
real,11
secnum,40
ngen,2,20000,2004,2014,10,0,0,-11.5
ngen,2,20000,1904,1914,10,0,0,-11.5
ngen,2,30000,2004,2014,10,0,0,-13
ngen,2,30000,1904,1914,10,0,0,-13
                                   !生成承台单元
e,12004,22004,90003
e,11904,21904,90003
e,12014,22014,90004
e,11914,21914,90004

e,22004,32004,90003
e,21904,31904,90003
e,22014,32014,90004
e,21914,31914,90004

cm,chtai,elem
esel,none

!******************************
!横系梁
type,9
mat,9
r,9
real,9
n,90007,12.5,0,10
n,90008,-12.5,0,10
sectype,90,beam,rect,hxl
secdata,3,3
secnum,90
e,22004,22014,90007
e,21904,21914,90008

cm,hxl,elem
esel,none


!******************************
!桩
                      !定义桩节点
n,50000,6.15+1.85,2.7,-13
n,50001,6.15+1.85,-2.7,-13
n,50002,6.15+1.85+4.5,2.7,-13
n,50003,6.15+1.85+4.5,-2.7,-13
n,50004,6.15+12.7-1.85,2.7,-13
n,50005,6.15+12.7-1.85,-2.7,-13

nsel,s,,,50000,50005,1
nsym,x,10,50000,50005,1
allsel
ngen,2,100,50000,50015,1,0,125,0     !首先生成平面内的24个节点

!每3m生成一层节点(共8层)
ngen,10,1000,50000,51000,1,0,0,-2.5

!生成桩单元

esel,none
type,11
mat,11
r,11
real,11
sectype,20,beam,csolid,pier
secdata,0.9
secnum,20
*do,i,0,8,1
  *do,k,0,5,1
    e,50000+k+i*1000,50000+k+(i+1)*1000,90003
    e,50000+k+10+i*1000,50000+k+10+(i+1)*1000,90003
    e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003
    e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003
  *enddo
*enddo

cm,zhuang,elem
esel,none

!*************************************
!在桩上施加弹簧单元
!第一层距离承台-3m
!计算弹簧的弹性系数
d=1.8
Kf=0.9
k0=1+1/d
b1=0.6                 !和每排桩有关的系数,当为2排时取0.6
L1=3.6                 !桩间净距
h1=3*(d+1)             !桩在地面或冲刷线下的计算深度
k=b1+(1-b1)/0.6*L1/h1
B=Kf*k0*k*d            !换算宽度

!第1层土 粉质量粘土 m=8 厚度5
!第2层土 粘土(细砂) m=15 厚度5
!第3层土 中砂 m=20 厚度5
!第4层土 粘土 m=25 厚度5
!第5层   基岩 c0=300 厚度10m

*dim,m,array,5
m(1)=8,15,20,25,4000             !其中4000为基岩系数

*dim,h,array,5                   !土的厚度
h(1)=5,5,5,5,10

*dim,c,array,4
c(1)=m(1)*h(1),m(2)*(h(1)+h(2)),m(3)*(h(1)+h(2)+h(3)),m(4)*(h(1)+h(2)+h(3)+h(4))

*dim,e,array,4                 !弹簧弹性系数
e(1)=0.5*c(1)*h(1)*B*1e6
*do,i,2,4,1
  e(i)=0.5*(c(i-1)+c(i))*h(i)*B*1e6
*enddo

                     !弹簧单元
et,20,combin14
r,20,e(1)                            !第1层土弹簧弹性系数
r,30,e(2)                            !第2层土弹簧弹性系数
r,40,e(3)                            !第3层土弹簧弹性系数
r,50,e(4)                            !第4层土弹簧弹性系数
r,60,m(5)*3.1415926*B*B/4*1e6        !基础岩石弹性系数

!第1层底部距离承台-5m
ngen,2,10000,51000,51115,1,1,0,0     !x方向
ngen,2,20000,51000,51115,1,0,1,0     !y方向

type,20
real,20
*do,i,0,5,1
 e,51000+i,51000+i+10000
 e,51000+i+10,51000+i+10+10000
 e,51100+i,51100+i+10000
 e,51100+i+10,51100+i+10+10000
 e,51000+i,51000+i+20000
 e,51000+i+10,51000+i+10+20000
 e,51100+i,51100+i+20000
 e,51100+i+10,51100+i+10+20000
*enddo

cm,spring1,elem
esel,none


!第2层底部距离承台-10m
ngen,2,10000,53000,53115,1,1,0,0     !x方向
ngen,2,20000,53000,53115,1,0,1,0     !y方向

type,20
real,30
*do,i,0,5,1
 e,53000+i,53000+i+10000
 e,53000+i+10,53000+i+10+10000
 e,53100+i,53100+i+10000
 e,53100+i+10,53100+i+10+10000
 e,53000+i,53000+i+20000
 e,53000+i+10,53000+i+10+20000
 e,53100+i,53100+i+20000
 e,53100+i+10,53100+i+10+20000
*enddo

cm,spring2,elem
esel,none


!第3层底部距离承台-15m
ngen,2,10000,55000,55115,1,1,0,0     !x方向
ngen,2,20000,55000,55115,1,0,1,0     !y方向

type,20
real,40
*do,i,0,5,1
 e,55000+i,55000+i+10000
 e,55000+i+10,55000+i+10+10000
 e,55100+i,55100+i+10000
 e,55100+i+10,55100+i+10+10000
 e,55000+i,55000+i+20000
 e,55000+i+10,55000+i+10+20000
 e,55100+i,55100+i+20000
 e,55100+i+10,55100+i+10+20000
*enddo

cm,spring3,elem
esel,none


!第4层底部距离承台-20m
ngen,2,10000,57000,57115,1,1,0,0     !x方向
ngen,2,20000,57000,57115,1,0,1,0     !y方向

type,20
real,50
*do,i,0,5,1
 e,57000+i,57000+i+10000
 e,57000+i+10,57000+i+10+10000
 e,57100+i,57100+i+10000
 e,57100+i+10,57100+i+10+10000
 e,57000+i,57000+i+20000
 e,57000+i+10,57000+i+10+20000
 e,57100+i,57100+i+20000
 e,57100+i+10,57100+i+10+20000
*enddo

cm,spring4,elem
esel,none


!第5层底部距离承台-22.5m
ngen,2,10000,58000,58115,1,1,0,0     !x方向
ngen,2,20000,58000,58115,1,0,1,0     !y方向

type,20
real,60
*do,i,0,5,1
 e,58000+i,58000+i+10000
 e,58000+i+10,58000+i+10+10000
 e,58100+i,58100+i+10000
 e,58100+i+10,58100+i+10+10000
 e,58000+i,58000+i+20000
 e,58000+i+10,58000+i+10+20000
 e,58100+i,58100+i+20000
 e,58100+i+10,58100+i+10+20000
*enddo

cm,spring5,elem
esel,none

!第6层底部距离承台-25m
ngen,2,10000,59000,59115,1,1,0,0     !x方向
ngen,2,20000,59000,59115,1,0,1,0     !y方向

type,20
real,60
*do,i,0,5,1
 e,59000+i,59000+i+10000
 e,59000+i+10,59000+i+10+10000
 e,59100+i,59100+i+10000
 e,59100+i+10,59100+i+10+10000
 e,59000+i,59000+i+20000
 e,59000+i+10,59000+i+10+20000
 e,59100+i,59100+i+20000
 e,59100+i+10,59100+i+10+20000
*enddo

cm,spring6,elem

allsel

!******************************************
!剔除不需要进行对称复制的单元
allsel
nsel,s,loc,z,-50,-0.01
nsel,r,loc,y,-10,10
cm,zdnode,node
esln,s,,all

cm,zdelem,elem
allsel
nsel,u,,,zdnode
nsym,y,100000,all
allsel
esel,u,,,zdelem
esym,,100000,all
allsel

!*****************************
!定义系杆单元(全桥通长)
!系杆单元
!定义系杆参数(每根系杆都采用8束17*7*d5高强钢丝)
xgarea=8*27*7*3.1415926*0.25*0.005*0.005      !系杆集束面积
exg=1.95e11
*dim,fxg,array,2
*dim,Llast,array,2
*dim,Listrn,array,2
*dim,detL,array,2
xgyl=0.3*1860*1e6
xgl=xgyl*xgarea
fxg(1)=xgl
fxg(2)=xgl
Lzero=250
*do,i,1,2,1
  detL(i)=Lzero*fxg(i)/(exg*xgarea)
  Llast(i)=Lzero+detL(i)
  Listrn(i)=detl(i)/Llast(i)
 r,i+100,xgarea,Listrn(i)
*enddo

type,10
mat,10
real,101
e,30001,130001            !系杆1(250m)
e,40001,140001            !系杆2(250m)
cm,xigan,elem
esel,none
allsel

!*****************************
!将盖梁和拱轴交界点进行耦合,放松uy,rotx
cp,1,ux,2004,30000
cp,next,uz,2004,30000
cp,next,roty,2004,30000
cp,next,rotz,2004,30000

cp,next,ux,2014,40000
cp,next,uz,2014,40000
cp,next,roty,2014,40000
cp,next,rotz,2014,40000

cp,next,ux,1904,30001
cp,next,uz,1904,30001
cp,next,roty,1904,30001
cp,next,rotz,1904,30001

cp,next,ux,1914,40001
cp,next,uz,1914,40001
cp,next,roty,1914,40001
cp,next,rotz,1914,40001

cp,next,ux,101914,140001
cp,next,uz,101914,140001
cp,next,roty,101914,140001
cp,next,rotz,101914,140001

cp,next,ux,101904,130001
cp,next,uz,101904,130001
cp,next,roty,101904,130001
cp,next,rotz,101904,130001

allsel
nsel,s,,,30000,130000,100000
nsel,s,,,40000,140000,100000

nummrg,node
allsel

!***************************************
*do,i,1,3,1
  cp,next,all,1800+i,1900+i              !盖梁和0号横梁耦合
*enddo
allsel
*do,i,5,6,1
  cp,next,all,1800+i,1900+i
*enddo

allsel
*do,i,1,2,1
  cp,next,all,1800+i+11,1900+i+11
  cp,next,all,i+11,2000+i+11,i+100000+11
*enddo
allsel
*do,i,4,5,1
  cp,next,all,1800+i+11,1900+i+11
  cp,next,all,i+11,2000+i+11,i+100000+11
*enddo

allsel
*do,i,1,3,1
  cp,next,all,i,2000+i,100000+i
*enddo
allsel
*do,i,5,6,1
  cp,next,all,i,2000+i,100000+i
*enddo

allsel
*do,i,1,3,1
  cp,next,all,101800+i,101900+i              !盖梁和0号横梁耦合
*enddo
allsel
*do,i,5,6,1
  cp,next,all,101800+i,101900+i 
*enddo

allsel
*do,i,1,2,1
  cp,next,all,101800+i+11,101900+i+11
*enddo
allsel
*do,i,4,5,1
  cp,next,all,101800+i+11,101900+i+11
*enddo


!**************************************
!桩节点和承台节点耦合
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-3,3
nsel,r,loc,x,8,20
cp,next,all,all

allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-3,3
nsel,r,loc,x,-20,-8
cp,next,all,all

allsel
nsel,s,loc,z,-13
nsel,r,loc,y,122,128
nsel,r,loc,x,-20,-8
cp,next,all,all

allsel
nsel,s,loc,z,-13
nsel,r,loc,y,122,128
nsel,r,loc,x,8,20
cp,next,all,all

allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-128,-122
nsel,r,loc,x,8,20
cp,next,all,all

allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-128,-122
nsel,r,loc,x,-20,-8
cp,next,all,all
allsel

!***********************************
!将吊杆转动进行放松
*do,i,1,17,1
  cp,next,ux,i*100+4,i*100+4+3000
  cp,next,ux,i*100+4+6000,i*100+4+10000
  cp,next,uy,i*100+4,i*100+4+3000
  cp,next,uy,i*100+4+6000,i*100+4+10000
  cp,next,uz,i*100+4,i*100+4+3000
  cp,next,uz,i*100+4+6000,i*100+4+10000

  cp,next,ux,i*100+4+10,i*100+4+3000+10
  cp,next,ux,i*100+4+6000+10,i*100+4+10000+10
  cp,next,uy,i*100+4+10,i*100+4+3000+10
  cp,next,uy,i*100+4+6000+10,i*100+4+10000+10
  cp,next,uz,i*100+4+10,i*100+4+3000+10
  cp,next,uz,i*100+4+6000+10,i*100+4+10000+10

  cp,next,ux,i*100+4+100000,i*100+4+3000+100000
  cp,next,ux,i*100+4+6000+100000,i*100+4+10000+100000
  cp,next,uy,i*100+4+100000,i*100+4+3000+100000
  cp,next,uy,i*100+4+6000+100000,i*100+4+10000+100000
  cp,next,uz,i*100+4+100000,i*100+4+3000+100000
  cp,next,uz,i*100+4+6000+100000,i*100+4+10000+100000

  cp,next,ux,i*100+4+100000+10,i*100+4+3000+100000+10
  cp,next,ux,i*100+4+6000+100000+10,i*100+4+10000+100000+10
  cp,next,uy,i*100+4+100000+10,i*100+4+3000+100000+10
  cp,next,uy,i*100+4+6000+100000+10,i*100+4+10000+100000+10
  cp,next,uz,i*100+4+100000+10,i*100+4+3000+100000+10
  cp,next,uz,i*100+4+6000+100000+10,i*100+4+10000+100000+10
*enddo
allsel

!***********************************
!将土弹簧的一端约束
nsel,s,loc,x,13.5
nsel,a,loc,x,9
nsel,a,loc,x,18
nsel,a,loc,x,-11.5
nsel,a,loc,x,-7
nsel,a,loc,x,-16
nsel,a,loc,y,3.7
nsel,a,loc,y,-1.7
nsel,a,loc,y,128.7
nsel,a,loc,y,123.3
nsel,a,loc,y,-128.7
nsel,a,loc,y,-123.3
nsel,a,loc,z,-40,-13
d,all,all
allsel

nsel,s,loc,z,-36.5
d,all,all

/color,pbak,off
/replot
allsel
nsel,s,loc,y,0
nsel,r,loc,z,0
nummrg,node,all
esel,none
allsel

  评论这张
 
阅读(232)| 评论(0)

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018