地震接觸命令流 ansys

1、本人現(xiàn)在在做一個(gè)地震波的時(shí)程分析，經(jīng)過計(jì)算，是收斂的，有位移時(shí)程結(jié)果，但是沒有加速度及速度的時(shí)程結(jié)果。請各位高手指點(diǎn)。謝謝！望大家一起交流，共同進(jìn)步。命令流如下，大概需要計(jì)算1個(gè)小時(shí)。 fini /clear /COM,Structural /prep7 !單元定義 !土及樁所用單元 et,1,plane42 !平面單元 keyopt,1,3,2 !打開plane42的平面應(yīng)變選項(xiàng) r,1 !邊界所用彈簧單元 et,2,14 keyopt,2,3,2 r,2,826 !*場地土的材料模型* !*上層4.5m粉質(zhì)砂土層 mp,dens,1,2000 mp,ex,1,2.5e8 mp,nuxy,1

2、,0.4 mp,mu,1,0.3 TB,DP,1 TBdata,1,1.2e4,13 !*下層24m基巖 mp,dens,2,2600 mp,ex,2,1.8e9 mp,nuxy,2,0.25 mp,mu,2,0.4 TB,DP,2 TBdata,1,8e4,32 !*下部換土層 mp,dens,3,2600 mp,ex,3,1.8e9 mp,nuxy,3,0.25 mp,mu,3,0.4 TB,DP,3 TBdata,1,8e4,32 !*四周換上層 mp,dens,4,2000 mp,ex,4,2.5e8 mp,nuxy,4,0.4 mp,mu,4,0.3 TB,DP,4 TBdata,1

3、,1.2e4,13 !*四周換下層 mp,dens,5,2600 mp,ex,5,1.8e9 mp,nuxy,5,0.25 mp,mu,5,0.4 TB,DP,5 TBdata,1,8e4,32 !*混凝土樁、承臺及墩* mp,dens,6,2500 mp,ex,6,3e10 mp,nuxy,6,0.2 !*平面自由場地層幾何模型* !*第一層土* RECTNG,-9,-4.5,28.5,24 RECTNG,-23,-11,28.5,24 !*第二層土* RECTNG,-2,0,24,12 RECTNG,-4.5,-4,24,12 RECTNG,-9,-4.5,24,12 RECTNG,-23

4、,-11,24,0 RECTNG,-2,0,12,8 RECTNG,-4,-2,12,8 RECTNG,-9,-4,12,8 RECTNG,-11,-9,12,0 RECTNG,-9,-6,8,6 RECTNG,-9,0,4,0 RECTNG,-23,0,0,-16 !*下部換土層* RECTNG,-6,0,8,6 RECTNG,-9,0,6,4 !*四周換土層* !四周上部換土層 RECTNG,-11,-9,28.5,24 !四周下部換土層 RECTNG,-11,-9,24,12 !*樁* rectng,-4,-2,24,12 !*承臺* rectng,-4.5,-4,28.5,24 rec

5、tng,-4,-3,28.5,24 k,100,-3,28.5,0 k,101,-3,24,0 k,102,-2,24,0 k,103,-1.65,28.5,0 a,100,101,102,103 k,104,-1.65,28.5,0 k,105,-2,24,0 k,106,0,24,0 k,107,0,28.5,0 a,104,105,106,107 !*墩* rectng,-3,-1.65,38.5,28.5 !生成對稱面 FLST,3,23,5,ORDE,2 FITEM,3,1 FITEM,3,-23 ARSYM,X,P51X, , , ,0,0 !布爾運(yùn)算 !土的布爾 allsel F

6、LST,2,34,5,ORDE,4 FITEM,2,1 FITEM,2,-17 FITEM,2,24 FITEM,2,-40 AGLUE,P51X !樁、承臺及墩的布爾 allsel FLST,2,12,5,ORDE,4 FITEM,2,18 FITEM,2,-23 FITEM,2,41 FITEM,2,-46 AGLUE,P51X !有限網(wǎng)格劃分 !劃分網(wǎng)格 lsel,s,loc,y,29,38 !墩 lesize,all,10 lsel,s,loc,y,25,28 !包括部分承臺 lesize,all,5 lsel,s,loc,y,13,22 lesize,all,12 lsel,s,lo

