24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhudongshan

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 427.1
红花: 2
帖子: 48
在线: 107.3小时
虫号: 3264241
注册: 2014-06-09
性别: GG
专业: 岩土与基础工程

[求助] ANSYS做地震时程分析，为什么地震作用下的位移反应要比静力自重条件下要小？

用ANSYS做隧道的地震时程分析，先做静力（即自重）条件下分析，后输入地震波分析，但是结果发现：动力条件下的位移与应力都要比静力条件下小，原因是什么？
命令流如下（为简化模型，下面这是模拟二级围岩体）：
finish
/clear
/filname,soild
/prep7
et,1,plane42
tt=0.01
N=2320 !数据记录间隔以及数据个数
mp,ex,1,2.5e10
mp,prxy,1,0.25
mp,dens,1,2600
tb,dp,1
tbdata,1,1.8e6,50 !二级围岩参数，D-P模型
k,1
k,2,150
k,3,150,120
k,4,,120
l,1,2
l,2,3
l,3,4
l,4,1
al,all
lesize,all,,,16
amesh,all
dl,1,,uy
dl,2,,ux
dl,4,,ux !模型两侧X向约束，底部Y向约束
finish

/solu
/config,nres,20000
alphad,0.815
betad,0.00049 !定义阻尼系数（取模态分析前两阶频率）
antype,trans
timint,off !首先做重力荷载作用下的响应，关闭时间积分
acel,,9.8
time,1e-6 !很小的时间步长
solve

timint,on !打开时间效应
*creat,qiananread,mac !创建宏文件，读入迁安地震波，峰值-1.58621
*dim,qianan,,2320
*vread,qianan(1,1),qianan,txt,,ijk
(g9.6)
*end
qiananread
*do,i,1,N,1
acel,,qianan(i) !地震时程计算
time,tt*i
solve
*enddo
Finish

图1.模型图（150m*120m）

图2.静力位移图（最大位移0.00688m）

图3.迁安波（峰值-1.58621m/s2出现在2.38s）

图4.输入迁安波（截取前3s）

图5.前3s模型顶部节点位移时程（设置了很小时间步长）

图6.前3s模型顶部节点位移时程（未设置很小时间步长）
模型顶部节点位移最大-0.00472474m,出现在2.79s，顶部中间附近节点位移也在-0.0047m左右，现在不明白为什么动力位移要比静力条件下要小？
是重力没有考虑进去吗？可是看论坛里和王新敏书上写可以先设置一个很小的时间积分考虑重力，在打开时间积分计算，结果还是不行（一开始是先在静力条件下算重力，然后再在瞬态下输入地震波，对比图5和图6，发现与设置一个很小时间步长的结果没有区别）。
是不是命令流有错？请教下大家，被这个问题卡了好久，谢谢~

图1.png模型图