时程分析实例

时程分析实例 finish/clear/filename,ex6-1 !设定工作名/title, The Dynamic Analysis of the Bridge !设定工作标题/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)=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,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)=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)=(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)=3.978873577e-3 !定义初始应变（预应力）2500cable_strain_bs(2)=3.978873577e-3 !定义初始应变（预应力）2500!R,118,0.031416,cable_strain_bs(1) !定义背索单元实常数R,119,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 !定义主索初始应变（预应力）360cable_strain_zs(5)=2.228169203e-3 !定义主索初始应变（预应力）350cable_strain_zs(6)=2.228169203e-3 !定义主索初始应变（预应力）350cable_strain_zs(7)=1.909859317e-3 !定义主索初始应变（预应力）300cable_strain_zs(8)=1.65521141e-3 !定义主索初始应变（预应力）260cable_strain_zs(9)=1.65521141e-3 !定义主索初始应变（预应力）260cable_strain_zs(10)=1.40056340e-3 !定义主索初始应变（预应力）220cable_strain_zs(11)=1.273239545e-3 !定义主索初始应变（预应力）200!*do,i,120,130 !定义主索实常数R,i,0.007853982,cable_strain_zs(i-119) *enddo!定义桥面与塔连接单元实常数r,131,16,21.3333,21.3333,4,4 !定义桥面与塔的连接梁单元实常数!定义材料参数mp,ex,1,3.4e10 !定义塔上部材料特性mp,prxy,1,0.3mp,dens,1,2.5e3!MP,EX,2,3.5e10 !定义桥面材料特性mp,prxy,2,0.3MP,DENS,2,2.5e3!MP,EX,3,2e11 !定义背索材料特性mp,prxy,3,0.3MP,DENS,3,7.85e3!*!生成关键点!*!定义桥塔的关键点*do,i,54.5,33.5,-0.5 !生成塔上部关键点k,110-i*2,-(82.1495+i)/7.403322,i,0*enddo!*do,i,33,0,-0.5 !生成塔下部关键点k,110-i*2,-15.6213*i/33.5,i,7.4839*i/33.5-7.4839*enddo*do,i,-0.5,-2.5,-0.5k,110-i*2,-2.7978*i/6,i,-i/2-7.4839*enddo*do,i,33,0,-0.5k,182-i*2,-15.6213*i/33.5,i,-7.4839*i/33.5+7.4839*enddo*do,i,-0.5,-2.5,-0.5k,182-i*2,-2.7978*i/6,i,i/2+7.4839*enddo!定义桥面关键点(其中配跨为曲线)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, !指定线的材料编号、实常数编号、单元类型号le size,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,-5.6695,-70,76,90,64,118.09989655,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!*!反应谱分析模块!*/solu Antype,spectr !设定谱分析类型Spopt,sprs,20,yes !设定谱分析选项Svtype,2 !设定单点响应谱类型Sed,1,1, !设定激励方向Freq,0.1181,0.1228,0.1436,0.1576,0.1629,0.1639,0.1761,0.2197,0.2273,0.2358,0.2588,0.2674,0.3209,0.3478,0. 3739,0.4002,0.4938,0.7503,1.3831,1.5370 !设定频率数值Sv,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.07996,0.0662,0.0454,0.0262,0.0238 !设定频谱数值SolveSave,rst2,dbFinish!*模态扩展*/solu Antype,modal !设定模态分析类型Expass,on !打开模态扩展Mxpand,20,yes,0.005 !设定模态扩展选项SolveSave,rst3,dbFinish!*合并模态*/solu Antype,spectr !设定谱分析类型Srss,0.02,disp !设定模态合并选项SolveSave,rst4,dbFinish/post26 !进入时间历程后处理器查看结果Nsol,2,2,u,x,ux_2 !定义节点2 X方向位移为变量2Nsol,3,2,u,y,uy_3 !定义节点2 Y方向位移为变量3NSOL,4,400,U,X, UX_4 !定义节点400 X方向位移为变量4NSOL,5,400,U,Y, UY_5 !定义节点400 Y方向位移为变量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.txt文件数据到horizon数组(e9.3,e11.3) !设定数据读入格式*END !宏结束标识 /INPUT,ansuitmp *CREATE,ansuitmp *VREAD,vertical(1,1),