7、c,y,8.5,11.5 lesize,all,4 lsel,s,loc,y,6.5,7.5 lsel,a,loc,y,4.5,5.5 lesize,all,2 lsel,s,loc,y,0.5,3.5 lesize,all,4 lsel,s,loc,x,-23 lsel,a,loc,x,23 lsel,r,loc,y,0,24 lesize,all,24 lsel,s,loc,x,-11 lsel,a,loc,x,11 lsel,r,loc,y,0,12 lesize,all,12 lsel,s,loc,x,-11.6,-22.5 lsel,a,loc,x,11.6,22.5 lesize,

8、all,12 lsel,s,loc,x,-9.6,-10.5 lsel,a,loc,x,9.6,10.5 lesize,all,2 lsel,s,loc,y,12,28.5 lsel,r,loc,x,4.6,8.5 lesize,all,5 lsel,s,loc,y,12,28.5 lsel,r,loc,x,-4.6,-8.5 lesize,all,5 lsel,s,loc,y,12,28.5 lsel,r,loc,x,4.1,4.4 lesize,all,1 lsel,s,loc,y,12,28.5 lsel,r,loc,x,-4.1,-4.4 lesize,all,1 lsel,s,loc

9、,y,24,28.5 lsel,r,loc,x,3.1,3.9 lesize,all,1 lsel,s,loc,y,24,28.5 lsel,r,loc,x,-3.1,-3.9 lesize,all,1 lsel,s,loc,y,24,38.5 lsel,r,loc,x,2.1,2.9 lesize,all,2 lsel,s,loc,y,24,38.5 lsel,r,loc,x,-2.1,-2.9 lesize,all,2 lsel,s,loc,y,8,28.5 lsel,r,loc,x,0.1,1.6 lesize,all,2 lsel,s,loc,y,8,28.5 lsel,r,loc,x

10、,-0.1,-1.6 lesize,all,2 lsel,s,loc,y,8,12 lsel,r,loc,x,2.1,3.9 lesize,all,3 lsel,s,loc,y,8,12 lsel,r,loc,x,-2.1,-3.9 lesize,all,3 lsel,s,loc,y,8,11 lsel,r,loc,x,4.1,5.9 lesize,all,3 lsel,s,loc,y,8,11 lsel,r,loc,x,-4.1,-5.9 lesize,all,3 lsel,s,loc,y,6,8 lsel,r,loc,x,6.1,8.9 lesize,all,3 lsel,s,loc,y,

11、6,8 lsel,r,loc,x,-6.1,-8.9 lesize,all,3 lsel,s,loc,y,6 lsel,r,loc,x,0.1,5.9 lesize,all,8 lsel,s,loc,y,6 lsel,r,loc,x,-0.1,-5.9 lesize,all,8 lsel,s,loc,y,0,4 lsel,r,loc,x,0.1,5.9 lesize,all,11 lsel,s,loc,y,0,4 lsel,r,loc,x,-0.1,-5.9 lesize,all,11 lsel,s,loc,y,-0.1,-15.4 lesize,all,8 lsel,s,loc,y,-16

12、lesize,all,25 !第一層土劃分網(wǎng)格 allsel asel,s,1,2 asel,a,24,25 aatt,1,1,1,0 !aatt,mat,real,type,esys amesh,all !第二層土劃分網(wǎng)格 allsel,all asel,s,3 asel,a,47 asel,a,49,55 asel,a,57 asel,a,58,60,2 asel,a,61,63,2 asel,a,64,68,2 asel,a,69,71,2 asel,a,72,74,2 asel,a,75 aatt,2,1,1,0 !aatt,mat,real,type,esys amesh,all !

13、下部換土層劃分網(wǎng)格 allsel,all asel,s,67,73,3 asel,a,65 aatt,3,1,1,0 !aatt,mat,real,type,esys amesh,all !四周換土層劃分網(wǎng)格 !四周上層換土層劃分網(wǎng)格 allsel,all asel,s,48 asel,a,56 aatt,4,1,1,0 !aatt,mat,real,type,esys amesh,all !四周下層換土層劃分網(wǎng)格 allsel,all asel,s,59 asel,a,62 aatt,5,1,1,0 !aatt,mat,real,type,esys amesh,all !樁、承臺及墩劃分網(wǎng)格

