地震接触命令流 ansys

1、本人现在在做一个地震波的时程分析，经过计算，是收敛的，有位移时程结果，但是没有加速度及速度的时程结果。请各位高手指点。谢谢！望大家一起交流，共同进步。命令流如下，大概需要计算1个小时。 fini /clear /COM,Structural /prep7 !单元定义 !土及桩所用单元 et,1,plane42 !平面单元 keyopt,1,3,2 !打开plane42的平面应变选项 r,1 !边界所用弹簧单元 et,2,14 keyopt,2,3,2 r,2,826 !*场地土的材料模型* !*上层4.5m粉质砂土层 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 !布尔运算 !土的布尔 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 !有限网格划分 !划分网格 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 !第一层土划分网格 allsel asel,s,1,2 asel,a,24,25 aatt,1,1,1,0 !aatt,mat,real,type,esys amesh,all !第二层土划分网格 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、下部换土层划分网格 allsel,all asel,s,67,73,3 asel,a,65 aatt,3,1,1,0 !aatt,mat,real,type,esys amesh,all !四周换土层划分网格 !四周上层换土层划分网格 allsel,all asel,s,48 asel,a,56 aatt,4,1,1,0 !aatt,mat,real,type,esys amesh,all !四周下层换土层划分网格 allsel,all asel,s,59 asel,a,62 aatt,5,1,1,0 !aatt,mat,real,type,esys amesh,all !桩、承台及墩划分网格

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 !设置实常数,最大接触摩擦26kPa !定义桩的刚性接触面（目标面） 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 !选择所有在已选线上的节点 esurf,all !在已选择的单元的自由面上产生新单元 !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 !选择所有在已选线上的节点 esurf,all !在已选择的单元的自由面上产生新单元 !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 !*以上几何建模完成* !*以下开始求解* !*模态分析*为了计算基频！ /solu antype,modal !pstres,off modopt,lanb,20,0,0,off !选择模态分析方法 mxpand,20,yes !扩展振型 outpr,basic,all !控制求解屏幕输出 outres,all,all !输出所有项，每一步都输出 allsel solve !*通用后处理* !/post1 !set,list !显示各阶频率 !set

18、,last !第一阶模态 !pldisp,1 !fini !*地震时程分析*（只有x方向的振动） !*将地震数据存在ff文件中* fini /config,nres,20000 !定义最大子步（只在还没有进行求解时才能设置） /config,fsplit,750 !由于结果文件太大(>4G，指定每个结果文件大小为750*4=3000M !/config,statu !显示修改后的值 /solu allsel antype,trans !指定为时程分析（采用荷载步法，而非连续法） dampratio=0.239 pi=3.1415 freq1=6.0590*2*pi freq2=6.838

19、1*2*pi !得到阻尼系数（质量系数和刚度系数） alphad,2*dampratio*freq1*freq2/(freq1+freq2 betad,2*dampratio/(freq1+freq2 trnopt,full !指定为完全瞬态分析 timint,off !关闭时间积分开关 outres,basic,all !输出基本项，每一项都输出 nlgeom,on !打开大变形开关 !ncnv,0 !即使不收敛也不中断分析 !sstif,off time,1e-8 !设置一个十分小的时间步 kbc,1 !设置加载方式为阶跃式 nsubst,2 !设置两个子步 acel,0,9.8,0 !施加重力加速度 alls lswr,1 !把这个写入第一子步 !solve !求解 !计算在重力作用下，结构的地震响应 timint,on !打开时间积分效应 !outres,basic,all !输出基本项，每一项都输出 kbc,0 !设置加载方式为