时程分析实例

时程分析实例时程分析实例时程分析实例资料仅供参考文件编号：2022年4月时程分析实例版本号：A修改号：1页次：1.0审核：批准:发布日期：时程分析实例finish/clear/filename,ex6-1!设定工作名/title,TheDynamicAnalysisoftheBridge!设定工作标题/prep7!进入前处理模块!**********************************!定义单元类型及单元实常数!**********************************Et,1,beam4!定义单元类型Et,2,shell63Et,3,link10!定义梁单元实常数*dim,BBS,,43!初始化描述塔上部截面几何特征的数组*dim,HHS,,43*dim,SSS,,43*dim,IIYYS,,43*dim,IIZZS,,43!*do,i,1,43!通过DO循环对数组赋值BBS(i)=*(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,,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,80!通过DO循环对数组赋值BBX(i)=HHX(i)=(i-1)*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)=(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!定义壳单元实常数r,115,1,,,,,!定义桥面-1单元实常数r,116,4,,,,,!定义桥面-2单元实常数r,117,1,,,,,!定义桥面-3单元实常数!定义索单元实常数*dim,cable_strain_bs,,2!定义背索单元实常数数组cable_strain_bs(1)=!定义初始应变（预应力）2500cable_strain_bs(2)=!定义初始应变（预应力）2500!R,118,,cable_strain_bs(1)!定义背索单元实常数R,119,,cable_strain_bs(2)!*dim,cable_strain_zs,,11!定义主索单元实常数数组cable_strain_zs(1)=!定义主索初始应变（预应力）380cable_strain_zs(2)=!定义主索初始应变（预应力）380cable_strain_zs(3)=!定义主索初始应变（预应力）360cable_strain_zs(4)=!定义主索初始应变（预应力）360cable_strain_zs(5)=!定义主索初始应变（预应力）350cable_strain_zs(6)=!定义主索初始应变（预应力）350cable_strain_zs(7)=!定义主索初始应变（预应力）300cable_strain_zs(8)=!定义主索初始应变（预应力）260cable_strain_zs(9)=!定义主索初始应变（预应力）260cable_strain_zs(10)=!定义主索初始应变（预应力）220cable_strain_zs(11)=!定义主索初始应变（预应力）200!*do,i,120,130!定义主索实常数R,i,,cable_strain_zs(i-119)*enddo!定义桥面与塔连接单元实常数r,131,16,,,4,4!定义桥面与塔的连接梁单元实常数!定义材料参数mp,ex,1,!定义塔上部材料特性mp,prxy,1,mp,dens,1,!MP,EX,2,!定义桥面材料特性mp,prxy,2,MP,DENS,2,!MP,EX,3,2e11!定义背索材料特性mp,prxy,3,MP,DENS,3,!**********************************!生成关键点!**********************************!定义桥塔的关键点*do,i,,,!生成塔上部关键点k,110-i*2,-+i)/,i,0*enddo!*do,i,33,0,!生成塔下部关键点k,110-i*2,*i/,i,*i/定义桥面关键点(其中配跨为曲线)k,192,-2,0,-6!生成左配跨关键点k,193,0,0,-6k,194,2,0,-6k,195,9,0,-6!*do,i,196,205!生成左桥面关键点k,i,(i-195)*9+9,0,-6*enddok,206,105,0,-6!k,207,-2,0,6!生成右配跨关键点k,208,0,0,6k,209,2,0,6k,210,9,0,6!*do,i,211,220!生成右桥面关键点k,i,(i-210)*9+9,0,6*enddok,221,105,0,6!定义背索锚碇关键点k,222,-31,0,14k,223,-31,0,-14!**********************************!划分网格，生成单元!***********************************do,i,1,42!划分塔上部网格l,i,i+1!在两个关键点之间生成线lsel,u,type,,1,131!不选已定义属性的线latt,1,i,1,,,,!指定线的材料编号、实常数编号、单元类型号lesize,all,,,1,,1!设定线的分割数lmesh,all,all!对线划分网格alls!选择全部图元*enddo!DO循环结束标志!*do,i,43,114!划分塔下部－近朝网格l,i,i+1lsel,u,type,,1,131latt,1,i,1,,,,lesize,all,,,1,,1lmesh,all,allalls*enddo!l,43,116!划分塔下部－远朝网格lsel,u,type,,1,131latt,1,43,1,,,,lesize,all,,,1,,1lmesh,all,allalls!*do,i,116,186l,i,i+1lsel,u,type,,1,131latt,1,i-72,1,,,,lesize,all,,,1,,1lmesh,all,allalls*enddo!l,9,222!划分背索网格lsel,u,type,,1,131latt,3,118,3,,,,lesize,all,,,1,,1lmesh,all,allalls!l,9,223lsel,u,type,,1,131latt,3,119,3,,,,lesize,all,,,1,,1lmesh,all,allalls!*do,i,1,11!划分主索网格l,i*4+1,209-il,i*4+1,221-ilsel,u,type,,1,131latt,3,i+119,3,,,,lesize,all,,,1,,1lmesh,all,allalls*enddo!l,110,193!划分桥面与塔的连接处网格l,182,208lsel,u,type,,1,131latt,3,131,1,,,,lesize,all,,,1,,1lmesh,all,allalls!wproat,0,-90,0!划分桥面网格cyl4,,-70,76,90,64,,0wproat,0,90,0!a,188,192,207,191!划分左桥面网格lsel,u,type,,1,131aatt,2,115,2,,,,lesize,all,2,,,,2lesize,213,2,,,,2lesize,214,2,,,,2lesize,215,2,,,,2lesize,216,2,,,,2lesize,217,,,,4,2lesize,218,2,,,,2lesize,219,,,,4,2allsamesh,all!a,192,193,208,207!划分桥面交界处网格a,193,194,209,208lsel,u,type,,1,131aatt,2,116,2,,,,lesize,all,2,,,,1allsamesh,all!*do,i,1,12!划分右桥面网格a,i+193,i+194,i+209,i+208*enddolsel,u,type,,1,131aatt,2,117,2,,,,lesize,all,2,,,,1allsamesh,allsavefinish!**********************************!加载、求解!**********************************/SOLUd,115,all!在节点施加位移约束d,187,alld,188,alld,189,alld,235,uyd,239,uyd,280,uyd,292,uyd,735,uyd,739,uyacel,,10!施加重力加速度!**********************************!模态分析模块!**********************************antype,modal!设定模态分析类型modopt,subsp,20!设定模态分析选项solvesave,rst11,dbfinish!**********************************!反应谱分析模块!**********************************/soluAntype,spectr!设定谱分析类型Spopt,sprs,20,yes!设定谱分析选项Svtype,2!设定单点响应谱类型Sed,1,1,!设定激励方向Freq,,,,,,,,,,,,,,,0.3739,,,,,!设定频率数值Sv,,,,,,,,,,,,,,,,,,,,!设定频谱数值SolveSave,rst2,dbFinish!******模态扩展******/soluAntype,modal!设定模态分析类型Expass,on!打开模态扩展Mxpand,20,,,yes,!设定模态扩展选项SolveSave,rst3,dbFinish!*******合并模态******/soluAntype,spectr!设定谱分析类型Srss,,disp!设定模态合并选项SolveSave,rst4,dbFinish/post26!进入时间历程后处理器查看结果Nsol,2,2,u,x,ux_2!定义节点2X方向位移为变量2Nsol,3,2,u,y,uy_3!定义节点2Y方向位移为变量3NSOL,4,400,U,X,UX_4!定义节点400X方向位移为变量4NSOL,5,400,U,Y,UY_5!定义节点400Y方向位移为变量5STORE,MERGE!存储定义的变量PLVAR,2,3,4,5,!图形显示变量2、3、4、5Finish!**********************************!地震波瞬态分析模块!***********************************DIM,horizon,ARRAY,2,50,0,,,!初始化horizon数组*DIM,vertical,ARRAY,2,50,0,,,!初始化vertical数组*CREATE,ansuitmp!宏开始标识*VREAD,horizon(1,1),'horizon','txt','',50,,,,,,!读入文件数据到horizon数组,!设定数据读入格式*END!宏结束标识/INPUT,ansuitmp*CREATE,ansuitmp*VREAD,vertical(1,1),'vertical','txt','',50,,,,,,!读入文件数据到vertica