24小时热门版块排行榜

>论坛更新日志 (308)
>虫友互识 (24)
>考博 (12)
>考研 (3)
>访问学者 (1)
>休闲灌水 (1)

返回列表

当前主题已经存档。

kantman

应助: (幼儿园)
在线:
虫号: 274411

[交流] 【求助】用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 ]

回复此楼

» 猜你喜欢

1楼 2008-03-26 23:49:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 kantman 的主题更新

返回列表