ansys谱分析结果过小，请高手帮忙指点！！

qingchen99

大家好，我在用ansys分析一个拱结构，跨度60米，谱分析的结果最大应力才0.6MPa，位移0.02mm，当我查看单个振型的计算结果的时候，结果却是正常的，请高手帮忙指点！！以下是我的命令流：
finish
/clear
/filname,arch

!建模，单位为N，m
/prep7
et,1,solid65,,,,,,,1
et,2,solid45
et,3,plane42
mp,ex,1,29254e6                !混凝土弹性模量29254MPa
mp,prxy,1,0.2
mp,dens,1,2500                !混凝土密度2500Kg/m3

mp,ex,2,173900e6        !钢材弹性模量173900MPa
mp,prxy,2,0.283
mp,dens,2,7800                !钢材密度7800Kg/m3

k,1
k,2,97e-3
k,3,4e-3,4e-3
k,4,93e-3,4e-3
k,5,4e-3,207e-3
k,6,93e-3,207e-3
k,7,,211e-3
k,8,97e-3,211e-3

a,3,4,6,5
a,1,3,5,7
a,4,2,8,6
a,1,2,4,3
a,5,6,8,7

lsel,s,tan1,x
lesize,all,,,10
lsel,s,tan1,y
lesize,all,,,5
lsel,s,length,,4*sqrt(2)
lesize,all,,,1
allsel,all

aatt,,,3       
amesh,all


k,10000,0,-3130.5e-3
k,10001,97e-3,-3130.5e-3
l,10000,10001

!生成混凝土圆弧体
asel,s,,,1
type,1                !选择1号面，并将其单元类型定义为1号即混凝土
extopt,esize,150        !定义体扫掠方向的单元尺寸为150
mat,1
vrotat,1,,,,,,10000,10001,74.5        !1号面绕轴旋转生成体
vrotat,1,,,,,,10000,10001,-74.5

!生成钢管圆弧体
asel,s,,,2,5,1
type,2
extopt,esize,150
mat,2
vrotat,2,3,4,5,,,10000,10001,74.5
vrotat,2,3,4,5,,,10000,10001,-74.5

allsel,all

!施加端部约束
asel,s,,,10,25,5
asel,a,,,28,31,3
asel,a,,,36
asel,a,,,41,47,3
da,all,all
dtran                !将几何模型上的约束转换到有限元模型上

!施加平面外约束
nsel,s,loc,x,0
nsel,a,loc,x,97e-3
nsel,r,loc,z,-1540e-3,-1200e-3
d,all,ux

nsel,s,loc,x,0
nsel,a,loc,x,97e-3
nsel,r,loc,z,1200e-3,1540e-3
d,all,ux

nsel,s,loc,x,0
nsel,a,loc,x,97e-3
nsel,r,loc,z,-640e-3,-300e-3
d,all,ux

nsel,s,loc,x,0
nsel,a,loc,x,97e-3
nsel,r,loc,z,300e-3,640e-3
d,all,ux

allsel,all

nummrg,node
numcmp,node


!加速度谱计算（按规范2001的设计反应谱计算）
tg=0.65                !特征周期
refmax=0.08        !水平地震影响系数
gar=9.8                !重力加速度
kes=0.05        !阻尼比
eta1=0.02+(0.05-kes)/(4+32*kes)                !直线下降段的下降斜率调整系数计算
eta2=1+(0.05-kes)/(0.08+1.6*kes)        !阻尼调整系数计算
gama=0.9+(0.05-kes)/(0.3+6*kes)                !曲线下降段的衰减系数计算
*dim,ttt,,20                !定义存放周期的数组，取20个点描述普曲线
*dim,fre,,20                !定义存放频率的数组
*dim,ace,,20                !定义存放加速度的数组

ttt(1)=0                !定义20个点的周期值
ttt(2)=0.1
ttt(3)=0.65       
*do,i,4,6               
ttt(i)=0.8+(i-4)*0.1
*enddo
*do,i,7,16               
ttt(i)=1.2+(i-7)*0.2
*enddo
ttt(17)=3.25
*do,i,18,20               
ttt(i)=4+(i-18)
*enddo

ace(1)=0.45*refmax*gar        !计算第1点对应的加速度（0到0.1为斜直线）
ace(2)=eta2*refmax*gar        !计算第2点对应的加速度（0.1到0.65为水平直线）
ace(3)=ace(2)                !计算第3点对应的加速度
*do,i,4,17                !循环计算曲线段对应的加速度
ace(i)=(tg/ttt(i))**gama*eta2*refmax*gar
*enddo
*do,i,18,20                !循环计算斜直线段对应的加速度
tempi=eta2*0.2**gama-eta1*(ttt(i)-5*tg)
ace(i)=tempi*refmax*gar
*enddo
fre(20)=1e5                !将频率求出，并采用升序排列
*do,i,2,20                !求频率
fre(21-i)=1/ttt(i)
*enddo
*do,i,1,20                !将加速度数组倒置并临时存放在ttt数组中
ttt(21-i)=ace(i)
*enddo
*do,i,1,20                !再将ttt数组导入与频率对应的加速度
ace(i)=ttt(i)
*enddo



!获得模态解
/solu
antype,modal
modopt,lanb,20
solve
save
finish


!获得谱解
/solu
antype,spectr
spopt,sprs,20,yes        !定义谱分析类型为单点谱分析
dmprat,0.05                !定义常阻尼比
sed,0,0,1                !定义激励方向为z向，即平面内
svtyp,2                        !定义单点响应谱类型为加速度谱
*do,i,1,20                !定义普曲线
freq,fre(i)
sv,,ace(i)
*enddo
solve
save
finish


!扩展模态
/solu
antype,modal
expass,on
mxpand,20,0,0,yes,0.001        !重要性超过0.001则进行扩展
solve
save
finish



!模态合并
/solu
antype,spectr
srss,0.01,disp
solve
save
finish