资源描述：
《时程分析实例.docx》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、时程分析实例finish/clear/filename,ex6-1!设定工作名/title,TheDynamicAnalysisoftheBridge!设定工作标题/prep7!进入前处理模块!**********************************!定义单元类型及单元实常数!**********************************Et,1,beam4!定义单元类型Et,2,shell63Et,3,link10!定义梁单元实常数*dim,BBS,,43!初始化描述塔上部截面几何特征的数组*dim,HHS,,4
2、3*dim,SSS,,43*dim,IIYYS,,43*dim,IIZZS,,43!*do,i,1,43!通过DO循环对数组赋值BBS(i)=1.5*(i-1)/42+3HHS(i)=(i-1)/21+3SSS(i)=BBS(i)*HHS(i)IIzzS(i)=BBS(i)*HHS(i)**3/12IIyyS(i)=HHS(i)*BBS(i)**3/12*enddo!DO循环赋值结束!*dim,BS,,42!初始化描述塔上部截面平均特征的数组*dim,HS,,42*dim,SS,,42*dim,IYYS,,42*dim,IZZS,,
3、42*do,i,1,42!通过DO循环对数组赋值hS(i)=(hhS(i)+hhS(i+1))/2bS(i)=(bbS(i)+bbS(i+1))/2SS(i)=(SSS(i)+SSS(i+1))/2IYYS(i)=(IIYYS(i)+IIYYS(i+1))/2IZZS(i)=(IIZZS(i)+IIZZS(i+1))/2*enddo!DO循环赋值结束!*dim,BBX,,80!初始化描述塔下部截面几何特征的数组*dim,HHX,,80*dim,SSX,,80*dim,IIYYX,,80*dim,IIZZX,,80!*do,i,1,8
4、0!通过DO循环对数组赋值BBX(i)=2.5HHX(i)=(i-1)*4.7164/79+5SSX(i)=BBX(i)*HHX(i)IIzzX(i)=BBX(i)*HHX(i)**3/12IIyyX(i)=HHX(i)*BBX(i)**3/12*enddo!DO循环赋值结束!*dim,BX,,79!初始化描述塔下部截面平均特征的数组*dim,HX,,79*dim,SX,,79*dim,IYYX,,79*dim,IZZX,,79*do,i,1,79!通过DO循环对数组赋值BX(i)=(BBX(i)+BBX(i+1))/2HX(i)=
5、(HHX(i)+HHX(i+1))/2SX(i)=(SSX(i)+SSX(i+1))/2IYYX(i)=(IIYYX(i)+IIYYX(i+1))/2IZZX(i)=(IIZZX(i)+IIZZX(i+1))/2*enddo!DO循环赋值结束!*do,i,1,42!定义塔上部单元实常数r,i,SS(i),IZZS(i),IYYS(i),bS(i),hS(i)*enddo!*do,i,43,114!定义塔下部单元实常数r,i,SX(i-42),IZZX(i-42),IYYX(i-42),bX(i-42),hX(i-42)*enddo!
6、定义壳单元实常数r,115,1,,,,,!定义桥面-1单元实常数r,116,4,,,,,!定义桥面-2单元实常数r,117,1,,,,,!定义桥面-3单元实常数!定义索单元实常数*dim,cable_strain_bs,,2!定义背索单元实常数数组cable_strain_bs(1)=3.978873577e-3!定义初始应变（预应力）2500cable_strain_bs(2)=3.978873577e-3!定义初始应变（预应力）2500!R,118,0.031416,cable_strain_bs(1)!定义背索单元实常数R,1
7、19,0.031416,cable_strain_bs(2)!*dim,cable_strain_zs,,11!定义主索单元实常数数组cable_strain_zs(1)=2.419155135e-3!定义主索初始应变（预应力）380cable_strain_zs(2)=2.419155135e-3!定义主索初始应变（预应力）380cable_strain_zs(3)=2.2918311805e-3!定义主索初始应变（预应力）360cable_strain_zs(4)=2.2918311805e-3!定义主索初始应变（预应力）360
8、cable_strain_zs(5)=2.228169203e-3!定义主索初始应变（预应力）350cable_strain_zs(6)=2.228169203e-3!定义主索初始应变（预应力）350cable_strain_zs(7
