24小时热门版块排行榜    

查看: 4394  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279分求调剂 一志愿211 +16 chaojifeixia 2026-03-19 18/900 2026-03-23 19:44 by stone_128
[考研] 一志愿陕师大生物学071000,298分,求调剂 +3 SYA! 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 328求调剂,英语六级551,有科研经历 +7 生物工程调剂 2026-03-17 12/600 2026-03-23 18:18 by YMU施老师
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +5 jiajunX 2026-03-22 5/250 2026-03-23 18:07 by YMU施老师
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +5 葵梓卫队 2026-03-18 7/350 2026-03-23 16:26 by lingjue
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-17 10/500 2026-03-23 13:05 by 醉在风里
[考研] 276求调剂。有半年电池和半年高分子实习经历 +8 材料学257求调剂 2026-03-23 9/450 2026-03-23 13:01 by ztnimte
[考研] 307求调剂 +3 余意卿 2026-03-21 3/150 2026-03-23 10:32 by Iveryant
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 0856材料专硕353求调剂 +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
信息提示
请填处理意见