14、 allsel,all asel,s,4,15 aatt,6,1,1,0 !aatt,mat,real,type,esys amesh,all !*接觸分析* et,3,169 et,4,171 KEYOPT,4,12,0 !選擇接觸類型(CONTA171為STANDARD r,3,1,0,0.039,0.1 !RMORE,2.6E4 !設(shè)置實(shí)常數(shù),最大接觸摩擦26kPa !定義樁的剛性接觸面（目標(biāo)面） lsel,s,line,76 lsel,a,line,15 lsel,a,line,72 lsel,a,line,70 lsel,a,line,69 lsel,a,line,23 lsel,a

15、,line,26 lsel,a,line,162 lsel,a,line,164 lsel,a,line,161 lsel,a,line,31 lsel,a,line,168 CM,target,LINE lplot type,3 real,3 nsll,s,1 !選擇所有在已選線上的節(jié)點(diǎn) esurf,all !在已選擇的單元的自由面上產(chǎn)生新單元 !esel,s,type,3 !選擇3類型的單元 !eplot !定義土的柔性接觸面（接觸面） lsel,s,line,2 lsel,a,line,194 lsel,a,line,14 lsel,a,line,212 lsel,a,line, ,12

16、 lsel,a,line,11 lsel,a,line,193 lsel,a,line,104 lsel,a,line,219 lsel,a,line,106 lsel,a,line,207 lsel,a,line,94 CM,contact,LINE lplot type,4 real,3 nsll,s,1 !選擇所有在已選線上的節(jié)點(diǎn) esurf,all !在已選擇的單元的自由面上產(chǎn)生新單元 !esel,s,type,4 !選擇4類型的單元 !eplot ALLSEL,ALL FINISH !*求解* !*施加約束* /SOLU nsel,s,loc,x,-23 nsel,a,loc,x,2

17、3 D,ALL,all,0 nsel,s,loc,y,-16 D,ALL,uy,0 ALLSEL /replot,fast fini save !*以上幾何建模完成* !*以下開始求解* !*模態(tài)分析*為了計(jì)算基頻！ /solu antype,modal !pstres,off modopt,lanb,20,0,0,off !選擇模態(tài)分析方法 mxpand,20,yes !擴(kuò)展振型 outpr,basic,all !控制求解屏幕輸出 outres,all,all !輸出所有項(xiàng)，每一步都輸出 allsel solve !*通用后處理* !/post1 !set,list !顯示各階頻率 !set

18、,last !第一階模態(tài) !pldisp,1 !fini !*地震時(shí)程分析*（只有x方向的振動） !*將地震數(shù)據(jù)存在ff文件中* fini /config,nres,20000 !定義最大子步（只在還沒有進(jìn)行求解時(shí)才能設(shè)置） /config,fsplit,750 !由于結(jié)果文件太大(>4G，指定每個(gè)結(jié)果文件大小為750*4=3000M !/config,statu !顯示修改后的值 /solu allsel antype,trans !指定為時(shí)程分析（采用荷載步法，而非連續(xù)法） dampratio=0.239 pi=3.1415 freq1=6.0590*2*pi freq2=6.838

19、1*2*pi !得到阻尼系數(shù)（質(zhì)量系數(shù)和剛度系數(shù)） alphad,2*dampratio*freq1*freq2/(freq1+freq2 betad,2*dampratio/(freq1+freq2 trnopt,full !指定為完全瞬態(tài)分析 timint,off !關(guān)閉時(shí)間積分開關(guān) outres,basic,all !輸出基本項(xiàng)，每一項(xiàng)都輸出 nlgeom,on !打開大變形開關(guān) !ncnv,0 !即使不收斂也不中斷分析 !sstif,off time,1e-8 !設(shè)置一個(gè)十分小的時(shí)間步 kbc,1 !設(shè)置加載方式為階躍式 nsubst,2 !設(shè)置兩個(gè)子步 acel,0,9.8,0 !施加重力加速度 alls lswr,1 !把這個(gè)寫入第一子步 !solve !求解 !計(jì)算在重力作用下，結(jié)構(gòu)的地震響應(yīng) timint,on !打開時(shí)間積分效應(yīng) !outres,basic,all !輸出基本項(xiàng)，每一項(xiàng)都輸出 kbc,0 !設(shè)置加載方式為