| 查看: 209 | 回复: 0 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】用surf154单元施加面载荷进行的土基瞬态动力学分析遇到了问题
|
|||
|
参考一些相关surf154方面的资料和书上的例题，做了下面的一个例子，问题描述：对分层土基的某块面积上施加随时间变化的载荷，求解地面上任意一点的动力响应。 注：选取4800X16100X3000(mm^3)的土体作为对象进行分析，需在工作目录下建立两个txt文件，压力载荷文件fxyz5row.txt和时间文件time5row.txt，数据分别如下： 压力载荷fxyz5row.txt的内容： 6.708710E+04 3.600730E+04 4.569551E+04 5.364202E+04 6.217627E+04 时间time5row.txt的内容： 2.683363E-03 1.001789E-02 1.735242E-02 2.468694E-02 3.202147E-02 命令流如下： /PREP7 ET,1,SOLID45 et,2,mesh200,6 et,3,154 /units,user,1000 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,3.94 MPDATA,PRXY,1,,0.35 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,1.828e-6 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,2,,20.6 MPDATA,PRXY,2,,0.3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,2.16e-6 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,3,,500 MPDATA,PRXY,3,,0.33 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,3,,2.16e-6 *SET,length_x1,2400 *SET,length_x2,13700 *SET,width_y1 , -2000 *SET,WIDTH_Y2 , 2000 *SET,heigth_z1,3000 *SET,HEIGTH_Z0, 1576 *SET,LENGTH_X1 , -2400 *SET,length,length_x2-length_x1 PARSAV,ALL,'para_canshu1',' ',' ' K,1,length_x1,width_y2,heigth_z0, K,2,length_x1,width_y1,heigth_z0, K,3,length_x1,width_y1,heigth_z0+heigth_z1, K,4,length_x1,width_y2,heigth_z0+heigth_z1, A,1,2,3,4 meshkey,1 TYPE, 2 MAT,1 lesize,1,,,10 lesize,2,,,10 lesize,3,,,10 lesize,4,,,10 amesh,1 *do,i,1,50 k,,-2400+322*i,0,1576 *enddo K,1000,length_x1,0,heigth_z0, l,1000,5 *do,i,1,49 l,4+i,5+i *enddo *do,i,1,50 lesize,4+i,,,1 *enddo TYPE, 1 MAT, 1 extopt,aclear,1 *do,i,1,50 asel,s,loc,x,-2400+322*(i-1) vdrag,all,,,,,,4+i *enddo csys,0 nsel,s,loc,z,2477,3077 esln,s mpchg,2,all nsel,s,loc,z,3675,4576 esln,s mpchg,3,all allsel,all nsel,s,loc,x,11726.9,12946.1 nsel,r,loc,y,77.9,1478.1 nsel,r,loc,z,1576 !esln,s !这里不知道该不该要这条命令 type,3 esurf,all finish /solu csys,0 nsubst,10 nlgeom,on nropt,full pred,on lnsrch,on outres,all,last nsel,s,loc,x,length_x1 d,all,ux nsel,s,loc,x,length_x2 d,all,ux nsel,s,loc,y,width_y1 nsel,a,loc,y,width_y2 d,all,uy nsel,s,loc,z,heigth_z0+heigth_z1 d,all,uz allsel,all acel,0,0,-10 time,1 solve !土体自重应力分析，不知有没有必要 *dim,fxyz,TABLE,5,1 *CREATE,ansuitmp *vread,fXYZ(1,1),fxyz5row,txt,,JIK,1,5 (e14.6) *vread,fXYZ(1,0),time5row,txt (e14.6) fXYZ(0,1)=1 *end /INPUT,ansuitmp finish /SOLU ANTYPE,trans btime=2.683363E-03 etime=3.202147E-02 dtime=0.00733453 *DO,itime,btime,etime,dtime time,itime nsub,10 KBC,1 !esel,s,type,,3 <-加入这条命令会出现问题，不能正常求解， !不加则可进行。 sfe,all,5,pres,0,fxyz(itime),0,0,1 SOLVE *ENDDO 在post1中 Main Menu> General Postproc> Plot Results> Contour Plot> Nodal Solu查看节点的uz位移发现变形严重超出想象。 这里想请各位高人帮忙分析分析，看看我的问题到底出在哪里，是处理方法？还是数据参数设置？还是求解设置？还是什么别的原因 ，俺是彻底没招了 。。。拜谢 [ Last edited by laizuliang on 2008-3-30 at 11:20 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
