24小时热门版块排行榜    

查看: 4951  |  回复: 16

zhudongshan

铁虫 (初入文坛)

[求助] 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,发现与设置一个很小时间步长的结果没有区别)。
是不是命令流有错?请教下大家,被这个问题卡了好久,谢谢~

ANSYS做地震时程分析,为什么地震作用下的位移反应要比静力自重条件下要小?
图1.png模型图


ANSYS做地震时程分析,为什么地震作用下的位移反应要比静力自重条件下要小?-1
图2.静力位移图(最大位移0.00688m).png


ANSYS做地震时程分析,为什么地震作用下的位移反应要比静力自重条件下要小?-2
图3.迁安波.png


ANSYS做地震时程分析,为什么地震作用下的位移反应要比静力自重条件下要小?-3
图4.输入迁安波(截取前3s).png


ANSYS做地震时程分析,为什么地震作用下的位移反应要比静力自重条件下要小?-4
图5.前3s模型顶部节点位移时程(设置了很小时间步长).png


ANSYS做地震时程分析,为什么地震作用下的位移反应要比静力自重条件下要小?-5
图6.前3s模型顶部节点位移时程(未设置很小时间步长).png
回复此楼
生命不止,奋斗不息。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZK2194184

铁杆木虫 (正式写手)

底部的约束改为全约束,即约束x和y两个方向;另外就是要保证你读取的地震波没问题,可以进一步交流

发自小木虫Android客户端
2楼2017-06-25 09:24:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhudongshan

铁虫 (初入文坛)

底部约束了两个方向,还是不行,查看了底部节点加速度时程,与输入地震波大致相同
生命不止,奋斗不息。
3楼2017-06-25 11:43:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZK2194184

铁杆木虫 (正式写手)

引用回帖:
3楼: Originally posted by zhudongshan at 2017-06-25 11:43:41
底部约束了两个方向,还是不行,查看了底部节点加速度时程,与输入地震波大致相同

底部约束位置加速度时程应该为0的啊,它被约束了,咋还能有位移,加速度呢

发自小木虫Android客户端
4楼2017-06-26 09:27:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhudongshan

铁虫 (初入文坛)

引用回帖:
4楼: Originally posted by ZK2194184 at 2017-06-26 09:27:52
底部约束位置加速度时程应该为0的啊,它被约束了,咋还能有位移,加速度呢
...

是,位移加速度都是零,底部往上第一个点
生命不止,奋斗不息。
5楼2017-06-26 10:45:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhudongshan

铁虫 (初入文坛)

引用回帖:
4楼: Originally posted by ZK2194184 at 2017-06-26 09:27:52
底部约束位置加速度时程应该为0的啊,它被约束了,咋还能有位移,加速度呢
...

还有就是,我怀疑这个重力是不是没加进去,因为我又算了一个不加重力的,直接加地震加速度时程,发现与在静力条件下加地震时程的结果一模一样
生命不止,奋斗不息。
6楼2017-06-26 10:47:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

金亚

新虫 (初入文坛)

你的荷载步命令流不对

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

7楼2017-06-26 15:32:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhudongshan

铁虫 (初入文坛)

送红花一朵
引用回帖:
7楼: Originally posted by 金亚 at 2017-06-26 15:32:19
你的荷载步命令流不对

你好,请问下哪不对,应该怎么写?

发自小木虫Android客户端
生命不止,奋斗不息。
8楼2017-06-26 16:00:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feifei836

至尊木虫 (著名写手)

时间步有问题

发自小木虫IOS客户端

» 本帖已获得的红花(最新10朵)

9楼2017-06-27 07:59:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhudongshan

铁虫 (初入文坛)

送红花一朵
引用回帖:
9楼: Originally posted by feifei836 at 2017-06-27 07:59:23
时间步有问题

你好,您能说的详细些吗,应该怎么改?谢谢~
生命不止,奋斗不息。
10楼2017-06-27 09:02:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhudongshan 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 5sJFDKp2kwuI 2026-11-06 6/300 2026-11-08 18:48 by l48hFd8Nwae2
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-08 18:43 by l48hFd8Nwae2
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 2rKfVQg2jQww 2026-11-06 6/300 2026-11-08 18:26 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:17 by Xc2FxHsHkBNc
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
信息提示
请填处理意见