ansys谱分析结果过小，请高手帮忙指点！！
大家好，我在用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
页:
[1]