地震接触命令流-ansys

本人现在在做一个地震波的时程分析，经过计算，是收敛的，有位移时程结果，但是没有加速度及速度的时程结果。请各位高手指点。谢谢！望大家一起交流，共同进步。命令流如下，大概需要计算1个小时。fini/clear/COM,Structural/prep7!单元定义!土及桩所用单元et,1,plane42!平面单元keyopt,1,3,2!打开plane42的平面应变选项r,1!边界所用弹簧单元et,2,14keyopt,2,3,2r,2,826!***场地土的材料模型**!***上层4.5m粉质砂土层mp,dens,1,2000mp,ex,1,2.5e8mp,nuxy,1,0.4mp,mu,1,0.3TB,DP,1TBdata,1,1.2e4,13!***下层24m基岩mp,dens,2,2600mp,ex,2,1.8e9mp,nuxy,2,0.25mp,mu,2,0.4TB,DP,2TBdata,1,8e4,32!***下部换土层mp,dens,3,2600mp,ex,3,1.8e9mp,nuxy,3,0.25mp,mu,3,0.4TB,DP,3TBdata,1,8e4,32!***四周换上层mp,dens,4,2000mp,ex,4,2.5e8mp,nuxy,4,0.4mp,mu,4,0.3TB,DP,4TBdata,1,1.2e4,13!***四周换下层mp,dens,5,2600mp,ex,5,1.8e9mp,nuxy,5,0.25mp,mu,5,0.4TB,DP,5TBdata,1,8e4,32!***混凝土桩、承台及墩****mp,dens,6,2500mp,ex,6,3e10mp,nuxy,6,0.2!***平面自由场地层几何模型**!***第一层土***RECTNG,-9,-4.5,28.5,24RECTNG,-23,-11,28.5,24!***第二层土***RECTNG,-2,0,24,12RECTNG,-4.5,-4,24,12RECTNG,-9,-4.5,24,12RECTNG,-23,-11,24,0RECTNG,-2,0,12,8RECTNG,-4,-2,12,8RECTNG,-9,-4,12,8RECTNG,-11,-9,12,0RECTNG,-9,-6,8,6RECTNG,-9,0,4,0RECTNG,-23,0,0,-16!***下部换土层***RECTNG,-6,0,8,6RECTNG,-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,24rectng,-4,-3,28.5,24k,100,-3,28.5,0k,101,-3,24,0k,102,-2,24,0k,103,-1.65,28.5,0a,100,101,102,103k,104,-1.65,28.5,0k,105,-2,24,0k,106,0,24,0k,107,0,28.5,0a,104,105,106,107!**墩*******rectng,-3,-1.65,38.5,28.5!生成对称面FLST,3,23,5,ORDE,2FITEM,3,1FITEM,3,-23ARSYM,X,P51X,,,,0,0!布尔运算!土的布尔allselFLST,2,34,5,ORDE,4FITEM,2,1FITEM,2,-17FITEM,2,24FITEM,2,-40AGLUE,P51X!桩、承台及墩的布尔allselFLST,2,12,5,ORDE,4FITEM,2,18FITEM,2,-23FITEM,2,41FITEM,2,-46AGLUE,P51X!有限网格划分!划分网格lsel,s,loc,y,29,38!墩lesize,all,,,10lsel,s,loc,y,25,28!包括部分承台lesize,all,,,5lsel,s,loc,y,13,22lesize,all,,,12lsel,s,loc,y,8.5,11.5lesize,all,,,4lsel,s,loc,y,6.5,7.5lsel,a,loc,y,4.5,5.5lesize,all,,,2lsel,s,loc,y,0.5,3.5lesize,all,,,4lsel,s,loc,x,-23lsel,a,loc,x,23lsel,r,loc,y,0,24lesize,all,,,24lsel,s,loc,x,-11lsel,a,loc,x,11lsel,r,loc,y,0,12lesize,all,,,12lsel,s,loc,x,-11.6,-22.5lsel,a,loc,x,11.6,22.5lesize,all,,,12lsel,s,loc,x,-9.6,-10.5lsel,a,loc,x,9.6,10.5lesize,all,,,2lsel,s,loc,y,12,28.5lsel,r,loc,x,4.6,8.5lesize,all,,,5lsel,s,loc,y,12,28.5lsel,r,loc,x,-4.6,-8.5lesize,all,,,5lsel,s,loc,y,12,28.5lsel,r,loc,x,4.1,4.4lesize,all,,,1lsel,s,loc,y,12,28.5lsel,r,loc,x,-4.1,-4.4lesize,all,,,1lsel,s,loc,y,24,28.5lsel,r,loc,x,3.1,3.9lesize,all,,,1lsel,s,loc,y,24,28.5lsel,r,loc,x,-3.1,-3.9lesize,all,,,1lsel,s,loc,y,24,38.5lsel,r,loc,x,2.1,2.9lesize,all,,,2lsel,s,loc,y,24,38.5lsel,r,loc,x,-2.1,-2.9lesize,all,,,2lsel,s,loc,y,8,28.5lsel,r,loc,x,0.1,1.6lesize,all,,,2lsel,s,loc,y,8,28.5lsel,r,loc,x,-0.1,-1.6lesize,all,,,2lsel,s,loc,y,8,12lsel,r,loc,x,2.1,3.9lesize,all,,,3lsel,s,loc,y,8,12lsel,r,loc,x,-2.1,-3.9lesize,all,,,3lsel,s,loc,y,8,11lsel,r,loc,x,4.1,5.9lesize,all,,,3lsel,s,loc,y,8,11lsel,r,loc,x,-4.1,-5.9lesize,all,,,3lsel,s,loc,y,6,8lsel,r,loc,x,6.1,8.9lesize,all,,,3lsel,s,loc,y,6,8lsel,r,loc,x,-6.1,-8.9lesize,all,,,3lsel,s,loc,y,6lsel,r,loc,x,0.1,5.9lesize,all,,,8lsel,s,loc,y,6lsel,r,loc,x,-0.1,-5.9lesize,all,,,8lsel,s,loc,y,0,4lsel,r,loc,x,0.1,5.9lesize,all,,,11lsel,s,loc,y,0,4lsel,r,loc,x,-0.1,-5.9lesize,all,,,11lsel,s,loc,y,-0.1,-15.4lesize,all,,,8lsel,s,loc,y,-16lesize,all,,,25!第一层土划分网格allselasel,s,,,1,2asel,a,,,24,25aatt,1,1,1,0!aatt,mat,real,type,esysamesh,all!第二层土划分网格allsel,allasel,s,,,3asel,a,,,47asel,a,,,49,55asel,a,,,57asel,a,,,58,60,2asel,a,,,61,63,2asel,a,,,64,68,2asel,a,,,69,71,2asel,a,,,72,74,2asel,a,,,75aatt,2,1,1,0!aatt,mat,real,type,esysamesh,all!下部换土层划分网格allsel,allasel,s,,,67,73,3asel,a,,,65aatt,3,1,1,0!aatt,mat,real,type,esysamesh,all!四周换土层划分网格!四周上层换土层划分网格allsel,allasel,s,,,48asel,a,,,56aatt,4,1,1,0!aatt,mat,real,type,esysamesh,all!四周下层换土层划分网格allsel,allasel,s,,,59asel,a,,,62aatt,5,1,1,0!aatt,mat,real,type,esysamesh,all!桩、承台及墩划分网格allsel,allasel,s,,,4,15aatt,6,1,1,0!aatt,mat,real,type,esysamesh,all!******接触分析******et,3,169et,4,171KEYOPT,4,12,0!选择接触类型(CONTA171)为STANDARDr,3,1,0,0.039,0.1!RMORE,,,2.6E4!设置实常数,最大接触摩擦26kPa!定义桩的刚性接触面（目标面）lsel,s,line,,76lsel,a,line,,15lsel,a,line,,72lsel,a,line,,70lsel,a,line,,69lsel,a,line,,23lsel,a,line,,26lsel,a,line,,162lsel,a,line,,164lsel,a,line,,161lsel,a,line,,31lsel,a,line,,168CM,target,LINElplottype,3real,3nsll,s,1!选择所有在已选线上的节点esurf,all!在已选择的单元的自由面上产生新单元!esel,s,type,,3!选择3＃类型的单元!eplot!定义土的柔性接触面（接触面）lsel,s,line,,2lsel,a,line,,194lsel,a,line,,14lsel,a,line,,212lsel,a,line,,12lsel,a,line,,11lsel,a,line,,193lsel,a,line,,104lsel,a,line,,219lsel,a,line,,106lsel,a,line,,207lsel,a,line,,94CM,contact,LINElplottype,4real,3nsll,s,1!选择所有在已选线上的节点esurf,all!在已选择的单元的自由面上产生新单元!esel,s,type,,4!选择4＃类型的单元!eplotALLSEL,ALLFINISH!****求解*******!**施加约束**/SOLUnsel,s,loc,x,-23nsel,a,loc,x,23D,ALL,all,0nsel,s,loc,y,-16D,ALL,uy,0ALLSEL/replot,fastfinisave!****以上几何建模完成*******!****以下开始求解******!***模态分析***为了计算基频！/soluantype,modal!pstres,offmodopt,lanb,20,0,0,,off!选择模态分析方法mxpand,20,,,yes!扩展振型outpr,basic,all!控制求解屏幕输出outres,all,all!输出所有项，每一步都输出allselsolve!***通用后处理***!/post1!set,list!显示各阶频率!set,last!第一阶模态!pldisp,