時程分析實例

1、時程分析實例finish/clear/filename,ex6-1!設定工作名/title,TheDynamicAnalysisoftheBridge!設定工作標題/prep7塊F*!定義單元類型及單元實常數(shù)!*Et,1,beam4Et,2,shell63Et,3,link10!定義梁單元實常數(shù)*dim,BBS,43幾何特征的數(shù)組*dim,HHS,43*dim,SSS,43!進入前處理模!定義單元類型!初始化描述塔上部截面*dim,IIYYS,43*dim,IIZZS,43*do,i,1,43BBS=1.5*(i-1)/42+3HHS(i)=(i-1)/21+3SSS(i)=BBS(i)*HH

2、S(i)IIzzS(i)=BBS(i)*HHS(i)*3/12IIyyS(i)=HHS(i)*BBS(i)*3/12*enddo!*dim,BS,42特征的數(shù)組*dim,HS,42*dim,SS,42*dim,IYYS,42*dim,IZZS,42*do,i,1,42hS(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)/2!通過DO循環(huán)對數(shù)組賦值!DO循環(huán)賦值結(jié)束!初始化描述塔上部截面平均!通過DO循環(huán)對數(shù)組賦值IZZS(i)=(IIZZS(i)+

3、IIZZS(i+1)/2*enddo!DO循環(huán)賦值結(jié)束*dim,BBX,80特征的數(shù)組*dim,HHX,80*dim,SSX,80*dim,IIYYX,80*dim,IIZZX,80!*do,i,1,80值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!*dim,BX,79征的數(shù)組*dim,HX,79*dim,SX,79!初始化描述塔下部截面幾何!通過DO循環(huán)對數(shù)組賦!DO循環(huán)賦值結(jié)束!初始化描述塔下部截面平均特*

4、dim,IYYX,79*dim,IZZX,79* do,i,1,79!通過DO循環(huán)對數(shù)組賦值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循環(huán)賦值結(jié)束!* do,i,1,42!定義塔上部單元實常數(shù)r,i,SS(i),IZZS(i),IYYS(i),bS(i),hS(i)*enddo!* do,i,43,114!定義塔下部單元實常數(shù)r,i,SX(i-42),IZ

5、ZX(i-42),IYYX(i-42),bX(i-42),hX(i-42)*enddo!定義殼單元實常數(shù)r,115,1,!定義橋面-1單元實常數(shù)r,116,4,!定義橋面-2單元實常數(shù)r,117,1,!定義橋面-3單元實常數(shù)!定義索單元實常數(shù)*dim,cable_strain_bs,2cable_strain_bs(1)=3.978873577e-3cable_strain_bs(2)=3.978873577e-3!R,118,0.031416,cable_strain_bs(1)R,119,0.031416,cable_strain_bs(2)!*dim,cable_strain_zs,11

6、數(shù)組cable_strain_zs(1)=2.419155135e-3380cable_strain_zs=2.419155135e-3380cable_strain_zs(3)=2.2918311805e-3360cable_strain_zs(4)=2.2918311805e-3360cable_strain_zs(5)=2.228169203e-3!定義背索單元實常數(shù)數(shù)組!定義初始應變（預應力）2500!定義初始應變（預應力）2500!定義背索單元實常數(shù)!定義主索單元實常數(shù)!定義主索初始應變（預應力）!定義主索初始應變（預應力）!定義主索初始應變（預應力）!定義主索初始應變（預應力）!定

7、義主索初始應變（預應力）350!定義主索初始應變（預應力）!定義主索初始應變（預應力）!定義主索初始應變（預應力）!定義主索初始應變（預應力）!定義主索初始應變（預應力）!定義主索初始應變（預應力）!定義主索實常數(shù)!定義橋面與塔的連接梁單!定義塔上部材料特性cable_strain_zs(6)=2.228169203e-3350cable_strain_zs(7)=1.909859317e-3300cable_strain_zs(8)=1.65521141e-3260cable_strain_zs(9)=1.65521141e-3260cable_strain_zs(10)=1.4005634

8、0e-3220cable_strain_zs(11)=1.273239545e-3200!*do,i,120,130R,i,0.007853982,cable_strain_zs(i-119)*enddo!定義橋面與塔連接單元實常數(shù)r,131,16,21.3333,21.3333,4,4元實常數(shù)!定義材料參數(shù)mp,ex,1,3.4e10mp,prxy,1,0.3!定義橋面材料特性mp,dens,1,2.5e3!定義背索材料特性MP,EX,2,3.5e10mp,prxy,2,0.3MP,DENS,2,2.5e3!MP,EX,3,2e11mp,prxy,3,0.3MP,DENS,3,7.85e3產(chǎn)

9、*!生成關(guān)鍵點I*!定義橋塔的關(guān)鍵點*do,i,54.5,33.5,-0.5!生成塔上部關(guān)鍵點k,110-i*2,-(82.1495+i)/7.403322,i,0*do,i,33,0,-0.5*enddo!生成塔下部關(guān)鍵點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

10、,i,-0.5,-2.5,-0.5k,182-i*2,-2.7978*i/6,i,i/2+7.4839*enddo!定義橋面關(guān)鍵點(其中配跨為曲線)k,192,-2,0,-6!生成左配跨關(guān)鍵點k,193,0,0,-6k,194,2,0,-6k,195,9,0,-6!*do,i,196,205!生成左橋面關(guān)鍵點k,i,(i-195)*9+9,0,-6*enddok,206,105,0,-6!k,207,-2,0,6!生成右配跨關(guān)鍵點k,208,0,0,6k,209,2,0,6*do,i,211,220!生成右橋面關(guān)鍵點k,i,(i-210)*9+9,0,6*enddok,221,105,0,6!

11、定義背索錨碇關(guān)鍵點k,222,-31,0,14k,223,-31,0,-14產(chǎn)*!劃分網(wǎng)格，生成單元!*!劃分塔上部網(wǎng)格!在兩個關(guān)鍵點之間生成線!不選已定義屬性的線!指定線的材料編號、實常數(shù)編號、單*do,i,1,42l,i,i+1lsel,u,type,1,131latt,1,i,1,元類型號lesize,all,1,1!設定線的分割lmesh,all,all!對線劃分網(wǎng)格!選擇全部圖元!DO循環(huán)結(jié)束標志!劃分塔下部-近朝網(wǎng)格!劃分塔下部-遠朝網(wǎng)格alls*enddo!*do,i,43,114l,i,i+1lsel,u,type,1,131latt,1,i,1,lesize,all,1,1l

12、mesh,all,allalls*enddo!l,43,116lsel,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!劃分背索網(wǎng)格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

13、,all,1,1lmesh,all,all*do,i,1,11alls!劃分主索網(wǎng)格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!劃分橋面與塔的連接處網(wǎng)格l,182,208lsel,u,type,1,131latt,3,131,1,lesize,all,1,1lmesh,all,allalls!wproat,0,-90,0!劃分橋面網(wǎng)格cyl4,-5.6695,-70,76,90,64,118.09989655,0wproat,0

14、,90,0a,188,192,207,191lsel,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處網(wǎng)格a,193,194,209,208lsel,u,type,1,131aatt,2,116,2,lesize,all,2,1alls!劃分左橋面網(wǎng)格!劃分橋面交界amesh,all!劃分右橋面網(wǎng)格a,i+1

15、93,i+194,i+209,i+208*enddolsel,u,type,1,131aatt,2,117,2,lesize,all,2,1allsamesh,allsavefinish產(chǎn)*!加載、求解!*/SOLUd,115,all!在節(jié)點施加位移約束d,187,alld,188,alld,189,alld,235,uyd,239,uyd,280,uyd,292,uyd,735,uyd,739,uyacel,10產(chǎn)*!模態(tài)分析模塊!*antype,modalmodopt,subsp,20solvesave,rst11,dbfinish!*!反應譜分析模塊!*/soluAntype,spect

16、rSpopt,sprs,20,yesSvtype,2!施加重力加速度!設定模態(tài)分析類型!設定模態(tài)分析選項!設定譜分析類型!設定譜分析選項!設定單點響應譜類型!設定激勵方向Sed,1,1,Freq,0.1181,0.1228,0.1436,0.1576,0.1629,0.1639,0.1761,0.2197,0.22730.2358,0.2588,0.2674,0.3209,0.3478,0.3739,0.4002,0.4938,0.7503,1.3831,1.5370設定頻率數(shù)值Sv,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08,

17、0.08,0.08,0.08,0.08,0.07996,0.0662,0.0454,0.0262,0.0238!設定頻譜數(shù)值SolveSave,rst2,dbFinish*模態(tài)擴展*/soluAntype,modalExpass,onMxpand,20,yes,0.005Solve!設定模態(tài)分析類型!打開模態(tài)擴展!設定模態(tài)擴展選項Save,rst3,dbFinishI*合并模態(tài)*/solu!設定譜分析類型!設定模態(tài)合并選項!進入時間歷程后處理器!定義節(jié)點2X方向位移為變量2!定義節(jié)點2Y方向位移為變量3!定義節(jié)點400X方向位移為變量4!定義節(jié)點400Y方向位移為變量5!存儲定義的變量!圖形顯

18、示變量2、3、4、5!初始化horizon數(shù)組!初始化vertical數(shù)組!宏開始標識Antype,spectrSrss,0.02,dispSolveSave,rst4,dbFinish/post26查看結(jié)果Nsol,2,2,u,x,ux_2Nsol,3,2,u,y,uy_3NSOL,4,400,U,X,UX_4NSOL,5,400,U,Y,UY_5STORE,MERGEPLVAR,2,3,4,5,FinishF*!地震波瞬態(tài)分析模塊!*DIM,horizon,ARRAY!,50,0,*DIM,vertical,ARRA,2,50,0,*CREATE,ansuitmp*VREAD,horizon(1,1),'horizon','txt','',50,犢入horizon.txt文件數(shù)據(jù)到horizon數(shù)組(e9.3,e11.3)*END!設定數(shù)據(jù)讀入格式!宏結(jié)束標識/INPUT,ansuitmp*CREATE,ansuitmp*VREAD,vertical(1,1),'vertical&a