求助:谱分析的位移太大

hmily2000

这是命令流，单位都是按国际单位来的，为什么算出来的位移有200呢:'( ，请高手指点一下吧:'(
/clear,start
/fileame,cassion

/prep7
et,1,plane42!划分网格
et,2,solid65!混凝土单元
et,3,solid45!岩石、回填土单元

!沉箱材料
mp,ex,1,2.9e10
mp,prxy,1,0.235
mp,dens,1,4870

!基床材料
mp,ex,2,2.8e10
mp,prxy,2,0.3
mp,dens,2,1800

!岩基材料
mp,ex,3,2.9e10
mp,prxy,3,0.3
mp,dens,3,2600

!回填土材料
mp,ex,4,10e8
mp,prxy,4,0.35
mp,dens,4,1920

k,1!k1
k,2,0,32!k2
k,3,0,32!k3
k,4,0,33!k4
k,5,0,49!k5
k,6,37,49!k6
k,7,37,33!k7
k,8,37,33!k8
k,9,35.5,32!k9
k,10,35.5,32!k10
k,11,35.5,32!k11
k,12,35.5,0!k12
k,13,40,49!k13
k,14,40,33!k14
k,15,40,33!k15
k,16,40,32!k16
k,17,40,32!k17
k,18,40,0!k18
k,19,40,49!k19
k,20,40,33!k20
k,21,52.8,49!k21
k,22,52.8,33!k22
k,23,52.8,33!k23
k,24,52.8,32!k24
k,25,52.8,32!k25
k,26,52.8,0!k26
k,27,56.8,33!k27
k,28,58.3,32!k28
k,29,58.3,32!k29
k,30,58.3,0!k30
k,31,84.8,32!k31
k,32,84.8,0!k32

!a1
a,1,2,11,12
!a2
a,3,4,7,9
!a3
a,4,5,6,7
!a4
a,7,6,13,14
!a5
a,10,8,15,16
!a6
a,12,11,17,18
!a7
a,18,17,25,26
!a8
a,16,15,23,24
!a9
a,20,19,21,22
!a10
a,24,23,27,28
!a11
a,26,25,29,30
!a12
a,30,29,31,32



!划分沉箱网格
lsel,s,,,29,31,2
lesize,all,,,20
lsel,s,,,28,30,2
lesize,all,,,20
amesh,9

!划分正下方基床
lsel,s,,,25,27,2
lesize,all,,,20
lsel,s,,,17,26,9
lesize,all,,,2
amesh,8

!划分左下方基床
lsel,s,,,16,18,2
lesize,all,,,10
lsel,s,,,15
lesize,all,,,2
amesh,5

!划分右下方基床
lsel,s,,,32,34,2
lesize,all,,,10
lsel,s,,,33
lesize,all,,,2
amesh,10

!划分左下方岩基
lsel,s,,,19,21,2
lesize,all,,,10
lsel,s,,,3,20,17
lesize,all,,,20
amesh,6

!划分正下方岩基
lsel,s,,,22,24,2
lesize,all,,,20
lsel,s,,,23
lesize,all,,,20
amesh,7

!划分右下方岩基
lsel,s,,,35,37,2
lesize,all,,,10
lsel,s,,,36
lesize,all,,,20
amesh,11

!划分左方岩基
lsel,s,,,2,4,2
lesize,all,,,20
lsel,s,,,1
lesize,all,,,20
amesh,1

!划分右方岩基
lsel,s,,,38,40,2
lesize,all,,,20
lsel,s,,,39
lesize,all,,,20
amesh,12

!划分划分回填土2
lsel,s,,,6,8,2
lesize,all,,,20
lsel,s,,,5,7,2
lesize,all,,,2
amesh,2

!划分划分回填土3
lsel,s,,,6,10,4
lesize,all,,,20
lsel,s,,,9,11,2
lesize,all,,,20
amesh,3

!划分划分回填土4
lsel,s,,,12,14,2
lesize,all,,,5
lsel,s,,,11,13,2
lesize,all,,,20
amesh,4

!拉伸沉箱
extopt,esize,10,0
type,2
mat,1
vext,9,,,,,-12

!拉伸基床和基床后填土
extopt,esize,10,0
type,3
mat,2
asel,s,,,5,8,3
asel,a,,,10
vext,all,,,,,-12

!拉伸岩基
extopt,esize,10,0
type,3
mat,3
asel,s,,,1,6,5
asel,a,,,7,11,4
asel,a,,,12
vext,all,,,,,-12


!拉伸回填土
extopt,esize,10,0
type,3
mat,4
asel,s,,,2,4,1
vext,all,,,,,-12


!定义沉箱与基床接触面
allsel
et,4,170!刚性目标面的单元类型
et,5,173!柔性接触面的单元类型
keyopt,5,9,0
keyopt,5,12,0
r,10
real,10
mp,mu,2,0.3
asel,s,,,17
cm,target,area
type,4
nsla,s,1
esurf,all
esel,s,type,,4
allsel

asel,s,,,24
cm,contact,area
type,5
nsla,s,1
esurf,all
esel,s,type,,5
allsel


!定义基床与回填土接触面
r,11
real,11
mp,mu,4,0.2
asel,s,,,20
cm,target,area
type,4
nsla,s,1
esurf,all
esel,s,type,,4
allsel

asel,s,,,64
cm,contact,area
type,5
nsla,s,1
esurf,all
esel,s,type,,5



!定义基床与岩基接触面
allsel
r,12
real,12
mp,mu,2,0.3
asel,s,,,37,41,4
asel,a,,,45
cm,target,area
type,4
nsla,s,1
esurf,all
esel,s,type,,4

allsel
asel,s,,,22,26,4
asel,a,,,30
cm,contact,area
type,5
nsla,s,1
esurf,all
esel,s,type,,5


!沉箱与回填土接触
allsel
r,13
real,13
mp,mu,4,0.3
asel,s,,,14
cm,target,area
type,4
nsla,s,1
esurf,all
esel,s,type,,4

allsel
asel,s,,,63
cm,contact,area
type,5
nsla,s,1
esurf,all
esel,s,type,,5

!回填土与岩基接触
allsel
r,14
real,14
mp,mu,4,0.2
asel,s,,,33
cm,target,area
type,4
nsla,s,1
esurf,all
esel,s,type,,4

allsel
asel,s,,,56
cm,contact,area
type,5
nsla,s,1
esurf,all
esel,s,type,,5

!施加边界条件
/sol
allsel
asel,s,loc,x,0!左侧边界
da,ux,all

asel,s,loc,x,84.8!右侧边界
da,ux,all

asel,s,loc,y,0!底边界
da,all,all

allsel
asel,s,,,16
sfa,all,,press,-1!水平面
asel,s,,,17
sfa,all,,press,1!垂直面

!施加重力加速度
allsel
acel,0,9.8,0,
allsel,all
save
finish

!模态分析
/sol
antype,modal
MODOPT,subsp,12
save
solve
finish

!模态扩展
/sol
antype,modal
expass,on
mxpand,12,,,yes
solve
finish

!谱分析
/solu
antype,spectr
spopt,psd,,yes
psdunit,1,pres,,
sed,-1,0,0
psdfrq,1,,0.03,0.04,0.05,0.07,0.09
psdfrq,1,,0.10,0.11,0.13,0.15,0.17
psdfrq,1,,0.25,0.28,,,
psdval,1,1.06e4,1.86e6,9.03e7,7.25e9,1.89e10
psdval,1,1.48e10,9.98e9,5.54e9,3.35e9,2.16e9
psdval,1,5.55e8,3.72e8,,,
asel,s,,,16
nsla,s,1
LVSCALE,1
allsel,all
pfact,1,node
psdunit,2,pres,,
sed,0,1,0
psdfrq,2,,0.03,0.04,0.05,0.07,0.09
psdfrq,2,,0.10,0.11,0.13,0.15,0.17
psdfrq,2,,0.25,0.28,,,
psdval,2,2.80e2,4.92e4,2.39e6,1.92e8,5.01e8
psdval,2,3.91e8,2.64e8,1.47e8,8.88e7,5.71e7
psdval,2,1.47e7,9.87e6,,,
asel,s,,,17
nsla,s,1
lvscale,1
allsel,all
PFACT,2,node
PSDRES,DISP,ABS
PSDRES,VELO,OFF
PSDRES,ACEL,OFF  

!模态合并  
/SOLU
ANTYPE,SPECTR
psdcom,0.005,12
solve
finish

hmily2000

没有高人指点一下吗:'( :'( :'(